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Magnetic  Resonance  Imaging  (MRI)  provides  radiologists  with  the  tools  to  better 

diagnose  disease  and  view  internal  injuries.  As  field  strengths  increase,  so  does  resolution 

and  so  do  potential  heahh  effects.  This  dissertation  focused  on  tissue  heating  generated  by 

the  radio-fi-equency  (RF)  energy.  Thermal  measurements  were  first  made  using 

thermocouple  probes.  Experimental  results,  which  were  affected  by  the  metallic  nature  of 

the  thermocouple,  were  discarded.  The  experiments  were  repeated  using  fiber  optic 

probes.  Tests  reveled  uneven  thermal  energy  deposition  throughout  the  head  phantom. 

Temperatures  nearer  the  outer  edge  increased  sooner.  Additionally,  there  was  thermal 

movement  resulting  in  increased  temperatures  to  the  center  of  the  head  phantom  up  to  an 

hour  after  the  RF  fields  were  removed.  Finally,  a  dififlision  imaging  technique  was 

performed  prior  to  and  after  a  scan.  Analysis  of  the  resultant  images  indicated  diflRision 

imaging  may  be  able  to  monitor  temperatures  less  than  the  1  °C  as  stated  in  literature. 


CHAPTER  1 
INTRODUCTION 


The  applied  theory  of  magnetic  resonance  provided  imaging  and  spectroscopy. 
Increased  field  strengths  enhanced  resolution  providing  radiologists  improved,  non- 
invasive diagnostic  tools.  Increased  field  strengths  also  brought  increased  potential  for 
health  effects.  The  MRI  system  utilizes  the  nuclear  spins  of  isotopes  to  define  the 
presence  of  the  isotope.  At  the  three  Tesla  (3T)  MRI  unit,  jointly  owned  by  the  Veterans 
Administration  and  Shands  Teaching  Hospital,  the  isotope  of  interested  during  the 
research  of  this  project  was  the  proton  (hydrogen).  Since  the  human  body  contains  a  high 
percentage  of  water,  proton  imaging  would  provide  some  type  of  image  fi"om  every  tissue 
type  and  every  location  within  the  body.  MRI  produces  three  fields;  the  static  magnetic 
field,  the  gradient  magnetic  field  and  the  RF  field. 

MRI  Fields 
Static  Magnetic  Fields 

According  the  to  Food  and  Drug  Administration  (FDA),  clinical  MRI  systems  in 
the  United  States  were  permitted  to  fiinction  on  a  routine  clinical  basis  at  static  magnetic 
field  strengths  of  up  to  2T.  This  was  recently  increased  to  4T.  Above  this  level,  evidence 
of  safety  must  be  provided  by  the  sponsor  or  device  manufacturer  prior  to  routine  clinical 
use.  Static  magnetic  fields  have  no  effect  on  skin  or  body  temperature  (Budinger,  1992). 
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Heat  was  only  deposited  in  tissue  when  current  flowed  through  the  tissue  and  current  was 
only  induced  by  a  changing  a  magnetic  field.  In  this  respect  tissue  was  invisible  to  static 
magnetic  fields.  Patient  movement  within  a  static  magnetic  field  was  considered  a 
dynamic  magnetic  field  fi'om  the  tissue's  viewpoint.  However,  this  movement  was  slow 
and  was  insufficient  in  creating  a  noticeable  current. 

Problems  that  have  been  reported  by  patients  after  entering  the  MRI  chamber  but 
before  the  gradients  fields  were  applied  were  dizziness  and  nausea.  These  efifects  mimic 
motion  sickness  and  have  a  similar  root  cause.  Balance  has  always  been  perceived  by  the 
circulating  inner  ear  fluids.  The  calcium  carbonate  particles,  called  otoliths,  moved  against 
the  hair  cells  in  the  saccule  and  utricle  sending  false  signals  to  the  brain  that  the  head  was 
moving.  The  brain  received  conflicting  signals  fi^om  the  inner  ear  and  the  other  senses. 
Increased  discomfort  compounded  the  problem.  For  individuals  easily  susceptible  to 
motion  sickness,  very  little  head  movement  was  required.  Since  this  problem  dealt  with 
the  motion  of  the  otoliths,  this  effect  was  also  associated  with  gradient  magnetic  fields. 

Gradient  Magnetic  Fields 

Effects  have  been  reported  by  patients  subjected  to  changing  magnetic  fields.  At 
switching  rates  of  two  Tesla  per  second  (T/s),  patients  have  reported  light  flashes  called 
magnetophosphenes.  Switched  magnetic  field  excited  the  optic  nerve  inducing  electrical 
currents.  The  optic  nerve  responded  as  if  it  received  a  light  pulse  fi"om  the  eye.  Since 
these  were  random  and  isolated,  the  result  appeared  as  flashes  of  light.  This  effect  was 
also  be  achieved  in  a  static  magnetic  field  by  rapid  movement  of  the  head  (Kanal  1990). 
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At  switching  rates  of  60  T/s  patients  have  experienced  neuromuscular  excitation. 
Faraday's  Law  of  Induction  stated  that  a  conductor,  in  the  presence  of  a  time- varying 
magnetic  field,  would  induce  a  voltage  in  the  conductor  that  was  oriented  perpendicular  to 
the  rate  of  change  of  the  magnetic  field.  When  a  voltage  difference  along  a  nerve 
approached  6  V/m,  neurological  impulses  trigger  a  peripheral  nerve  stimulation.  This 
peripheral  nerve  stimulation  was  an  uncontrolled,  involuntary  skeletal  muscle  contraction 
and/or  twitching  (Reilly  1989,  Schaefer  1992).  One  FDA  guideline  to  limit  patient 
exposure  to  time-varying  magnetic  fields  specifically  required  demonstration,  with  "valid 
scientific  evidence,"  that  the  rate  of  change  of  magnetic  field  for  the  system  was  not 
sufficient  to  cause  peripheral  nerve  stimulation.  At  switching  rates  of  500  T/s  cardiac 
excitation  became  possible.  A  perturbation  of  an  electrocardiogram  was  identified  as  a 
superimposed  signal  of  the  gradient  pulse  (Budinger  1992).  No  long  term  effects  were 
demonstrated  by  the  superimposed  signal. 

While  biological  tissue  was  predominantly  transparent  to  dynamic  magnetic  fields, 
there  were  still  elements  within  the  tissue  that  responded  to  magnetic  fields.  As  magnetic 
fields  increased,  these  affected  elements  began  causing  noticeable  effects.  Currents, 
induced  by  the  changing  magnetic  field,  flowed  through  the  resistive  body  fluids  to 
produce  heat.  The  amount  of  heat  or  energy  produced  was  a  function  of  the  fi-equency, 
intensity  and  duration  of  the  exposure  and  the  coupling  between  the  RF  coil  and  the 
subject. 
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Radio-Frequency  Fields 

The  predominant  efifect  of  RF  fields  was  RF  heating.  The  biological  factors  of 
importance  were  conductivity,  specific  gravity  of  the  tissue,  size  of  the  current  loop, 
circulation  and  the  ability  of  the  body  to  loose  heat.  An  issue  that  often  arose  was  not 
whether  localized  heating  existed  but  how  much  and  what  was  its  biological  impact.  The 
majority  of  damage  to  cells  and  molecules  was  through  ionization  caused  by  heat.  RF 
fields  were  too  high  in  fi^equency  to  stimulate  excitable  tissues  electrically  however  these 
fields  deposited  energy  into  tissues  which  raised  their  temperatures.  The  normal  body 
temperature  remained  approximately  37  °C.  The  highest  safe  core  temperature  for 
individuals  was  39.4  "C.  The  skin  pain  threshold  was  43  °C  and  proteins  denatured  at 
about  45  °C.  This  small  range  of  temperatures  became  critical  at  higher  field  strengths. 
The  potential  for  localized  heating  and  subsequent  tissue  damage  increased  as  RF 
fi'equencies  increased  to  the  point  where  the  patient's  size  approached  one-half 
wavelength  of  the  incident  excitation  fi^equency  (Gordon  1992). 

An  area  of  major  concern  for  MRI  of  the  head  was  the  eye.  The  vitreous  body  of 
the  eye  contained  a  colorless,  jelly  mass  between  the  lens  and  the  retina.  Heat  removal 
resulted  fi^om  the  formation  and  circulation  of  the  aqueous  humor  as  well  as  from  the 
ophthalmic  artery.  The  eye  had  a  reduced  capability  for  heat  dissipation  as  a  resuh  of  a 
lack  of  vascularization.  Acute  near-field  exposures  of  RF  radiation  to  the  eye  of 
laboratory  animals  created  cataracts  as  a  resuh  of  the  thermal  disruption  of  ocular  tissues. 
The  threshold  for  cataractogenesis  had  been  estimated  at  100  W/kg.  An  investigation 
revealed  no  discernible  effects  on  the  eyes  of  rats  resulting  from  MR  procedures  at 
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exposures  that  far  exceeded  typical  levels  used  in  the  clinical  setting  (Shellock  1992). 

Patients  exposed  to  SAR  levels  of  0.5  to  1.3  W/kg  had  shown  significant 
elevations  in  temperatures  in  the  cornea.  The  cornea  was  anterior  to  the  lens  and 
contained  no  blood  vessels.  Heatmg  of  the  cornea  had  the  potential  to  lead  to  keratitis,  an 
inflammation  of  the  cornea  while  heating  of  the  lens  accelerated  the  onset  of  cataracts. 
The  cornea  was  pain-sensitive  and  a  conscious  patient  would  be  able  to  feel  this  heating 
before  it  became  severe. 

Various  underlying  health  conditions  affected  an  individual's  ability  to  tolerate  RF 
heating.  The  state  of  a  patient's  health,  cardiovascular  conditioning  and  obesity  were 
primary  factors.  Various  medications  (e.g.  diuretics,  beta-blockers,  calcium  blockers, 
amphetamines,  muscle  relaxers  and  sedatives)  also  altered  the  thermoregulatory  responses 
to  RF  heating.  The  heat  loss  characteristics  of  the  body  were  not  applicable  to  head 
phantom  studies  so  tissue  heating  characteristics  from  MRI  to  the  head  phantom, 
therefore,  provided  a  worst-case  scenario.  Such  conditions  would  exist  for  a  stroke  victim 
with  limited  thermoregulatory  responses  in  the  head. 

Static  and  dynamic  magnetic  fields  produced  temporary  effects  to  patients  with  no 
long  term  consequences.  The  effects  fi-om  RF  fields  resulted  in  increased  temperatures  to 
patients  which  could  be  dangerous.  A  head  phantom  provided  a  versatile  platform  for 
conducting  heating  tests  and  defining  methods  to  measure  and  monitor  thermal  changes. 
In  order  to  determine  what  method  to  use  to  measure  temperature  changes,  the  range  of 
acceptable  temperature  changes  were  derived  fi-om  standards. 


Current  Standards 
Radiation  falls  into  one  of  two  categories:  ionizing  and  non-ionizing.  Ionizing 
radiation  is  defined  as  having  suflBcient  energy  to  ionize  an  atom.  MRI  falls  into  the 
category  of  non-ionizing  radiation  since  its  electromagnetic  fields  create  insufiBcient 
energy  to  ionize  molecules,  but  sufficient  energy  for  thermalization.  To  quantify  this 
thermal  action  the  energy  absorption  over  a  small  mass  of  material  is  defined  as  the 
specific  absorption  rate  (SAR)  and  is  given  in  units  of  watts  per  kilogram.  Although  the 
SAR  defines  the  energy  absorbed  within  tissue,  it  has  always  been  impractical  to  measure 
this  energy  change  in  live  subjects.  Instead  the  surface  power  density,  given  in  W/cm^,  has 
been  measured.  Utilizing  the  tissue  density  and  tissue  attenuation  characteristics,  the  SAR 
can  be  inferred.  The  SAR  has  been  the  currently  accepted  means  of  indicating  RF 
exposure  levels  and  several  organizations  have  established  guidelines  for  permissible  levels 
of  exposure  to  static  and  varying  magnetic  fields  and  to  RF  fields. 

INIRC-IRPA 

The  International  Non-Ionizing  Radiation  Committee  of  the  International 
Radiation  Protection  Association  (INIRC-IRPA)  is  an  international  association  providing 
universal  guidance  on  a  variety  of  non-ionizing  radiation  standards  (ESflRC-IRPA  1991). 
For  static  magnetic  fields,  a  patient's  cardiovascular  fixnction  should  be  monitored  during 
exposures  above  2T  to  the  head  or  trunk  and  above  5T  for  limbs.  For  time-varying 
magnetic  fields,  a  patient's  cardiovascular  function  should  be  monitored  during  exposures 
above  6  T/s.  Tune-varying  magnetic  fields  should  not  exceed  20  T/s.  For  RF  fields. 
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exposures  up  to  one  hour  should  not  exceed  a  total  energy  deposition  of  120  W-min/kg 
which  equates  to  a  SAR  of  2  W/kg.  To  avoid  overheating,  the  following  product  of  time 
and  SAR  should  not  exceed: 

60  W-min/kg  (1  W/kg  for  1  hour)  averaged  over  the  head, 

120  W-min/kg  (2  W/kg  for  1  hour)  averaged  over  the  trunk  or 

180  W-min/kg  (3  W/kg  over  1  hour)  averaged  over  the  extremities. 
Instantaneous  SAR  should  not  exceed: 

4  W/kg  in  the  head, 

8  W/kg  in  the  trunk, 

12  W/kg  in  the  extremities  or 

10  W/kg  in  0.01  kg  or  more  than  10  minutes  for  the  eyes. 

The  National  Council  on  Radiation  Protection  and  Measurements  (NCRP) 
Scientific  Committee  39,  Report  No.  67,  is  entitled  Radiofrequency  Electromagnetic 
Fields:  Properties,  Quantities  and  Units,  Biophysical  Interaction  and  Measurements 
(NCRP  67).  Report  No.  67  is  a  primary  source  upon  which  American  National  Standards 
Institute  (ANSI)  C95. 1-1982  is  based  and  it  concluded  there  is  little  possibility  of  directly 
measuring  the  absorption  of  energy  by  biological  bodies  at  the  cellular  level.  This 
becomes  important  when  dealing  with  hot  spot  generation  which  is  the  concentration  of 
energy  into  a  small  volume.  Cell  death,  especially  in  the  brain  area,  can  cause  irreversible 
damage. 
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The  3T  MRI  system,  for  proton  imaging,  operates  at  128  MHz  limiting  the 
equivalent  power  density  to  one  mW/cm^  (Table  1-1).  These  guidelines  make  no 
provisions  for  near-field  radiation.  In  far-field  radiation,  generally  one  wavelength  ft^om 
the  antenna,  the  electric  and  magnetic  fields  become  well  behaved  and  remain  orthogonal 
to  each  other  and  the  direction  of  propagation.  This  permits  power  density  calculations 
using  either  electric  or  magnetic  field  measurements.  In  the  near-field,  close  to  the 
antenna,  the  electric  field  and  magnetic  fields  have  not  aligned  with  each  other. 
Additionally,  there  is  a  reactive  field.  This  field  is  the  inductive  or  capacitive  storage  of 
energy  by  the  antenna  itself  It  does  not  propagate  far  fi-om  the  anterma  but  still  can 
interact  with  tissue  within  its  influence.  ANSI  C95. 1-1982  provides  results  fi'om  either  the 
electric  field  or  the  magnetic  field  and  uses  far-field  assumptions  to  calculate  SAR.  The 
3T  MRI  system  operates  only  in  the  near  field.  The  equivalent  power  density  must  then 
become  a  function  of  both  the  electric  and  magnetic  fields. 

ffiC 

The  International  Electrotechnical  Commission  (lEC)  provides  extensive  standards 
for  RF  exposures.  A  complete  listing  of  the  lEC  standards  are  provided  m  Appendix  A. 
The  following  requirements  are  applicable  for  room  temperatures  below  24  °C  and  where 
the  relative  humidity  does  not  exceed  60%.  The  requirements  regarding  SAR  fi'om  RF 
power  fijrther  assumes  that  the  SAR  averaged  over  any  10  second  period  does  not  exceed 
5  times  the  stated  temporal  average  SAR  limit. 

a.         Whole-body  SAR.  Three  operating  modes  with  regard  to  whole-body 
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SAR  are  defined,  namely  Normal  Operating  Mode,  First  Level  Controlled 
Operating  Mode  and  Second  Level  Controlled  Operating  Mode 
(Figure  1-1). 
i.  The  Normal  Operating  Mode  comprises  values  of  whole-body  SAR 

not  higher  than  1.5  W/kg  averaged  over  any  period  of  15  minutes. 
ii.         The  First  Level  Controlled  Operating  Mode  comprises  values  of 

whole-body  SAR  not  higher  than  4  W/kg  averaged  over  any  period 

of  1 5  minutes. 
iii.        The  Second  Level  Controlled  Operating  Mode  comprises  values  of 

whole-body  SAR  that  may  exceed  4  W/kg  averaged  over  any 

period  of  1 5  min. 
Head  SAR.  Two  operating  modes  with  regard  to  head  SAR  are  defined, 
namely  Normal  Operating  Mode  and  Second  Level  Controlled  Operating 
Mode  (Figure  1-2). 
i.  The  Normal  Operating  Mode  comprises  values  of  head  SAR  not 

higher  than  3  W/kg  averaged  over  any  period  of  10  minutes. 
ii.         The  Second  Level  Controlled  Operating  Mode  comprises  values  of 

head  SAR  that  may  exceed  3  W/kg  averaged  over  any  period  of  10 

minutes. 
Local  Tissue  SAR  With  the  Use  of  Special  Coils.  Two  operating  modes 
with  regard  to  local  tissue  SAR  are  defined,  namely  Normal  Operating 
Mode  and  Second  Level  Controlled  Operating  Mode. 
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The  Normal  Operating  Mode  comprises  values  of  local  tissue  SAR 
in  any  one  gram  of  tissue  not  exceeding  8  W/kg  in  the  head  or  torso 
or  12  W/kg  in  the  extremities  averaged  over  any  period  of  5  min. 
The  Second  Level  Controlled  Operating  Mode  comprises  values  of 
local  tissue  SAR  that  may  exceed  the  upper  limit  for  the  Normal 
Operating  Mode. 


FDA 


The  Food  and  Drug  Administration  (FDA)  provides  the  approval  for  MRI  systems 
in  the  U.S.  Currently,  static  magnetic  fields  up  to  2T  are  permitted  for  clinical  magnetic 
resonance  systems.  Above  that  level,  the  system  must  show  evidence  of  safety  before 
routine  clinical  use.  Therefore  the  3T  system  has  not  been  approved  by  the  FDA  for 
patient  use  but  has  been  permitted  as  a  research  unit.  The  FDA  provided  three 
alternatives  for  patient  exposure  to  time- varying  magnetic  fields. 

1.  Demonstrate  that  the  maximum  change  in  magnetic  field  strength  (dB/dt) 
of  the  system  is  6  T/s  or  less.  This  is  1/10  the  threshold  of  peripheral  nerve 
stimulation  and  almost  1/100  the  possible  threshold  for  cardiac  excitation. 

2.  A.        Demonstratethat  for  axial  gradients  (Figure  1-3): 

L  dB/dt  <  20  T/s  for  T  >  1 20  us; 

u.         dB/dt  <2400/(t(^s))  T/s  for  12  ^s<T<  120  ns  or 

iii.        dB/dt  <  200  T/s  for  T<  12  jis. 
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B.        Demonstrate  that  for  transverse  gradients,  dB/dt  is  considered  to  be 
below  the  level  of  concern  when  less  than  three  times  the  above 
limits  for  axial  gradients. 
3.         Demonstrate  with  "valid  scientific  evidence"  that  the  rate  of  change  of 
magnetic  fields  for  the  system  are  not  sufficient  to  cause  peripheral  nerve 
stimulation  by  an  adequate  margin  of  safety  (which  the  FDA  establishes  as 
at  least  a  factor  of  three). 
Exposure  to  RP  energy,  which  the  FDA  considers  below  the  level  of  concern,  is  a 
SARof 

1.  0.4  W/kg  or  less  averaged  over  the  body; 

2.  8.0  W/kg  or  less  spatial  peak  in  any  1  gram  of  tissue,  or 

3.  3.2  W/kg  or  less  averaged  over  the  head. 

The  FDA  has  stated  that  the  exposure  to  RF  energy  must  not  increase  the  body's 
core  temperature  by  1  °C  and  that  the  temperature  not  exceed  39  °C  in  the  trunk  and  40  "C 
in  the  extremities.  Normal  body  temperatures  range  fi-om  35.5  °C  to  40  °C  with  an 
averagejtjf^^.  ^^^^^  ^  document  on  September  29,  1997  (Appendix  B)  titled 
"Guidance  for  Magnetic  Resonance  Diagnostic  Devices  -  Criteria  for  Significant  Risk 
Investigations."  This  document  was  designed  to  provide  guidance  only  and  did  not 
constitute  a  regulatory  change.  This  guidance  allowed  exposures  to  static  magnetic  fields 
up  to  four  Tesla.  Exposures  to  RF  energy  were  also  modified  to  SARs  of 
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1.  4  W/kg  averaged  over  the  whole  body  for  any  15  minute  period; 

2.  3  W/kg  averaged  over  the  head  for  any  10  minute  period;  or 

3.  8  W/kg  in  any  gram  of  tissue  in  the  head  or  torso  or  1 2  W/kg  in  any  gram 
of  tissue  in  the  extremities  for  any  5  minute  period. 

For  this  research,  which  deah  with  the  head,  the  SAR  limits  were  lowered  from  3.2  W/kg 
to  3  W/kg.  The  time  limits,  which  were  not  previously  a  part  of  the  guidance,  limits 
applications  of  brief  periods  of  higher  energy  depositions.  Quantitative  standards  for 
gradient  magnetic  fields  was  removed  and  replaced  with  the  non-definitive  provision  that 
gradient  fields  should  not  "produce  severe  discomfort  or  painful  nerve  stimulation." 

Contrast 

The  ANSI  standards  only  address  SAR  in  the  far-field  region  where  the  electric 
and  magnetic  fields  have  become  aligned  with  each  other.  These  standards  provide  no  real 
guidance  for  the  MRI  system.  The  INIRC-IRPA,  FDA  and  lEC  standards  address  time- 
varying  magnet  fields  with  the  predominant  reaction  of  concern  being  peripheral  nerve 
stimulation  which  can  be  most  damaging  when  affecting  the  heart  muscle.  Each  standard 
provides  a  safety  factor  of  10  with  the  new  FDA  guideline  removing  quantitative  values. 
Most  peripheral  nerve  stimulation  was  not  observed  or  anticipated  until  changing  magnetic 
fields  exceeded  60  T/s.  The  3T  MRI  can  switch  magnetic  fields  at  this  rate.  The  lEC 
standards  addressed  a  uniform  method  of  measuring  power  deposition.  Two  methods  of 
SAR  measurement  are  available.  One  method  is  designed  specifically  for  whole  body 
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SAR.  The  other  method  measures  the  average  power  absorbed  by  an  object  of  known 
mass  to  provide  a  representative  indicator  of  SAR. 

Until  ratification  of  the  new  FDA  guidance,  the  FDA  standards  must  still  be  met 
before  patient  use  of  the  3T  system.  There  was  no  guidance  on  any  uniform  testing 
procedure  available.  It  becomes  the  responsibility  of  the  particular  facility  to  provide 
proof  of  compliance.  Therefore  the  results  obtained  during  this  research  may  be  used  to 
provide  proof  of  the  safety  of  the  3T  MRI  system. 

Specific  Absorption  Rate 
Of  the  three  parameters  of  concern  during  MRI;  static  magnetic  fields,  gradient 
magnetic  fields  and  RF  fields,  the  RF  fields  provide  effects  resulting  in  possible  permanent 
harm  to  patients.  By  definition,  the  specific  absorption  rate  is  "the  time  derivative  of  the 
incremental  energy  (dW)  absorbed  by,  or  dissipated  in  an  incremental  mass  (dm)  contained 
in  a  volume  element  (dV)  of  a  given  density."  (NCRP  67).  Mathematically  this  is 
represented  by 


SAR    =    — 
dt 


i  dW] 
dm] 


d 
dt 

\    dW\ 
[pidV)] 

Realistically  the  amount  of  energy  deposited  within  a  unit  of  tissue  can  not  be 
measured.  However,  other  characteristics  can  be  measured.  Power  deposition  can  be 
inferred  if  the  amount  of  energy  provided  to  the  target  tissue  were  known  and  certain 
characteristics  of  the  target  tissue  were  also  known.  With  the  power  deposition  defined 
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and  tissue  mass  measured,  the  SAR  can  be  defined  over  any  period  of  time.  This  section 
describes  three  methods  for  the  development  of  SAR  equations:  (1)  using  loading  factors, 
(2)  using  Poynting  vectors  with  far-field  assumptions  and  (3)  using  Poynting  vectors 
without  far-field  assumptions. 

Loading  Factors 

Loading  factors  provided  RF  deposition  measurements  by  defining  the  loading 
factor  of  an  empty  coil  (Q^)  and  the  loading  factor  of  a  coil  loaded  with  the  patient  (Qj) 
(Werner  1993).  The  patient's  loading  factor,  Q^,  was  then  calculated. 

_L  =  _L  -  J_ 

Qt  ~   Qx~  Qo 

The  RF  power  P^  deposited  m  tissue  in  the  presence  of  a  linearly  polarized  B,  field  became 


^  2 


'V 


S) 


(jiL 


where  o)  -  27T;f  was  the  angular  fi-equency,  L  was  the  coil  inductance  and  S  was  a  coil 
constant  giving  the  magnetic  field  produced  per  unit  current  in  the  coil.  The  required 
value  of  Bj  depended  on  the  desired  flip  angle  a  and  the  RF  pulse  length  x  (Roschmann 
1987).  If  the  RF  pulses  had  a  rectangular  envelope,  the  required  linearly  polarized  B,  was 
to 
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B^    =    2a/Y,T 

where  Yi  was  the  gyromagnetic  ratio  of  the  considered  nucleus  i  (for  protons  Yi  =  42.58 

MHz/T). 

The  SAR  was  the  resuh  of  the  scan  time  and  mass  of  the  target. 


SAR    =    -^— 


For  the  phantom  this  was  a  valid  assumption.  For  the  human  head  this  was  not 
vaUd  but  could  still  be  safely  assumed  at  lower  field  strength  MRI  system.  The  amount  of 
tissue  heating  becomes  unique  to  each  part  of  the  brain  as  a  function  of  its  electrical 
characteristics.  Some  portions  absorb  more  energy  and  would  display  increased 
temperatures.  Overall,  temperature  increases  would  remain  reasonable  for  lower  field 
strength  MRIs  and  would  not  have  caused  any  permanent  damage.  However,  at  3T, 
power  deposition  per  unit  mass  encroached  on  the  SAR  limits.  At  the  higher  energies,  the 
generation  of  hot  spots  became  an  issue.  The  tissue  material  with  high  conductive 
properties  would  absorb  more  RF  energy  while  those  with  lower  conductive  properties 
would  absorb  less.  The  total  power  deposition  per  unit  mass,  averaged  over  the  entire 
head,  would  meet  FDA  guidelines  but  provide  no  method  to  determine  hot  spots.  Using 
loading  factors  was  not  an  acceptable  method  for  safeguarding  patient  injuries  but  may, 
however,  still  provide  a  first  step  in  determming  the  range  of  energy  deposition. 


"i.... 


■'.^ 
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SAR  Development  Avith  Far-Field  Assumptions 

Another  method  described  for  calculating  SAR  was  the  Poynting  vector.  The 
Poynting  vector  theorem  was  developed  to  find  the  power  in  a  uniform  plane  wave  made 
up  of  an  electromagnetic  field  .  An  electromagnetic  field  consists  of  two  parts,  the  electric 
field  and  the  magnetic  field.  The  ANSI  C95. 1-1982  protection  guides  provide  an 
equivalent  power  density  equation  based  on  one  of  the  measurable  fields. 


SAR    =    —  |£.P 
2p      ' 


where  o  =  tissue  conductivity 
p  =  tissue  density 
Ej  =  electric  field 
The  Poynting  vector  was  defined  as  the  cross  product  between  the  electric  and 
magnetic  fields  and  interpreted  as  an  instantaneous  power  density  with  units  of  watts  per 
area  (W/m^). 


E\- 


-       =    E  X  H 


X  H 


\mj 


2 


m 


Under  the  conditions  of  the  SAR  equations  defined  by  NCRP  Report  67,  the  target  tissue 
was  assumed  in  the  far-field.  At  the  antenna,  both  the  electric  field  and  the  magnetic  field 
exist  independently.  Both  begin  propagating  outward.  Since  there  was  no  immediate 
relationship  between  the  two  fields,  their  phase  angles  were  random.  In  the  far  field,  with 
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a  predictable  90°  relationship,  an  intrinsic  impedance,  r|,  was  defined  as  the  square  root  of 
the  ratio  of  permeability  to  permittivity. 


T1    = 


permeability 
^  permittivity  \ 


8 


Additionally,  and  simply,  in  the  far-field,  the  intrinsic  impedance,  as  a  perfect  dielectric, 
was  defined  as  the  ratio  of  the  electric  field  to  the  magnetic  field 


E 


Using  this  identity,  and  averaging  the  Poynting  vector,  the  magnetic  field  component  of 
the  Poynting  vector  was  removed. 


P  -  lEl 

2    Tl 


The  human  body  can  be  considered  a  lossy  dielectric  in  which  case  the  intrinsic 
impedance  became  a  complex  value.  The  complex  value  was  written  in  polar  form  defined 
by  a  vector  length  and  a  vector  angle  representing  the  phase  difference  of  the  electric  and 
magnetic  fields.  As  will  be  demonstrated  later,  certain  assumptions  were  made  to  reduce 
the  complex  nature  of  the  intrinsic  impedance. 


n    " 


^ 


a  +  y  (oe 


=  n„/e. 
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Additionally,  as  the  electric  field  progressed  into  the  tissue,  it  attenuated  as  a  function  of 
the  tissue's  skin  depth,  e"**,  where  d  denoted  distance  into  the  tissue  and  6  denoted  skin 
depth.  The  human  body,  and  the  head  in  particular,  have  different  tissue  types  at  different 
depths  and  the  skin  depths  of  these  tissues  vary.  For  the  development  of  the  SAR 
equation  for  the  phantom  an  assumption  of  a  homogeneous  tissue  mass  was  used. 
The  Poynting  vector  was  modified  to 


"     2        Ti  "     2    n   /e        '     2       Ti„      ^ 


The  average  power  from  a  continuously  emitting  source  applied  to  a  tissue  was  derived 
from  the  total  power  and  the  total  area  of  concern. 


average  power    =    P  X  area    = /  -  0  xA 

2        Ti„ 


The  polar  form  for  average  power  implied  a  real  and  imaginary  portion  of  power.  It  was 
only  the  real  portion,  that  defined  the  reactive  aspect  of  power,  that  lead  to  energy 
deposition  in  tissue.  The  Poynting  vector  reduced  to 


1  ^  -7dlft  n  J 

average  power e        coso    ^A 

2  "^m 


To  further  simplify  the  intrinsic  impedance,  NCRP  Report  67  assumed  the  target 
tissue  met  the  definition  of  a  good  conductor.  The  larger  the  ratio  of  conduction  current. 
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o,  to  displacement  current,  we,  the  better  the  conductor.  In  a  good  conductor,  the  electric 
and  magnetic  fields  remain  close  to  90°  out  of  phase  which  then  supports  the  previous 
assumption  used  to  remove  the  magnetic  field  component  fi"om  the  SAR  equation.  This 
was  validated  by  calculating  the  phase  relation  when  o  »  oe  for  a  given  tissue. 


tan'    _^      .    90" 
I  we  J 


Tissue  high  in  water  content  met  this  requirement.  This  included  muscles  and  both  white 
and  gray  brain  matter  and  the  tissue  equivalent  material  used  in  the  head  phantom.  For 
adipose  laden  tissue  and  bone,  with  low  water  content,  the  tissue  conductivity  and 
displacement  current  are  closer.  However,  for  the  development  of  the  SAR  equation  o  » 
coe  and  the  following  modification  was  appUed  to  the  intrinsic  impedance. 


n 


ywji 


^  o  +  70)6     ^ 


o 


ab 


Z45!    = 


4cos(tan-'(l)) 
oo 


NOTE: 


b 


This  was  derived  firom  the  propagation  constant  equation 


Y    =   7wV^ie 


N 


H8 


For  a  good  conductor,  a/c^x  »  1  and 


Y    =   yWiie. 


N 


.  a 


-J —  =  j\Fj^ 

(0|A 


=    l/PO^^coJio/IsO^ 
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-    /«MoZ45^ 


=    /(ojlo  cos45°  +yYo|iosm45° 


=    — /w^io  +j — y^ojio 

v/2  ^2 


The  real  portion  of  this  equation  was  the  attenuation  constant  which  defined  the  amount  of 
attenuation  per  distance  traveled  in  the  medium.  This  distance  was  based  on  1/e 
attenuation  which,  by  definition,  also  became  the  foundation  of  skin  depth.  The 
attenuation  constant,  therefore,  was  the  inverse  of  the  skin  depth,  6. 


^ 


Using  only  the  real  portion  of  the  intrinsic  impedance 
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ab  ao 


The  simplified  intrinsic  impedance  was  applied  to  the  average  power  equation 


average  power    =    — E^e  ""°  x  A 
The  SAR  was  defined  as  the  average  power  per  mass  of  the  target  tissue. 


cAu    -     average  power 
mass 


Substituting  the  definition  of  the  average  power 


SAR    =    — 

p  X  volume 


To  simplify  the  equation,  the  volume  of  interest  was  defined  as  an  area.  A,  down 
to  one  skin  depth,  6.  Additionally,  the  power  calculation  defined  the  incident  power  at  the 
skin  surface  (d  =  0). 


SAR    -    -^E^  . 
2p 


SAR  Development  Without  Far-Field  Assumptions 

Two  assumptions  were  made  in  order  to  derive  the  SAR  equation.  First,  that  the 
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target  tissue  resided  in  the  far-field  and  second,  that  the  target  tissue  met  the  requirement 
of  a  good  conductor.  These  assumptions  are  reviewed  in  the  next  two  sections. 

Far-Field  Assumption.  The  far-field  region  begins  at  some  finite  distance  fi'om  the  antenna 
and  stretches  to  infinity.  In  the  far-field,  the  electric  field  strength,  magnetic  field  strength 
and  power  density  are  orthogonal  to  each  other  and  one  parameter  can  define  the  others. 

P    =    E  ^  H    =    -^—    =    nOiiH^ 
120  7U 

(where  P  was  expressed  in  W/m^,  E  in  V/m  and  H  in  A/m) 
The  area  from  the  antenna  to  the  beginning  of  the  far-field  defines  the  near-field 
which  consists  of  the  radiating  field  in  addition  to  a  reactive  (nonradiating)  field.  This 
energy  can  be  capacitive  or  inductive  and  can  affect  the  electric  and  magnetic  field  vectors. 
The  reactive  region  remains  adjacent  to  the  antenna  and  extends  to  about  one  wavelength 
from  the  antenna.  Like  the  far-field,  the  electric  and  magnetic  field  intensity  decreases 
with  distance  but  unlike  the  far-field  the  radiation  still  lacks  the  plane-wave  characteristics 
of  the  far-field  electric  and  magnetic  fields.  No  relative  assumptions  are  available  to 
characterize  one  field  by  the  other.  Measurement  of  both  fields  are  required  to  accurately 
define  the  SAR.  For  a  3T  system  operating  near  128  MHz,  one  wavelength  extends  to 
2.34  m  which  placed  the  target  tissue  well  in  the  near-field  region.  Far-field  assumptions 
for  SAR  were  invalid  but  unfortunately  "...recommendations  of  RF  exposure  limits  are 
based  on  the  experimental  data  for  the  far-field  exposures  and  are  usually  formulated  in 
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terms  of  the  far-field  parameters.  Only  general  guidance  regarding  the  SAR  was  given  for 
the  near-field  exposures,"  (Stuchley  1985). 

Good  Conductor.  The  second  assumption  made  in  the  SAR  development  was  that  the 
material  imaged  with  the  MRI  was  a  good  conductor.  The  intrinsic  impedance,  r\,  was 
simplified  under  the  assumption  that  o  »  toe,  or  that  the  ratio  of  conduction  current  to 
displacement  current  (o/ooe)  was  much  greater  than  unity. 

For  the  human  head,  o/we  for  bone  and  fat  were  0.47  and  0.52  respectively,  which 
do  not  approach  the  definition  for  a  good  conductor.  Since  these  were  the  incident  tissues 
of  the  head,  the  good  conductor  assumption  for  intrinsic  impedance  became  invalid. 
Likewise,  for  the  head  phantom,  the  polyethylene  casing  had  a  low  water  content  and  had 
a  conduction  current  that  was  smaller  than  the  displacement  current,  o/o)e  =  0.0002  .  The 
electric  and  magnetic  fields,  even  if  they  had  been  orthogonal  would  have  lost  cohesion 
after  transcending  the  bone  and  fat  layers  of  the  head  or  the  polyethylene  casing  of  the 
head  phantom.  Brain  matter,  having  a  higher  water  content  than  fat  proved  to  be  a  better 
conductor  and  fell  under  the  guidelines  of  a  conduction  current  greater  than  the 
displacement  current,  o/coe  =  2.71 .  The  tissue  equivalent  material  was  constructed  mainly 
of  water  which  had  a  conduction  current  that  was  four  times  the  displacement  current, 

O/(08  =  4. 

With  the  good  conductor  assumption,  the  skin  depth  was  defined  as 


1     _         (ji[La 
6     ~    \  ~T 
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or 


To  accurately  define  the  skin  depth,  regardless  of  the  tissue's  conduction  property  will 
require  recalculation  of  the  propagation  constant. 

The  propagation  constant  combined  an  attenuation  constant  with  a  phase  constant. 
The  attenuation  constant  defined  the  loss  of  signal  strength  as  a  wave  propagated  through 
a  medium.  The  phase  constant  defined  the  angular  shift  of  the  wave  per  distance  traveled. 
The  inverse  of  the  attenuation  constant  yielded  the  skin  depth.  Mathematically,  the 
equation  for  the  propagation  constant  was  y  =  a  +  jP  and  the  equation  for  skin  depth  was 
6  =  1/a. 

Differences  Between  Near-Field  and  Far-Field  Assumptions.  Power  deposition  estimation 
from  a  theoretical  basis  in  the  near-field,  as  required  by  the  MRI,  is  complex.  No 
correlation  between  the  electric  and  magnetic  fields  exists  allowing  only  assumptions  to  be 
made.  The  assumptions  required  to  formulate  a  generalized  estimate  eliminates  the  minor 
imperfections  that  give  rise  to  individual  power  deposition  values  for  each  individual 
imaged  (Mahony  1995,  Bottomley  1985,  Bottomley  1981,  Boesiger  1992). 

The  skin  depth  was  used  for  comparison  between  near-field  and  far-field 
assumptions.  Using  far-field  assumptions. 
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a    =    —  8„E  0)    =    0.00804 

0)8     ^    " 


6    =  ^  -    0.5  m 

^Jln  X  128  X  10^  X  47r  X  lo''  x  8.04  x  lo' 

Next  the  propagation  constant,  which  contained  no  far-field  assumptions,  was  used  to 
calculate  skin  depth. 


ywy^ioE, 


1-A 


Y     =     127.5  /89.97° 

6    =    127.5  cos(89.97'')    =    0.248  w 

The  far-field  assumption  tended  to  overestimate  the  skin  depth  of  the  polyethylene 
outer  husk  by  a  factor  of  two.  Therefore,  using  far-field  assumptions  would  cause  an 
under-estimation  of  power  deposition  in  "poor"  conductors.  For  the  head  phantom,  this 
would  include  the  out  husk.  For  the  human  head  this  would  include  adipose  tissue  and 
bone.  Since  both  of  these  are  near  the  outer  surfaces,  the  under-estimation  of  power 
deposition  would  be  similar,  still  allowing  use  of  the  head  phantom  in  human  head  power 
deposition  studies. 
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Figure  1-1.  EC  Whole-Body  SAR  Guidelines. 
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Figure  1-2.  EC  Head  SAR  Guidelines. 
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Figure  1-3.  FDA  Axial 
Gradient  Magnetic  Field 
Guidelines. 
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Table  1  - 1 .  ANSI  C95 .1-1982  Protection  Guide  for  Radio-Frequency  Electromagnetic 
Radiation. 


1      Frequency 
1          range 

Equivalent  power 
density^'' 

(Electrical  field)^ 

(Magnetic  field)^ 

MHZ 

mW/cm^ 

Wm^ 

AVm^ 

0.3-3 

100 

4X10' 

2.5 

3.0-30 

900/f 

4X10^  (900/f^) 

0.025  (900/f^) 

30  -  300 

1.0 

4X10^ 

0.025 

300-1500 

fi'300 

4X10'  (f7300) 

0.025  (fi'300) 

1500-100,000 

5.0 

2X10* 

0.125 

*  measured  5  cm  or  greater  from  any  object  in  the  field  and  averaged  for  any  0. 1  h  (6  min.). 
''  (electric  field)Vl200ii  or  12Ti(niagnetic  field)^  whichever  is  greater. 


CHAPTER  2 
MATERIALS  AND  METHODS 


Experiments  were  performed  using  invasive  and  non-invasive  temperature 
measurements  at  the  VA/Shands  3T  MRI  Unit.  Invasive  measurements  were  made  using 
thermocouples  and  fiber  optic  cables.  A  thermocouple  system  was  inexpensive  but  greatly 
affected  the  results.  A  fiber  optic  system  was  expensive  but  immune  to  RF  fields  which 
provided  better  results.  Temperatures  were  measured  non-invasively  using  a  diffusion 
MRI  technique.  Ail  measurements  were  performed  on  a  head  phantom  installed  in  the 
General  Electric  (GE)  Head  Coil. 

Head  Phantom 

The  head  phantom  used  through  these  experiments  was  created  by  Charles 
Webster  and  Jim  Scott  using  a  plaster  cast  of  a  volunteer.  A  polyethylene  covering  was 
form-fitted  over  the  plaster  cast  which  was  then  broken  and  removed.  The  outer  husk  was 
constructed  of  6.35  mm  (1/4")  polyethylene  which  had  known  electrical  properties  of 
permittivity,  intrinsic  impedance  and  skin  depth  (Table  2-1).  The  polyethylene  material 
contained  a  low  water  content  and  was  suspected  of  being  a  poor  conductor. 

The  polyethylene  mold  was  filled  with  tissue  equivalent  material  (Table  2-2).  The 
tissue  equivalent  material  was  formulated  by  Jim  Scott  to  provide  similar  loading 
characteristics  in  the  GE  Head  Coil  as  the  human  head  (Olsen  1992,  Chou  1984). 
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This  material  had  a  density  of  1 .01  g/cm^  and  a  total  volume  of  5250  ml. 

Equipment  was  not  available  to  determine  the  permittivity  and  conductivity  of  the  tissue 
equivalent  material  but  certain  assumptions  allowed  useable  test  resuhs.  Tissue,  such  as 
brain  matter,  has  a  specific  heat  of  3.5  J/cm^C.  The  tissue  equivalent  material  has  a 
0.00461  to  1  ratio  of  NaCl  to  water.  Using  this  value  with  known  values  of  specific  heat, 
the  tissue  equivalent  material  specific  heat  was  extrapolated  (Table  2-3)(Leonard  1984). 
Therefore,  the  tissue  equivalent  material  should  have  thermal  characteristics  close  to  that 
of  the  brain  (Table  2-4). 

ThermQcoupk  Measuremfflts 

Temperature  measurements  m  the  tissue  equivalent  material  were  made  using  a 
thermocouple  connected  to  an  Atkins  Model  #39641-T.  The  T-Type  indicated  the 
thermocouple  was  made  of  copper-constantan.  The  junction  of  the  copper  and  constantan 
was  exposed  allowing  for  fast  response  times.  The  Atkins  unit  provided  temperatures  to 
within  0.  rC  or  0.  IT.  Since  0.  IT  was  a  smaller  unit,  a  slight  increase  in  accuracy  was 
obtained  fi^om  measurements  made  in  Fahrenheit.  Temperature  measurements  were  later 
converted  to  Centigrade. 

A  6.35  mm  (1/4")  hole  had  been  drilled  into  the  location  of  the  right  eye  of  the 
head  phantom  and  used  as  the  insertion  point  for  temperature  reading  devices.  Originally, 
the  thermocouple  was  inserted  into  the  phantom  and  left  in  the  phantom  during  the  scan. 
Interference  fi-om  the  RF  signals  distorted  the  readings  on  the  Atkins  unit  so  a  low-pass 
filter  was  installed.  The  calibration  of  the  Atkins  unit  was  invalidated  with  the  addition  of 
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the  low-pass  filter.  The  unit  was  then  disconnected  fi^om  the  thermocouple  during  the 
scan.  A  pre-existing  phantom  protocol  had  been  programmed  into  the  Advanced  NMR 
(ANMR)  console.  This  protocol  was  modified  to  provide  a  higher  range  of  output  power. 
The  modified  protocol  is  listed  in  Appendix  C. 

A  48  minute  scan  at  18  watts  average  power  (as  indicated  on  the  Omega  process 
indicator)  exhibited  a  3.9  °C  rise  in  temperature  using  the  Atkins  thermocouple.  The 
temperature  immediately  dropped  an  average  of  1. 1  °C  within  two  minutes  of  the  end  of 
the  scan.  Further  testing  revealed  that  the  temperature  of  the  tissue  equivalent  material 
was  raised  by  the  presence  of  the  thermocouple  cable  during  the  scan.  This  invalidated 
results. 

The  experimental  method  was  adjusted.  Prior  to  the  scan  the  thermocouple  was 
inserted  to  one  of  three  depths:  the  top  surface,  middle  of  the  tissue  equivalent  material  or 
the  back  of  the  head  phantom  (Figure  2-1).  The  temperature  was  read  for  approximately 
two  minutes  or  until  the  reading  on  the  Atkins  unit  stabilized.  After  the  initial  temperature 
was  determined  the  thermocouple  was  removed  and  the  head  phantom  scanned  for  a 
predetermined  time  using  the  modified  Phantom  Protocol.  At  the  conclusion  of  the  scan, 
the  phantom  was  removed  fi'om  the  MRI  unit  and  the  thermocouple  reinserted. 
Temperature  measurements  were  generally  obtained  withm  30  seconds  of  termination  of 
the  scan  and  were  subsequently  recorded  at  0. 1  T  intervals  at  one  minute  intervals.  To 
improve  results,  fiber  optic  cables  with  ten  meter  extensions  were  purchased. 
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Fiber  Optic  Measurements 

Three  fiber  optic  probes  were  purchased  to  replace  the  thermocouple  cables.  The 
Fluoroptic®  Thermometer  used  a  patented  fiberoptic  probe  manufactured  by  Luxtron. 
The  tips  of  the  fiber  optic  cables  use  the  fluorescence  of  phosphors  to  sense  temperatures. 
The  probes  are  immune  to  electromagnetic  interference  and  can  operate  with  the  -200  °C 
to  450  °C  range.  Since  the  fiber  optic  probes  were  immune  to  RF  effects,  the  probes 
remained  in  the  phantom  during  the  scan.  The  probes  were  linearly  positioned 
approximately  four  centimeters  fi"om  each  other.  The  exact  location  of  the  probes  was 
determined  for  each  experiment  through  imaging  the  head  phantom.  Using  the  landmark 
as  a  common  reference.  The  hole  at  the  right  eye  was  plotted  sagitally  and  axially.  The  tip 
of  each  probe  was  then  plotted.  To  account  for  any  rotation  of  the  head  phantom  between 
testing  dates,  the  outer  periphery  was  plotted. 

Using  the  drill  hole  location  and  the  outer  bounds  of  the  head  phantom,  the 
distance  from  the  surface  to  the  probe  head  was  calculated.  This  distance  was  defined  in 
millimeters  and  represented  the  "concentric  ring"  (Figure  2-2). 

Once  the  locations  of  the  probes  were  determined,  the  temperature  measurements 
were  recorded.  Some  experiments  were  performed  on  weekends  allowing  the  head 
phantom  twenty-four  hours  to  come  into  thermal  equilibrium  with  the  room  which  ranged 
form  17.2  to  17.8  °C.  Even  allowing  forty-eight  hours  for  the  head  phantom  to  come  into 
equilibrium,  the  head  phantom  would  still  have  diflfering  internal  temperatures.  On  the 
average,  the  temperature  varied  0.2  °C  between  points  within  the  head  phantom.  To 
compensate  for  variations  between  probes  and  starting  temperatures,  each  probe  was 
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evaluated  independently.  Equations  were  derived  using  absolute  temperature  variations 
only. 

Heating  Protocol 

A  head  phantom  imaging  protocol  had  been  previously  programmed  into  the 
ANMR  console.  Appendix  D  lists  the  program's  parameters.  The  output  power  was 
adjusted  by  increasing  or  decreasing  the  Transmission  Repetition  (TR)  value.  The  TR  was 
the  length  of  time  between  subsequent  starts  of  a  pulse  sequence.  The  phantom  protocol 
used  a  spin  echo  which  had  a  90°  pulse  followed  by  a  180°  pulse  (Figure  2-3).  The  TR 
was  the  length  of  time  between  subsequent  90°  pulses. 

By  decreasing  the  TR  value,  the  pulses  cycled  in  faster  succession  which  increased 
the  average  power  being  transmitted.  The  majority  of  experiments  were  conducted  at  an 
average  power  of  1 8  watts  as  indicated  on  the  Omega.  The  head  phantom  was  placed  into 
the  GE  Head  Coil  in  the  same  orientation  as  a  patient.  The  "landmark"  of  the  imager  was 
set  on  the  bridge  of  the  nose  of  the  head  phantom.  This  represented  the  iso-center  of  the 
head  phantom  and  represented  the  spatial  location  of  x=0,  y=0  and  z=0  of  the  imaging 
chamber.  The  patient  carriage  containing  the  GE  Head  Coil  and  head  phantom  was  then 
sent  into  the  imaging  chamber.  The  Phantom  Protocol  was  selected  on  the  ANMR 
console  and  the  TR  adjusted  to  obtain  a  desired  average  power  deposition.  A  8  mm 
camcorder  was  positioned  to  record  the  fiber  optic  probe  results,  the  average  power  as 
indicated  on  the  Omega  and  the  current  time.  The  camcorder  images  were  used  in  post 
imaging  analysis  to  enter  time,  average  power  and  temperature  into  a  database  program. 
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Thermal  Imaging 
MRI  affects  tissue  by  depositing  thermal  energy.  This  thermal  energy  can  alter 
some  parameters  measurable  by  MRI.  Temperature  resolution  using  MRI  was  first 
applied  to  imaging  hypothermia  patients.  For  hypothermia  patients,  MRI  had  the  dual 
benefit  of  providing  heat  directly  to  the  body  core  and  imaging  the  temperature  changes 
within  the  body.  Unwanted  heat  deposition  was  also  the  problem  for  general  imaging. 
The  MRI  system  can  be  self-regulating  by  monitoring  heat  deposition  during  a  scan. 
Imaging  relies  on  the  detection  and  quantification  of  signals  fi"om  the  tissue.  Thermal 
imaging  described  in  literature  required  a  reference  thermal  image  (Delannoy  1991,  Hall 
1990,  Schwarzbauer  1995).  The  temperature  was  then  defined  for  this  reference  image. 
This  method  required  extensive  calibration  within  each  tissue  type.  Establishing  a 
reference  temperature  had  an  inherent  accuracy  problem  which  propagated  into  a  larger 
error  during  subsequent  imaging.  To  reduce  this  error,  the  first  difiusion  image  only 
defined  the  baseline.  Regardless  of  the  starting  temperature,  any  differences  in  subsequent 
images  would  provide  a  difference  in  temperature  fi"om  the  reference  image.  Quantifying 
temperature  changes  in  this  method  should  reduce  errors  and  increase  accuracies.  There 
were  at  least  three  characteristics  of  protons  that  lent  themselves  to  measuring  changes  in 
temperature:  Tl  relaxation  time,  difiusion  and  proton  chemical  shift. 

Tl  Relaxation  Time  Imaging 

The  Tl  relaxation  time,  also  known  as  the  Spin  Lattice  or  Longitudinal  Relaxation 
Time,  was  defined  as  the  characteristic  time  constant  for  the  return  of  the  proton's 
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longitudinal  axis  to  alignment  with  B^,  the  main  magnetic  field.  A  Tl  map  of  the  target 
tissue  can  be  generated  by  an  RF  pulse  of  sufficient  fi"equency  and  strength,  acting  in  a 
direction  90°  fi-om  B^  (Figure  2-4).  Shutting  off  the  transmitter  and  using  the  transmit 
antenna  as  a  receiver,  the  signal  strength  was  a  fiinction  of  the  quantity  of  the  nuclear 
spins  realigning  with  B^  and  provided  the  Tl  map.  After  absorption  of  the  RF  energy,  the 
absorbed  energy  dissipated  into  the  surrounding  molecules  at  an  exponential  rate.  After 
one  Tl  period,  63%  of  the  protons  are  realigned  with  B^,,  also  known  as  their  rest  value. 
After  two  Tl  periods,  86%  have  realigned  and  after  three  Tl  periods,  95%  have  realigned 
(Figure  2-5).  As  temperature  increased,  the  relaxation  time  decreased  (Nelson  1987). 
Normally  the  Tl  for  protons  in  water  was  two  seconds. 

For  thermal  imaging,  a  reference  image  would  be  taken  at  a  known  temperature. 
After  tissue  heating,  another  Tl  image  was  taken  and  the  images  compared.  The 
difference  in  Tl  times  defined  the  temperature  change.  One  aspect  of  thermal  imaging 
was  no  change  in  temperature  of  the  cerebral  spinal  fluid  and  blood  from  short,  localized 
heating  tests.  The  fluids  present  during  a  slice  selective  RF  pulse  would  have  moved  to 
another  physical  location  by  the  time  a  second  image  was  taken.  The  replacing  fluids,  not 
having  absorbed  any  RF  energy,  would  already  be  at  their  rest  value.  The  temperature 
would  appear  to  remain  constant.  Early  in  the  development  of  Tl  mapping  a  sensitivity  of 
about  2  °C  was  obtained  (Parker  1983).  Since  the  FDA  standards  that  the  MRI  system 
must  abide  by  restricted  temperature  increases  to  only  1  °C,  this  method  would  prove 
unuseable.  The  sensitivity  was  directly  related  to  the  scan  time  which  lasted  five  minutes. 
During  this  scan  period,  the  entire  tissue  continued  to  absorb  RF  energy  which  in  turn 
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continued  to  increase  the  temperature.  The  solution  to  making  Tl  imaging  a  viable 
thermal  indicator  was  to  shorten  the  scan  time. 

Using  an  inversion-recovery  Snapshot  FLASH  MRI  technique  reduced  acquisition 
times  to  approximately  two  seconds  (Deichmann  1992,  NekoUa  1992).  An  inversion 
recovery  technique  began  with  a  1 80°  pulse  which  inverted  the  magnetization  and  flipped 
the  proton  axial  spin  opposite  to  B^.  After  a  time  period,  t,  a  90°  pulse  would  be  applied 
(Figure  2-6).  Transverse  magnetization  would  be  detected  (Figure  2-7).  The  inversion- 
recovery  Snapshot  FLASH  MRI  technique  took  images  within  milliseconds  of  the  RF 
pulse  by  using  smaller  flip  angles.  Since  smaller  flip  angles  require  less  RF  energy, 
temperature  increases  were  minimal.  One  sacrifice  for  this  rapid  measurement  was  a 
reduction  of  signal-to-noise.  This  can  be  compensated  for  by  averaging  several  images  in 
quick  succession.  The  precision  of  the  Tl  measurement  had  improved  to  better  than  three 
percent  with  this  method  (Schwarzbauer  1995).  A  three  percent  precision  at  a  normal 
body  temperature  of  37  °C  was  1 . 1  °C.  This  method  would  still  not  satisfy  FDA 
requirements  since  a  reading  of  0  °C  change  could  in  reality  exceed  the  1  °C  limit. 

In  addition,  temperature  changes  were  not  the  only  factors  affecting  the  Tl 
relaxation  times.  The  presence  of  fi"ee  radicals  and  water  concentration  also  affected  Tl 
(Knuttel  1986).  These  factors  became  a  concern  when  their  concentrations  changed  fi'om 
the  time  the  reference  set  of  images  were  taken  untU  the  time  the  final  Tl  images  were 
taken.  A  larger  quantity  of  protons  in  the  second  Tl  image  would  increase  the  apparent 
percentageofrealignment  of  nuclear  spins  with  Bq.  This  increased  quantity,  and 
subsequent  apparent  increased  realignment  would  promote  a  misconceived  perception  of 
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increased  tissue  temperature.  The  MRI  system  only  generates  non-ionizing  radiation,  so  it 
was  unlikely  the  generator  of  it  created  additional  free  radicals  during  the  imaging  process. 
Since  the  free  radicals  do  possess  a  polarity  which  was  subject  to  magnetic  fields,  only  if 
the  free  radicals  were  uniformly  distributed,  whereby  the  movement  of  radicals  out  of  a 
particular  location  matched  the  movement  of  radicals  into  the  particular  location,  could 
the  effects  of  free  radicals  be  ignored. 

Two  additional  factors  that  affected  the  Tl  relaxation  times  were  perfiision  and 
diffusion.  Perfiasion  was  the  pseudo-random  flow  at  low  velocities  of  blood  moving  along 
finely  divided  structures  of  the  capillaries.  Diffusion  was  the  Brownian  motion  of 
individual  molecules  moving  with  large  random  thermal  velocities.  Brownian  motion  was 
the  constant  erratic  movement  of  particles  being  struck  by  other  particles.  Temperature 
changes  affected  the  rate  of  Brownian  motion  by  affecting  the  viscosity  and  diffusion  of 
water  in  tissue. 

Diffljsion  Imaging 

In  biological  tissue,  a  fraction  of  the  water  movement  in  capillaries  was  affected  by 
diffusion  (LeBihan  1988).  Water  flowing  down  a  linear  path  has  an  overall  direction 
down  the  path,  but  upon  closer  attention,  the  water  molecules  collide  and  deflect  off  each 
other  creating  a  random  motion.  This  interaction  was  calculated  by  differentiating  the 
Stokes-Einstein  relation  between  viscosity  and  diffusion. 
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where  k  =  Boltzmann  constant, 

E,  =  the  activation  energy  for  translational  molecular  difiusion, 
D  =  the  translational  molecular  diffusion  and 
T  =  temperature. 

A  temperature  change  map  can  be  derived  from  two  diffusion  images  (LeBihan 
1989). 
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An  initial  reference  diffusion  image,  D^,  was  required  as  well  as  the  initial  temperature. 
Using  E,  =  0.21  eV,  and  applying  a  1  °C  increase  in  temperature,  the  diffusion  coefficient 
increased  by  2.4%.  This  value  was  dependent  on  an  accurate  E,.  For  water  the  activation 
energy  remained  relatively  constant  at  0.21  eV.  For  tissue,  this  value  may  be  different  and 
should  be  determined  for  each  tissue  type  present  in  the  image.  If  this  activation  energy 
varied  as  a  function  of  either  Tl  and/or  T2,  an  image  weighted  in  favor  of  one  or  the  other 
would  be  required.  The  intensity  of  the  resuhant  image  would  sufiBce  to  assign  E,  values 
to  each  tissue  type.  Then  a  pre-  and  post-exam  diffusion  image  sequence  would  map  the 
temperature  changes. 

Both  perfusion  and  diffusion  images  can  be  obtained  using  the  Intravoxel 
Incoherent  Motion  (IVIM)  method.  The  velocities  of  molecules  in  a  particular  tissue  type 
were  directly  related  to  the  kinetic  energy  of  the  molecules.  This  kinetic  energy,  in  the 
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form  of  thermal  energy  can  be  imaged  using  a  spin-echo  sequence  followed  by  a  second 
sequence  with  increased  gradient  pulses.  The  increased  gradient  pulses  increase  the 
effects  of  the  perfiision  and  diffusion.  Since  the  pulse  sequences  are  identical,  except  for 
the  gradient  pulses,  Tl  and  T2  effects  were  negated.  This  leaves  only  an  apparent 
diffusion  coefficient  (ADC)  which  was  a  combination  of  both  perfusion  and  diffusion.  A 
third  sequence  with  an  even  stronger  gradient  pulse  greatly  reduced  perfiision  while  only 
moderately  reducing  diffusion  effects.  A  comparison  between  the  second  and  third 
sequences  defined  a  diflfiasion  image.  A  comparison  of  the  diffusion  image  with  the  ADC 
image  fi^om  the  first  two  sequences  produced  a  perfiision  image.  Temperature  changes  do 
not  greatly  affect  perfiision.  Perfiision,  the  blood  flow  through  capillaries,  may  increase 
from  increased  cardiac  action  caused  by  the  deposition  of  heat  elsewhere  in  the  body.  The 
body's  response  to  thermal  changes  through  a  long  imaging  session  were  possible  but 
would  only  reflect  an  overall  temperature  increase  throughout  the  body.  Appendix  E, 
"Heat  Stress  Index",  described  this  effect  in  greater  detail.  Localized  heating,  which  was  a 
consideration  for  imaging  at  3T,  would  not  be  accurately  represented  by  perfiision 
imaging. 

The  rVIM  method  was  sensitive  to  motion  artifacts.  For  head  imaging,  cardiac 
gating  would  reduce  artifacts  fi-om  blood  flow,  and  to  a  lesser  degree,  cerebral  spinal  fluid 
(CSF)  flow.  Eye  movement  artifacts  would  be  reduced  by  preplanning  the  direction  of  the 
phase-encoding  gradient  pulses.  Increased  IVIM  images  using  faster  techniques  would 
also  reduce  motion  artifacts.  The  brain  was  a  non-moving  organ  and  thermal  imaging 
resolution  would  be  possible  as  long  as  patient  movement  was  reduced.  The  eye,  as  a 
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result  of  its  limited  vascularization,  was  the  tissue  most  susceptible  to  thermal  damage  and 
the  most  difficult  to  image  using  the  IVIM  method.  An  alternative  would  be  to  use  echo- 
planar  imaging. 

Echo-planar  imaging  used  a  gradient  echo  technique  consisting  of  a  rapid  sequence 
of  phase  steps  between  each  repeated  pulse  sequence.  The  echo-planar  imaging  provides 
rapid  images,  from  15  to  30  images  a  second.  This  produced  virtually  real-time  images  for 
functional  studies  of  the  brain.  During  brain  activity  there  is  a  rapid,  momentary  increase 
in  the  blood  flow  to  the  active  area  of  the  brain.  Therefore,  echo-planar  imaging  allows 
mapping  of  brain  activity.  Echo-planar  imaging  also  becomes  a  benefit  to  limit  motion 
artifacts  by  imaging  in  a  fraction  of  a  second.  Since  the  echo-planar  imaging  gradients 
follow  in  rapid  succession,  the  signal-to-noise  decreased.  However,  the  large  quantity  of 
data  points  from  the  multiple  sequences  improved  the  estimation  of  the  perfiision  and 
difiusion  (Turner  1990).  With  decreased  imaging  times,  diffusion  imaging  without  regard 
for  perfusion  artifacts  may  also  become  possible. 

Proton  Chemical  Shift  Imaging 

A  third  method  of  non-invasive  temperature  measurements  used  chemical  shift 
imaging.  Fat  protons  and  water  protons  have  differing  chemical  shifts  when  imaged  in  a 
MRI.  Generally,  chemical  shift  imaging  was  used  as  fat  suppression  imaging  .  The  pure 
water  proton  chemical  shift  was  proportional  to  a  change  in  temperature  by  approximately 
-0.01  ppm/°C  (Ishihara  1995).  This  temperature  dependence  resuhed  from  the  stretching, 
rupturing  and/or  bending  of  the  hydrogen  bonds.  Temperature  changes  were  difficult  to 
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detect  using  the  chemical  shift  because  of  the  small  shift.  An  accuracy  of  1  °C  was  usually 
expected  which  was  at  the  FDA  allowable  threshold  of  1  °C  change  in  tissue  temperature. 
This  meant  that  a  reading  of  1  "C  from  the  chemical  shift  could  actually  be  from  0  °C  to  2 
°C.  The  chemical  shift  was  detected  by  measuring  the  time  between  the  90°  pulse  and 
180°  pulse  and  comparing  this  to  the  difference  in  time  between  the  180°  pulse  and  the 
spin  echo  (Figure  2-3).  The  time  diflFerence  between  sets  of  images  was  the  resuh  of 
temperature. 

Chemical  shift  imaging  used  the  echo  planar  imaging  technique  to  help  reduce 
motion  artifacts  from  difiusion-weighted  techniques.  However,  the  echo  planar  imaging 
had  a  strong  sensitivity  to  magnetic  susceptibility.  Even  though  human  tissue  was  mostly 
not  affected  by  the  magnetic  field,  at  3T  even  small  differences  in  the  electrical  properties 
of  tissue  caused  a  dephasing  of  spins  and  frequency  shifts.  This  made  the  chemical  shift  of 
a  particular  tissue  weighted  by  the  concentration  of  KCl  and  NaCl  in  the  tissue.  To  use  an 
absolute  temperature  difference  between  the  same  tissue  over  a  heating  period,  the 
concentration  of  electrolytes  would  have  to  be  assumed  constant  over  a  temperature 
range.  Thechemicalshiftofprotons  was  also  affected  by  the  pH  of  the  tissue.  Since  the 
MRI  was  non-ionizing,  hydrogen  ions  were  not  created,  but  since  they  are  polarized  they 
may  be  moved  from  a  region  of  high  density  to  a  region  of  low  density.  Again,  for 
absolute  temperature  measurements  between  two  time  frames,  the  pH  should  not  change 
rapidly  and  the  reference  set  of  images,  if  taken  just  prior  to  the  exam,  would  still  provide 
a  good  baseline  for  comparison  of  post-exam  chemical  shift  images. 
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Phantom  Thermal  Imaging 

Of  the  three  thermal  imaging  techniques  described  above,  difiusion  imaging  using 
echo  planar  imaging,  appeared  the  best  alternative.  Neither  Tl  nor  chemical  shift  imaging 
v^as  solely  affected  by  temperature  and  both  were  affected  by  diffusion.  This  would 
require  techniques  to  eliminate  diflEusion  effects.  The  best  method  to  do  that  was  from 
two  sets  of  difiusion  weighted  images  to  remove  the  diflRasion  affects.  This,  in  effect, 
relied  on  running  diffusion  techniques  in  addition  to  other  techniques  to  image  Tl  or  the 
chemical  shift  as  a  fimction  of  temperature.  Since  diffusion  changes  from  one  image  to 
another  was  solely  affected  by  temperature,  using  just  diffusion  techniques  reduced 
imaging  time  thereby  decreasing  the  effect  from  RF  deposition  in  tissue. 

The  benefit  of  using  diffusion  imaging  was  three-fold.  First,  with  less  techniques, 
less  RF  energy  deposition  took  place  during  temperature  measurements.  Second,  with 
shorter  imaging  techniques,  results  would  be  easier  to  post  process.  Third,  while  not 
applicable  to  phantom  studies,  with  shorter  imaging  techniques  the  effects  of  patient 
motion  would  be  reduced. 
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Figure  2-1.  Thermocouple 
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Figure  2-2.  Concentric  Rings  Used  to 
Define  Location  in  Head  Phantom. 


'  I 


*  V 


180^ 


44 


90' 


go*^ 


!■< TR 3H 


Figure  2-3.  Spin  Echo  Sequence. 
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Figure  2-4.  90°  RF  Pulse. 
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Figure  2-5.  Tl  Recovery  Time. 
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Figure  2-6.  Inversion  Recovery  Pulse. 
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Table  2-1 .  Electrical  Properties  For  Polyethylene  Husk  of  Head  Phantom. 


Loss  Tangent 

Permittivity 

Intrinsic  Impedance 

Skin  Depth 

o 

[   F) 

e« 

(F] 

j](Q) 

6{m) 

0.0005 

2.26 

27.9  /0.2 

0.248 

Table  2-2.  Ingredients  to  Create  Tissue  Equivalent  Material. 

7.5%TX-51 

290g/gal.  Water 

0.45%  NaCl 

17.44g/gal.  Water 

0.07  mole  Gadopentetate  Dimiglumine 

0.544  ml/gal.  Water 

Table  2-3.  Thermal  Properties 

5  of  Tissue  Equivalent  Material  in  Head  Phantom. 

NaCl  Ratio 

Thermal  Conductivity 
(W/m°K) 

Specific  Heat 
(J/cm'TC) 

0.000765 

0.674  ±  0.025 

3.52 

0.00306 

0.646  ±  0.097 

3.53 

0.00461 

0.629 

3.54 

Table  2-4 

Electrical  Properties  of  Tissue  Equivalent  Material  in  Head  Phantom. 

Loss  Tangent 

Permittivity 

Intrinsic  Impedance 

Skin  Depth 

o 

i   F) 

T1(Q) 

6(/n) 

2.71 

61 

34.3  /32,6 

0.06 

CHAPTER  3 
RESULTS 


Three  sets  of  results  were  obtained  from  this  project.  The  first  set  was 
thermocouple  resuhs.  These  results  were  recorded  just  prior  to  heating  with  the  MRI 
system  and  then  immediately  after  the  scan.  The  thermocouple  wire  generated 
interference  and  increased  temperatures  when  left  in  the  head  phantom  during  the  scan. 
The  second  set  of  results  were  from  the  fiber  optic  cables  which  were  immune  to  RF  fields 
and  provided  information  during  the  heating  of  the  head  phantom.  The  third  set  of  data 
was  from  a  difiusion  imaging  technique.  Previous  thermal  data  from  invasive  temperature 
measurements  determined  the  difiusion  technique  could  provide  accurate  temperature 
resuhs. 

Thermocouple  Resuhs 
All  plots  from  experiments  with  thermocouple  probes  displayed  a  combination  of  two 
curves  (Figure  3-1).  The  predominant  curve  during  the  first  minutes  resuhed  from  the 
creation  of  a  "depletion  zone."  After  insertion  the  thermocouple  wire  became  a  heat  sink 
for  the  surrounding  tissue  equivalent  material  and  accelerated  the  heat  loss  for  the  first 
two  to  four  minutes.  Using  the  plots  from  the  first  four  minutes,  an  equation  was  derived 
of  this  initial  cooling.  Since  the  first  recorded  temperature  was  approximately  thirty 
seconds  after  the  scan  ended  and  a  few  seconds  after  the  thermocouple  was  inserted,  the 
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initial  temperature  was  calculated  from  the  assumption  that  the  cooling  followed  an 

exponentially  decaying  function  over  time  with  the  following  general  format. 


Fo  is  the  unknown  initial  temperature,  F,i  is  the  known  temperature  approximately  thirty 
seconds  after  the  scan,  t,  is  the  thirty  seconds  of  thermal  decay  time  and  X  is  the  unknown 
thermal  decay  constant. 

From  data  gathered  after  the  initial  four  minutes,  a  second  exponential  curve  was 
derived  which  defined  the  expected  cooling  curve  of  the  phantom.  This  cooling  curve  was 
then  adjusted  to  compensate  for  the  initial  heat  loss  provided  by  the  thermocouple  wire. 
This  method  depended  on  the  validity  of  the  thermocouple  cable  acting  as  a  heat  sink. 

Temperature  measurements  were  made  just  below  the  top  surface  of  the 
phantom's  husk  and  through  the  tissue  equivalent  material  to  the  back  of  the  head 
phantom.  There  was  more  fluctuation  in  temperature  measurements  at  the  top  surface  of 
the  phantom  which  was  exposed  to  the  air.  The  back  of  the  head  phantom  was  placed  on 
a  cushion  which  provided  insulation  from  heat  loss.  Temperature  measurements  from  this 
region  provided  stable  results  (Table  3-1).  The  heating  of  the  phantom  from  RF 
transmissions  were  compiled  and  a  heating  equation,  as  a  function  of  the  average  power 
read  off  the  Omega  process  indicator  and  the  scan  duration,  was  derived. 

Temp.  Inc.  (T)  =  ((4E-4  X  Powei^  -  (l/9)E-4  X  Power')  X  Scan  Time)(5/9) 
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This  equation  provided  close  approximations  to  actual  thermocouple  results  up  to 
thirty  watts  average  power.  Complications  were  found  when  measuring  temperatures  in 
the  middle  of  the  head  phantom.  After  the  scan  was  complete  and  the  thermocouple 
installed,  temperatures  continued  to  rise.  This  indicated  that  thermal  movement  would 
become  an  issue  of  concern  when  mapping  temperatures.  By  suggestion  this  also 
indicated  the  RF  energy  did  not  penetrate  the  entire  phantom  during  a  scan.  If  it  had,  then 
core  temperatures  would  have  matched  the  outer  surfaces.  This  demonstrated  the 
phenomenon  of  "skin  effect."  At  high  frequencies,  a  current  carrying  conductor  maintains 
a  flux  on  the  outer  edges  only.  This  theory  is  used  by  the  power  industry  which  allows 
them  to  use  high  tensile  strength,  but  low  current  carrying  potential,  material  to  string  long 
power  lines.  Only  the  outer  sheath  requires  current  conductor  properties.  This  same 
effect  was  demonstrated  during  thermal  testing  of  the  head  phantom.  With  this 
complication  and  relying  on  the  subjective  interpretation  of  the  heat  sink  properties  of  the 
thermocouples,  accurate  results  with  thermocouples  were  suspect. 

Fiber  Optic  Results 
Three  fiber  optic  probes  were  simultaneously  inserted  into  the  head  phantom.  The 
fiber  optic  probes  provided  temperature  measurements  during  the  scan  to  better 
characterize  the  heating  of  the  head  phantom.  Appendix  F  contains  the  raw  data  and 
graphs  derived  from  phantom  testing  using  the  fiber  optic  probes.  The  average  power 
measurements  were  obtained  from  the  Omega  at  one  minute  intervals.  The  shaded 
portions  of  the  graphs  represent  the  period  the  imaging  was  performed. 
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Heating  Equation 

An  exponential  function  was  assumed  for  the  heating  of  the  tissue. 


-)a(2)t  ^ 


T     [l   -  e    '''' 

max  \ 


Two  assumptions  were  used  in  determining  this  function.  First,  the  tissue  equivalent 
material,  similar  to  human  tissue,  would  not  have  an  unlimited  capability  to  absorb  thermal 
energy.  At  some  point  the  tissue  would  stop  increasmg  in  temperature  as  energy  was 
applied.  This  implied  some  upper  threshold,  T,^.  Since  both  the  tissue  equivalent 
material  and  human  brain  contain  mostly  water,  this  upper  threshold  would  be  sunilar  and 
would  also  be  far  above  the  pain  threshold  of  an  individual.  Therefore,  the  determination 
of  this  threshold  was  not  significant  to  this  research.  Only  the  assumption  that  this 
threshold  was  above  the  current  operating  temperatures  was  required.  Second, 
temperature  increases  as  a  function  of  power  would  not  be  Imear  as  they  approached  the 
threshold.  As  the  temperature  increases  there  would  be  a  constant  decrease  in  the  amount 
and  rate  that  additional  energy  could  increase  temperatures.  This  implies  an  exponential 
function.  Since  the  function  changes  exponentially,  it  was  assumed  that  there  would  be  a 
thermal  half-life,  tj/j   Experimental  data  would  be  used  to  define  these  parameters. 

The  thermal  data  results  were  measured  using  three  fiber  optic  probes  and  the 
location  of  each  probe  tip  was  determined  using  the  scanned  image.  Table  3-2,  3-3  and  3- 
4  list  the  locations  of  the  probe  tips  in  relation  to  the  outer  edge  of  the  head  phantom  and 
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to  the  eye  hole  drilled  into  the  head  phantom.  There  was  no  probe  #1  cable  used  during 
the  testing.  Appendix  F  lists  the  results  of  the  temperature  readings  for  the  first  thirty 
minutes  of  the  experiment.  For  the  majority  of  the  experimental  scans,  the  thirty  minute 
period  included  the  entire  scanning  time  of  twenty-two  minutes  and  several  minutes 
afterward. 

There  were  four  effects  apparent  in  the  data.  Three  of  the  effects  dealt  with  a 
delay  time  which  was  the  time,  in  minutes,  since  the  beginning  of  the  scan  until  there  was  a 
noticeable  increase  in  temperature. 

1.  Thedelay  time  was  a  function  ofdistancefi'om  the  outer  surface.  This  was 
intuitive  once  the  other  effects  were  defined  and  validates  some  of  the 
thermocouple  results.  The  electric  field  attenuated  as  it  passed  through  the  tissue 
equivalent  material.  This  attenuation  deposited  energy  in  the  outer  layer  thereby 
heating  it  first.  The  remaining  energy,  deposited  farther  into  the  tissue,  heated  the 
tissue  equivalent  material  at  a  lower  rate. 

2.  The  temperature  gradient  within  the  head  phantom  caused  a  delay  by  "passing" 
thermal  energy  to  the  lower  temperature  tissue  equivalent  material.  Temperatures 
varied  0.5  °C  within  one  centimeter  of  the  probes  at  the  start  of  Test  #3.  The  plots 
showed  a  slower  increase  in  temperature  where  the  tissue  was  warmer.  This 
demonstrated  an  increased  thermal  conductivity  of  the  warmer  tissue  equivalent 
material.    The  resuhs  fi-om  Test  #9  also  demonstrated  this  pattern.  There  the 
temperature  gradient  between  probes  was  0.8  °C  over  two  centimeters.  The  higher 
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temperature  zone  increased  in  temperature  at  a  slower  rate  than  the  inner  core 
even  though  it  was  closer  to  the  skin  and  should  have  heated  first. 

3.  The  lower  the  starting  temperature  of  the  head  phantom  the  sooner  a  temperature 
increase  was  registered.  If  the  head  phantom  had  sufficient  time  to  remove  thermal 
energy  then  the  delay  time  in  seeing  a  response  from  a  scan  of  18  watts  average 
power  was  one  minute  per  centimeter  depth.  The  higher  the  starting  temperature, 
the  delay  time  increased  by  twenty-five  percent  for  every  one  degree  centigrade.  A 
comparison  of  the  experimental  results  fi-om  Test  #2,  which  had  a  higher  starting 
temperature,  to  that  of  Test  #1,  where  the  temperatures  were  stable,  demonstrated 
this  point.  With  just  under  a  one  degree  centigrade  change  in  starting 
temperatures,  delay  times  in  registering  a  temperature  rise  increased  fi^om  ten 
minutes  to  twelve  minutes  at  a  depth  of  74  mm. 

4.  After  the  completion  of  the  scan,  temperatures  at  the  inner  core  continued  to  rise. 
Test  #1  showed  a  continued  increase  in  temperatures  after  completion  of  the  scan 
at  62  mm  and  73  mm  depths.  The  temperature  rise  at  the  41  mm  depth  continued 
to  rise  for  approximately  5  minutes  before  stabilizing.  This  reinforced  the  concept 
that  the  electric  field  was  not  uniformly  distributed  throughout  the  head  phantom 
and  that  the  majority  of  the  thermal  energy  was  deposited  along  the  outer  edges. 
Several  experiments  were  run  at  lower  average  power.  These  scans  generated  less 

thermal  energy  and,  as  expected,  the  rate  of  temperature  increase  was  less.  Since  the 
temperature  increase  appeared  as  an  exponential  fiinction,  a  "thermal  half-life"  of 
temperature  increase  was  derived.  From  test  results,  this  value  appeared  as  a  fiinction  of 
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the  average  power  (Table  3-5).  These  values  are  subjective  and  derived  only  from  the 

appearance  of  the  curve.  A  plot  of  the  natural  logs  of  the  half-life  against  the  average 
power  allowed  a  slight  modification  of  the  numbers  to  create  an  exponential  function 
relating  average  power  to  half-life  (Figure  3-2). 


/  _       „  6.375  -  0.134  *  (avg  power) 

'l/2  *^ 


Using  this  relationship  between  thermal  half-life  and  average  power  reading,  as  displayed 
on  the  Omega,  a  heating  equation  was  defined. 
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Cooling  Equations 

Cooling  equations  would  provide  additional  information  on  the  thermal 
characteristics  of  the  head  phantom  tissue  equivalent  material.  Ideally  the  tissue 
equivalent  material  would  begin  cooling  at  the  end  of  the  scan.  This,  however,  did  not 
take  into  account  the  thermal  movement  within  the  tissue.  Thermal  measurements  taken 
with  fiber  optic  cables  (Figure  3-3  through  3-8)  show  temperatures  taken  after  the  scan 
was  complete  and  no  fiirther  energy  was  deposited  into  the  head  phantom.  The  tissue 
equivalent  material  closer  to  the  outer  edge  of  the  head  phantom  began  reducing 
temperatures.  The  tissue  equivalent  material  closer  to  the  core  continued  to  increase  in 
temperature. 
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The  increased  temperatures  made  creating  a  cooling  curve  impossible.  Some 
useful  information  was  still  obtained  from  these  results.  Thermal  energy  traveled  in  the 
direction  of  cooler  tissue  equivalent  material.  This  direction  was  towards  the  center  of  the 
head  phantom.  This  implied  that  the  polyethylene  outer  husk  of  the  phantom  provided 
insulation  to  the  tissue  equivalent  material  and  minimal  heat  transfer  to  the  atmosphere. 
This  was  also  demonstrated  when  temperatures  had  not  stabilized  within  the  head 
phantom  after  twenty-four  hours  of  sitting  at  room  temperature. 

Diflusion  Imaging 
Diffusion  imaging  was  preferred  over  other  imaging  techniques  since  the  only 
parameter  that  affected  it  was  temperature.  Other  imaging  techniques,  such  as  T2- 
weighted  and  perfusion,  also  required  a  diffusion  image  technique  to  remove  the  diffusion 
effects  from  the  image.  The  diffusion  technique  was  run  on  the  head  phantom  prior  to 
heating.  After  the  diffusion  technique  was  run  and  images  captured,  an  imaging  technique 
was  performed  on  the  head  phantom  that  generated  approximately  18  watts,  as  indicated 
on  the  Omega,  for  approximately  twenty-two  minutes.  Immediately  after  the  head 
phantom  was  heated,  another  diffusion  technique  was  performed.  A  comparison  of  the 
two  spatially  different  diffusion  techniques  would  indicate  changes  as  a  result  of 
temperature  alone. 

Procedures 

To  capture  the  difflision  images  required  the  InstaScan  computer  system 
connected  to  the  ANMR  unit.  Figure  3-9  displays  the  connections  between  the  InstaScan 
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computer,  the  ANMR  Console,  the  uGT  Sun  Workstation  and  the  Windows-based 
computer.  Also  listed  are  the  program  names  used  to  transfer  and  manipulate  images. 
Prior  to  difiusion  imagmg,  "ghosts"  needed  to  be  tuned  from  the  system.  Ghosts  were 
unwanted  artifacts  of  the  original  image.  Appendix  G  lists  the  steps  performed  on  the 
InstaScan  computer  and  the  ANMR  console  for  ghost  tuning. 

Difiusion  Imaging.  Once  the  ghost  images  were  attenuated,  the  diffusion  technique  was 
run.  Appendix  H  lists  the  steps  performed  on  the  InstaScan  computer  and  the  ANMR 
console  for  diffusion  imaging.  Appendix  I  lists  the  parameters  for  diffusion  imaging  on  the 
ANMR  console.  The  diffusion  technique  provided  seven  images.  The  first  image  was  a 
T2-weighted  image.  The  next  three  images  were  diffusion  images,  at  the  slice  location, 
for  the  three  vectors  X,  Y  and  Z.  The  final  three  images  were  another  set  of  diffusion 
images  with  increased  gradients.  For  phantom  studies  the  second  set  of  images  were  not 
required.  The  diffusion  technique  actually  provides  both  perfusion  and  diffusion  which 
was  described  earlier  as  the  ADC.  Perfijsion  is  blood  flow  through  small  capillaries  and  a 
second  diffusion  image  with  increased  gradients  reduces  the  perfiision  effects  more  than 
the  diffusion  effects.  Since  the  phantom  did  not  have  capillaries,  there  were  no  effects 
from  perfijsion.  However,  for  this  diffusion  technique  to  be  effective  on  humans,  a  second 
diffusion  technique  will  be  required.  If  desired,  the  second  set  of  images  can  be  subtracted 
from  the  first  set  of  images  to  provide  images  that  predominately  contain  perfiision  effects 
and  a  small  portion  of  diffusion  effects. 
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Image  Conversion.  The  diffusion  images  were  physically  located  on  the  InstaScan 

computer.  In  order  to  manipulate  them,  they  were  first  moved  back  to  the  A>fMR  console 

by  selecting  the  images  fi-om  the  database  on  the  InstaScan  computer  and  selecting 

TRANSFER.  From  this  location  they  were  retrieved  fi^om  the  uf3T  Sun  Workstation. 

Appendix  J  lists  the  steps  to  transfer  the  images  fi"om  the  ANMR  console  to  the  Sun 

Workstation.  Once  on  the  Sun  Workstation  the  files  were  converted  to  another  format 

that  was  transferred  to  a  Microsoft®  Windows-based  computer  for  fiirther  processing. 

Image  Manipulation.  Appendix  K  Usts  the  program  written  by  the  author  for  this  research. 

The  program  was  written  for  a  Windows®  95  platform  and  was  designed  to  take  a  bitmap 

image  and  digitize  each  pixel  into  a  long  integer  indicating  its  color.  The  location  of  the 

pixel,  given  in  X  and  Y  coordinates,  and  the  number  associated  with  its  color,  was  saved 

in  a  Microsoft®  Access  database.  Once  the  images  were  saved  into  the  database,  the 

program  would  take  the  images  from  the  first  diffusion  run  and  subtract  the  images  from 

the  second  diffusion  run.  The  result  represented  the  change  in  the  two  images.  Since  a 

change  in  diffusion  images  was  only  affected  by  temperature,  the  values  from  the  program 

characterized  temperature  changes. 

Figure  3-10  is  a  display  of  the  entire  set  of  images.  The  first  row  shows  the  seven 

images  from  the  diffusion  imaging  performed  prior  to  heating  the  head  phantom.  The 
second  row  shows  the  seven  images  from  the  diffusion  imaging  performed  after  a  twenty- 
two  minute  spin  echo  sequence  generating  approximately  1 8  watts  as  indicated  on  the 
Omega.  The  third  row  shows  the  comparison  of  the  two  images.  Since  the  image  transfer 
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process  involved  several  steps,  some  noise  was  expected  from  signal  degradation  in  the 
final  images.  As  a  check  for  undesired  degradation  of  the  images,  the  second  row  of 
images  was  saved  twice  from  the  InstaScan  computer.  This  was  accomplished  by 
transferring  the  same  set  of  images  twice  back  over  to  the  ANMR  system.  Each  set  was 
saved  as  the  next  sequential  series  number.  At  that  point  the  second,  duphcate  set  of 
diffusion  images  were  processed  identically  as  the  original.  The  images  were  transferred 
to  bitmaps  and  digitized.  The  two  sets  of  identical  images  were  then  subtracted  from  each 
other  and  displayed  in  the  fourth  row. 

Image  Analysis.  The  images  in  the  first  row  of  Figure  3-10  were  provided  by  the  difiusion 
technique  taken  prior  to  heating  the  head  phantom  (Figure  3-11).  The  same  heating 
technique  that  had  been  used  for  fiber  optic  tests  was  applied  to  the  head  phantom  to 
increase  its  temperature.  This  technique  applied  eighteen  watts  average  for  twenty-two 
minutes.  The  diflEusion  technique  was  run  for  a  second  time  and  the  results  displayed  in 
row  two  of  Figure  3-10.  The  first  diflftision  image  (the  second  image  in  the  series)  is 
displayed  in  Figure  3-12.  The  two  images  were  digitized  and  a  number  was  assigned  to 
the  color  of  each  pixel.  The  color  values  derived  for  each  pixel  of  each  image  were 
subtracted  to  provide  a  method  to  quantify  the  change  in  diffusion.  This  was  validated  by 
test  results  beginning  with  the  Atkins  thermocouple  unit  and  later  with  the  fiber  optic 
system.  The  fiber  optic  probes  indicated  that  the  temperature  increased  along  the  outer 
rim  of  the  head  phantom  before  temperatures  would  increase  deeper  in  the  head  phantom. 
This  would  translate  into  an  initially  greater  change  in  diffiision  along  the  outer  portions  of 
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the  head  phantom  but  that  by  the  end  of  the  scan  temperatures  would  be  similar  and  the 
overall  difiiision  image  would  be  more  uniform. 

The  numerical  value  which  defines  each  pixel  in  the  image  is  a  long  integer 
representation  of  its  color.  The  numerical  format  consisted  of  six  hex  digits  representing 
the  colors  red,  green  and  blue  with  each  color  having  a  magnitude  of  0  to  FF  (hex).  For 
instance,  a  pure  red  color  would  have  the  hex  number  OOOOFF,  a  pure  green  color  would 
have  the  hex  number  OOFFOO  and  a  pure  blue  color  would  have  the  hex  number  FFOOOO. 
Since  the  images  were  displayed  on  a  gray  scale,  the  values  for  red,  green  and  blue  were 
identical  (i.e.  white  equals  000000,  light  gray  equals  111111,  dark  gray  equals  AAAAAA 
and  black  equals  FFFFFF).  This  simplified  numerical  interpretation  by  allowing 
subtraction  of  only  one  color  which  had  a  range  of  0  to  255.  Changes  in  these  numbers 
would  indicate  the  magnitude  of  change  in  diffusion,  and  thereby,  change  in  temperature. 

Figures  3-13  through  3-15  represent  T2 -weighted  images  that  accompanied  the 
diflRision  image  technique.  The  images  had  noise  attached  to  the  images  as  they  were 
processed.  This  noise  transferred  as  darker  pixels  when  converted  to  mathematical  values. 
The  noise  was  later  reduced  by  programming  a  fiher  to  remove  numerical  values  below  50 
(on  a  scale  of  0  to  255).  Relevant  portions  of  the  scan  had  image  values  above  50.  The 
filter  value  was  derived  experimentally  through  trial  and  error.  Also  present  on  the  image 
was  a  horizontal  line.  The  line  represented  the  location  of  the  displayed  histogram. 

The  top  histogram  (Image  002  window)  was  taken  fi-om  the  pre-heating  diffusion 
image  and  the  middle  histogram  (Image  004  v^dndow)  was  taken  fi-om  the  post-heating 
diffusion  image.  Thescaleofthe  window  was  fi-om  0  to  255.  These  two  images  appeared 
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similar  so  a  composite  image  was  created.  The  histogram  from  the  post-heating  diffusion 
image  was  subtracted  from  the  pre-heating  diffusion  image. 

The  third  histogram  is  the  Composite  window  which  was  scaled  vertically  from 
100  to  -100.  Image  composites  where  the  post-heating  images  were  less  intense  than  the 
pre-heating  images  would  have  the  histogram  generated  upwards.  The  downward 
generated  histogram  indicated  that  the  intensity  of  the  pixel  increased  in  intensity,  which 
translated  into  an  increased  numerical  value  when  digitized  by  the  computer,  as 
temperatures  increased.  Three  slices,  as  indicated  by  the  horizontal  line  on  the  right  side 
image  of  Figures  3-13  through  3-24,  were  taken  from  each  axis  of  the  diffusion  image. 
Slices  were  taken  from  the  top,  which  represented  the  eye  and  nose  area  of  the  head 
phantom,  the  middle  and  the  lower  or  back  portion  of  the  head  phantom. 

The  Composite  images  from  all  subtracted  diffusion  results  displayed  histograms  in 
the  negative  direction.  This  indicated  the  intensity  of  the  images  increased  with  increased 
temperature.  Had  the  Composite  histogram  been  random,  where  equal  numbers  of  points 
would  have  been  positive  and  negative,  then  this  diffusion  imaging  technique  to  measure 
temperature  changes  could  be  considered  non-effective.  The  histograms  demonstrated 
diffusion  imaging  does  provide  a  means  of  measure  temperature  changes. 

The  diffusion  images  represent  diffusion  changes  in  three  orthogonal  vectors.  To 
quantify  one  numerical  results  for  each  pixel,  a  vector  sum  of  the  three  diffusion  images 
was  obtained.  For  this  research,  the  average  intensity  of  all  the  pixels  was  calculated  for 
each  vector.  The  average  was  obtained  from  adding  the  intensity  diflferences  of  the  three 
slices  taken  from  each  axis.  The  sum  of  the  Composite  images  from  Figure  3-16  through 
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3-18  provided  the  average  intensity  change  for  the  X-axis.  The  sum  of  the  Composite 
images  from  Figure  3-19  through  3-21  provided  the  average  intensity  change  for  the  Y- 
axis.  The  sum  of  the  Composite  images  from  Figures  3-22  through  3-24  provided  the 
average  intensity  change  for  the  Z-axis.  These  averages  were  summed  using  the  following 
equation. 

Pixel  Intensity    =    yjx'^  +  Y"^  +  Z'^ 

The  average  change  in  intensity  of  the  composite  images  displayed  in  Figures  3-16  through 
3-24  was  40.7  on  a  scale  of  0  to  255.  The  imaging  technique  used  for  heating  the  head 
phantom  generated  sufficient  thermal  energy  to  increase  temperatures  by  approximately 
0.5  °C.  These  results  indicated  temperature  changes  less  than  1  °C  can  be  measured. 


62 


Phantom  Cooling  Curve 


Time 


Figure  3-1 .  Cooling  Curve  Using  Thermocouple  Probes. 
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Figure  3-2.  Natural  Log  of  Thermal  Half-Life. 
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Figure  3-3.  Test  #4  Post  Scan  Temperature  Measurements. 
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Figure  3-5.  Test  #4  Post  Scan  Temperature  Measurements. 
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Figure  3-6.  Test  #6  Post  Scan  Temperature  Measurements. 
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Figure  3-7.  Test  #6  Post  Scan  Temperature  Measurements. 
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Figure  3-8.  Test  #6  Post  Scan  Temperature  Measurements. 
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Figure  3-9.  Difiusion  Imaging  Equipment  and  Software 
Flowchart. 
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IrAial  Dffution  Irrvages 


Figure  3-10.  Diffusion  Technique  Images  and  Changes  in  Images. 
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Figure  3-11.  Difiusion  Image  Prior  to  Heating  Head 
Phantom. 
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Figure  3-12.  Difiusion  Image  After  Heating  Head 
Phantom. 
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Figure  3-13.  Histogram  From  Slice  of  T2-Weighted  Image. 
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Figure  3-14.  Histogram  From  Slice  of  T2 -Weighted  Image. 
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Figure  3-15.  Histogram  From  Slice  of  T2-Weighted  Image. 


hnage  002-2 


l«n9sSe4-2 


tuiifmiiMhMi 


littAiiiiiiitiiiiiTHiiii 


fyr9*S  iijftiTf^if^i  Vfl*' 


Figure  3-16.  Histogram  From  Slice  of  Diffusion  Image  X-Axis. 
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Figure  3-17.  Histogram  From  Slice  of  Diffusion  Image  X-Axis. 
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Figure  3-18.  Mstogram  From  Slice  of  Diffusion  Image  X-Axis. 
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Figure  3-19.  Histogram  From  Slice  of  Diffusion  Image  Y-Axis. 
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Figure  3-20.  Histogram  From  Slice  of  Diffusion  Image  Y-Axis. 
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Figure  3-21 .  IBstogram  From  Slice  of  Diffusion  Image  Y-Axis. 
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Figure  3-22.  Histogram  From  Slice  of  Diffusion  Image  Z-Axis 
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Figure  3-23.  ffistogram  From  Slice  of  Diffusion  Image  Z-Axis. 
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Figure  3-24.  Histogram  From  Slice  of  Diffusion  Image  Z-Axis. 
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Table  3  - 1 .  Temperature  Measurements 

Using  Thermocouples. 

Date 

Power  (W) 

Location 

Start  Temp 
(F) 

Stop  Temp 
(F) 

Scan  Time 
(min) 

11/13/96 

30 

back 

70.6 

73.5 

48.22 

11/13/96 

30 

back 

73.1 

76.2 

48.22 

1 1/20/96 

31 

front 

69.5 

70.1 

10.02 

11/20/96 

31 

front 

68.7 

70.7 

10.02 

11/20/96 

31 

back 

69.8 

71.5 

10.02 

11/20/96 

31 

back 

71.1 

71.8 

10.02 

11/25/96 

30 

front 

70.5 

73.8 

48.42 

12/16/96 

17 

back 

68.8 

70.4 

32.07 

12/16/96 

17 

back 

70.2 

72.3 

32.07 

12/28/96 

14 

back 

67.7 

70.0 

42.73 

12/28/96 

14 

back 

69.4 

72.3 

42.73 

12/28/96 

14 

back 

70.9 

73.0 

42.73 

12/29/96 

14 

front 

68.1 

71.5 

42.73 

12/29/96 

14 

front 

70.9 

73.4 

42.73 

12/29/96 

14 

front 

72.5 

74.7 

42.73 

1/4/97 

9 

back 

67.2 

68.7 

42.73 

1/4/97 

9 

back 

68.3 

69.4 

42.73 

1/5/97 

37 

back 

69.7 

74.5 

42.45 

1/5/97 

37 

back 

74.0 

78.6 

42.45 

1/5/97 

37 

back 

77.9 

81.9 

42.45 
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Table  3-2.  Head  Phantom  Probe  #2  Location. 


Location  on  Drawing 

Concentric 
Ring  No. 

Test 

Dx 

Dy 

x 

y 

#1 

12.2 

26.3 

70 

52 

41 

#2 

14.1 

25.4 

72 

51 

40 

#3 

0 

22.5 

58 

49 

38 

#4 

0.9 

15 

59 

41 

30 

#5 

1.9 

13.1 

60 

39 

28 

#6 

3.7 

10.3 

62 

36 

25 

#7 

-8.1 

29.2 

50 

55 

30 

#8 

6.5 

10.3 

65 

36 

25 

#9 

-2.8 

19.6 

55 

46 

35 

Table  3-3.  Head  Phantom  Probe  #3  Location. 

Locat  on  on  Drawing 

Concentric 
Ring  No. 

Test 

Dx 

Dy 

X 

y 

#1 

24.4 

84.4 

82 

110 

62 

#2 

35.7 

65.7 

94 

92 

74 

#3 

9.4 

57.2 

67 

83 

47 

#4 

36.5 

J0.99 

94 

57 

46 

#5 

13.1 

62.8 

71 

89 

51 

#6 

4.7 

88.1 

63 

11A 

A1 

#7 

-14.1 

59.1 

44 

85 

24 

#8 

28.1 

50.6 

86 

77 

66 

#9 

1.9 

62.8 

60 

89 

40 

Table  3-4.  Head  Phantom  Probe  #4  Location. 

Location  on  Drawing 

Concentric 
Ring  No. 

Test 

Dx 

Dy 

X 

y 

#1 

34.7 

121.9 

93 

148 

73 

#2 

52.8 

95.3 

111 

121 

91 

#3 

18.8 

86.3 

77 

112 

57 

#4 

26.2 

67.5 

84 

94 

64 

#5 

15.0 

101.2 

73 

197 

K1 

#6 

3.7 

95.6 

62 

122 

42 

#7 

-30.0 

94.7 

28 

121 

8 

#8 

42.2 

80.6 

100 

107 

80 

#9 

16.0 

99.3 

74 

125 

54 

77 


Table  3-5.  Thermal  Half-Lif 

Average 
Power 

Half-Life 

5 

300 

10 

150 

20 

40 

«  Using  Fiber  Optic  Cables. 


CHAPTER  4 
DISCUSSION 


In  the  first  series  of  thermal  test,  thermocouples  were  considered  a  cost  eflFective 
method  of  testing.  Analysis  of  data  revealed  three  major  problems.  First,  the 
thermocouple  wires,  being  metallic,  interfered  with  the  imaging  of  the  head  phantom  by 
creating  artifacts.  For  use  with  difiusion  imaging  the  exact  locations  of  the  thermocouples 
are  required.  With  artifacts  on  the  images  it  was  difficuh  to  plot  the  exact  location  of  the 
thermocouple  tip.  Additionally,  the  artifacts  would  also  affect  the  difiusion  image.  This 
would  not  provide  an  accurate  indication  of  the  change  in  diffusion  as  a  function  of 
temperature,  especially  near  the  thermocouple  probes  which  were  the  temperature 
indicators. 

Second,  the  wires  conducted  the  RF  fields  and  increased  the  temperatures  in  the 
tissue  equivalent  materials  surrounding  the  wires.  This  created  a  false  indication  that  more 
energy  was  being  deposited  into  the  tissue  resulting  in  greater  temperature  increases.  The 
equation  of  average  power  as  indicated  on  the  Omega  and  length  of  the  scan  indicated  a 
greater  increase  in  temperatures  over  a  shorter  time  interval.  A  comparison  of  this 
equation  with  the  one  derived  fi-om  fiber  optic  system  results,  revealed  an  increase  in 
temperature  increases  in  the  tissue  equivalent  material  as  a  function  of  average  power  and 
scan  time.  From  the  standpoint  of  safety,  the  thermocouple  resuUs  had  an  inherent  safety 
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margin  by  over  estimating  temperature  increases.  This  margin,  however,  limited  the 
capabilities  of  3TMRI. 

Third,  after  the  thermocouples  were  removed  from  the  head  phantom  during  a 
scan  and  later  reinstalled,  the  wire  acted  as  a  heat  sink.  The  wires,  being  at  room 
temperature,  were  slightly  cooler  than  the  tissue  equivalent  material  of  the  head  phantom. 
Once  installed  into  the  head  phantom,  the  cooler  wire  acted  as  an  excellent  thermal 
conductor  and  rapidly  lower  temperatures.  Readmgs  of  1  °C  drop  in  temperature  within 
two  to  four  minutes  were  common.  One  method  to  reduce  this  would  be  to  shield  the 
thermocouple  end  in  plastic.  This  would  insulate  the  wire  from  the  tissue  equivalent 
material  but  would  also  slow  down  the  response  time  of  the  thermocouple.  Since  the 
cooling  of  the  head  phantom  is  slow  this  may  be  possible,  however,  the  thermocouple  wire 
would  still  conduct  heat  faster  until  it  reached  equilibrium  with  the  tissue  equivalent 
material. 

Fiber  optic  cables  were  immune  to  magnetic  and  RF  fields.  Three  fiber  optic 
probes  were  simultaneously  installed  in  the  head  phantom.  Measurements  of  temperatures 
were  possible  and  revealed  characteristics  that  would  never  have  been  seen  with 
thermocouple  wires.  Since  the  RF  output  power  was  constant  during  a  long  scan,  it  was 
assumed  that  temperatures  would  increase  Unearly.  The  majority  of  temperature  increases 
were  exponential.  This  exponential  curve  started  earlier  in  the  scan  the  closer  the  probe 
was  to  the  outer  edge  of  the  head  phantom.  This  demonstrated  that  RF  deposition,  and 
subsequently  temperature  increases,  were  greater  on  the  outer  edge.  Thermal  energy  did 
move  inward  and  at  the  end  of  a  scan  the  entire  head  phantom  had  increased  in 
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temperature.  Additionally,  test  results  from  Test  #9  showed  a  slower  rate  of  temperature 
increase  along  the  outer  edge  in  comparison  to  an  inner  core  probe.  This  validates  the 
assumption  that  the  capability  of  the  tissue  equivalent  material  to  increase  temperatures 
from  a  given  power  deposition  decreases  as  the  tissue  equivalent  material  increases  in 
temperature.  The  amount  of  energy  deposited  in  the  tissue  equivalent  material  nearer  the 
core  was  less  than  the  energy  deposited  near  the  outer  edge,  but  since  the  inner  tissue 
equivalent  material  was  lower  in  temperature  the  deposited  energy  increased  the 
temperature  more.  These  mechanisms  were  not  detectable  with  the  thermocouple  system. 
Fiber  optic  cables,  remaining  installed  during  a  scan,  will  be  instrumental  m  testing 
diffusion  images  at  different  stages  of  a  scan. 

Deriving  cooHng  curves  would  better  quantify  thermal  energy  deposition.  The 
temperature  of  any  volume  of  tissue  can  be  mathematically  defined  as  simultaneously 
receiving  energy  and  dissipating  energy  during  a  scan. 


^Volume  j\     IN        ^OUTJ 


After  a  scan,  with  no  source  of  energy,  the  temperature  of  a  volume  of  tissue  is  a  function 
of  energy  loss  only. 


Volume  j^OUT 
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Defining  the  energy  deposition  properties  of  the  3T  system  would  be  a  matter  of 

subtracting  the  cooling  equation  fi-om  the  heating  equation.  However,  during  these  tests, 

there  was  constant  thermal  movement  up  to  an  hour  after  the  scan. 

These  results  would  prove  useful  in  better  characterizing  the  tissue  equivalent 
material.  Unless  the  same  process  is  performed  on  human  tissue,  there  would  be  no 
correlation  and  the  information  would  not  be  useful.  Additionally,  the  human  head  has  a 
thermal  regulatory  mechanism  which  is  not  available  to  the  head  phantom.  A  better 
characterization  of  the  tissue  equivalent  material  would  not  improve  test  resuhs. 

Diffusion  imaging,  being  only  affected  by  temperature,  can  provide  non-invasive 
measurement  of  temperature  increases.  The  FDA's  guideline  limiting  temperature 
increases  to  any  1  gram  of  tissue  to  1  °C  can  only  be  met  through  imaging.  Any  other 
method  requires  mvasive  techniques  that  would  not  be  comfortable  to  patients.  Difiusion 
is  the  only  technique  developed  so  far  that  is  only  affected  by  temperature.  All  other 
imaging  techniques  are  affected  by  some  other  parameter  besides  temperature.  Until 
another  technique  is  developed  or  the  other  parameters  can  be  mitigated,  further  testing 
using  diffusion  imaging  is  warranted. 


CHAPTERS 
CONCLUSION 


For  the  results  obtained  during  this  research  to  be  valid  for  the  human  head,  the 
head  phantom  must  be  a  good  representative  of  the  human  head.  Physically  the  head 
phantom  is  an  identical  shape  of  a  human  head.  Its  electrical  loading  characteristics  were 
defined  through  an  average  loading  of  volunteers  from  3T  MRI.  Using  an  established 
recipe  and  retesting  in  the  3T  imager,  a  tissue  equivalent  material  to  simulate  the  human 
head  was  determined  (Olsen  1992,  Chou  1984).  The  primary  component  in  the  recipe  was 
water  which  is  also  the  primary  component  in  brain  tissue.  This  suggests  the  thermal 
characteristics  of  the  head  phantom  and  human  head  would  also  be  similar.  The  outer 
husk  of  the  head  phantom  was  an  insulator  similar  to  the  skull  and  skm  of  the  human  head. 
Both  contain  low  water  content  and  should  have  similar  thermal  characteristics. 
Fortunately,  with  the  use  of  difiusion  imaging,  exact  duplication  of  test  results  between  the 
head  phantom  and  human  head  may  not  be  required. 

Difiusion  imaging  had  proved  effective  in  determining  the  change  in  temperatures 
in-vivo  (LeBihan  1989,  Stehling  1991,  DePoorter  1994,  McCain  1995,  Mischler  1995). 
Results  mentioned  in  literature  indicated  an  accuracy  of  1.0  °C.  Temperatures  during  the 
difiusion  imaging  for  this  research  used  a  heating  technique  that  raised  temperatures  by 
0.5  °C  which  was  easily  noticeable  in  test  resuhs. 
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This  research  has  demonstrated  that  a  difilision  imaging  technique,  already 
available  on  MRIs,  can  be  used  to  measure  and  monitor  temperature  increases  in  the 
human  head  with  resolution  better  than  1  °C.  An  on-line  program,  similar  to  the  one 
written  for  this  research  would  provide  real  time  analysis  of  thermal  changes  by  utilizing  a 
quick,  initial  diffusion  image  and  then  subsequent  images  throughout  the  session  to 
compare  the  change  in  diffusion.  Further  testing  is  required  to  better  characterize  the 
temperature  increase  relationship  to  RF  fields. 

Future  Work 

Accuracy  in  determining  temperature  changes  through  diffusion  imaging  would 
entail  using  the  fiber  optic  probes  in  conjunction  with  the  diffusion  imaging  technique. 
The  result  would  be  a  description  of  the  change  in  pbcel  intensity  with  change  in 
temperature.  Temperature  gradients  existed  in  the  head  phantom  and  would  be  a  cause 
for  increased  errors  in  temperature  comparisons  between  the  diffusion  image  results  and 
the  fiber  optic  resuhs.  To  minimize  this  effect,  the  diffusion  image  should  be  run  on  the 
slice  containing  the  fiber  optic  probe  tips. 

A  variety  of  heating  techniques  resulting  in  a  variety  of  different  temperature 
increases  fi-om  0. 1  °C  to  2  °C  should  be  performed.  The  representative  diffusion  images 
would  then  be  transferred,  digitized  and  analyzed  using  software  similar  to  that  used 
during  this  research.  Since  noise  was  added  each  time  an  image  was  processed,  an 
analysis  program  should  be  written  for  the  InstaScan  unit  or  the  Sun  workstation.  This 
would  eliminate  the  need  to  transpose  the  images  fi-om  a  Unix  bitmap  to  a  portable  gray 
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map  to  a  Windows®  bitmap  and  finally  digitized.  A  program  written  for  the  InstaScan 
would  also  reduce  the  time  required  to  analyze  images. 

The  next  step  in  applying  this  technique  to  humans  would  be  to  run  the  diffusion 
imaging  technique  on  a  cadaver  head.  The  cadaver  head,  like  the  head  phantom,  should 
start  at  room  temperature.  Artificially  warming  the  cadaver  to  normal  body  temperature 
would  present  a  problem  of  maintaining  that  temperature.  Once  removed  fi^om  the 
warming  mechanism,  the  head  would  begin  to  cool.  The  change  in  temperature  would 
skew  the  temperature  changes  as  a  fimction  of  MRI  and  also  skew  difiusion  results. 
Temperature  increases  by  the  RF  fields  would  be  negated  by  the  cooling  of  the  head.  The 
amount  of  temperature  increase  as  a  fimction  of  the  length  and  strength  of  the  RF  fields 
would  have  to  be  assumed  as  similar  for  tissue  at  normal  body  temperatures  (37  °C)  and  at 
room  temperature  (18  °C). 

Fiber  optic  probes  would  be  inserted  into  the  various  bram  and  optic  tissues.  A 
difiusion  image  would  be  acquired  in  the  location  of  the  fiber  optic  probe  tips.  A  heating 
technique  would  then  be  applied  to  the  cadaver  head.  Post-heating  difiusion  images 
would  then  be  compared  m  the  same  manner  as  the  head  phantom  difiusion  images.  Only 
a  handfiil  of  images,  providing  different  temperature  changes,  would  generate  a 
relationship  between  tissue  type  and  pixel  intensity  as  a  fimction  of  temperature.  This 
pixel  intensity  verses  temperature  change  would  provide  a  conversion  factor  for  each 
tissue  type.  Most  brain  tissue  (white  matter  and  gray  matter)  have  a  high  water  content. 
Temperature  changes  as  a  fimction  of  the  RF  fields  would  most  hkely  be  similar.  As  long    . 
as  the  thermal  characteristics  of  the  brain  tissues  remain  similar  between  room  temperature 
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and  normal  body  temperature,  the  conversion  factor  would  be  valid  for  living  tissue  at 
normal  body  temperature. 


APPENDIX  A 
INTERNATIONAL  ELECTROTECHNICAL  COMMISSION  (lEC) 


Operating  Mode  Ranges  for  dB/dt 

Normal  Operating  Mode 

dB/dt  <  20  T/s  for  1 20  jis  <  T 

dB/dt  <  2400/T  T/s  for  2.5  us  <  T  ^  120  ^s 

dB/dt  <  960  T/s  for  T  <  2. 5  us 

Where  the  pulse  width,  T,  is  defined  as  the  duration  of  the  change  of  magnetic  flux 
density  expressed  in  microseconds  (|is).  For  sinusoidally  or  other  continuously 
varying  periodic  magnetic  fields  the  duration  of  change  shall  be  considered  to  be 
the  half  period  the  magnetic  flux  waveform. 

For  values  of  T  not  exceeding  2.5  us,  the  limit  on  peak  dB/dt  is  based  on  heating, 
assuming  a  sinusoidal  waveform.  Therefore,  the  limit  on  peak  dB/dt  of  960  T/s 
may  be  exceeded  for  non-sinusoidal  waveforms,  provided  that  the  r.m.s.  value  of 
dB/dt  is  less  than  960  T/s  and  the  peak  dB/dt  is  less  than  2400/T  T/s. 

First  Level  Controlled  Operating  Mode 

dB/dt  <  20  T/s  for  3000  ^is  <  T 

dB/dt  <  60000/T         T/s  for  45  us  <  T  ^  3000  jis 

dB/dt  <  1 330  T/s  for  T  <  45  jis 

For  values  of  T  not  exceeding  45  us,  the  limit  on  peak  dB/dt  is  based  on  heating, 
assuming  a  sinusoidal  waveform.  Therefore,  the  limit  on  peak  dB/dt  of  1330  T/s 
may  be  exceeded  for  non-sinusoidal  waveforms,  provided  that  the  r.m.s.  values  of 
dB/dt  is  less  than  1330  T/s  and  the  peak  dB/dt  is  less  than  60,000/T  T/s. 

Second  Level  Controlled  Operating  Mode 

Comprises  values  of  dB/dt  which  exceed  the  upper  limit  for  range  of  the  First 
Level  Controlled  Operating  Mode. 


86 


•-v3 


87 
Operating  Modes 

First  Level  Controlled  Operating  Mode 

Magnetic  Resonance  equipment  that  is  capable  of  values  of  dB/dt  and/or  SAR 
above  the  upper  limits  of  the  Normal  Operating  Mode,  shall  comply  with  the 
following: 

a)  Before  the  start  of  each  scan,  an  indication  of  the  operating  mode  defined 
by  the  maximum  dB/dt  value  for  the  scan  and  a  prediction  of  the  value  of  SAR  that 
will  actually  be  used  during  the  scan  shall  be  displayed  at  the  Operators  console. 

b)  If  the  value  of  dB/dt  or  SAR  is  such  as  to  enter  the  First  Level  Controlled 
Operating  Mode,  the  attention  of  the  Operator  shall  be  drawn  to  this  condition  by 
a  clear  indication  on  the  Operators  console  and  provision  shall  be  made  so  that  a 
deliberate  action  of  the  Operator  shall  be  necessary  to  start  the  scan.  A  record  of 
these  values  shall  be  inseparable  fi^om  the  image  data. 

c)  A  means  of  control  shall  be  provided  to  ensure  that  values  of  dB/dt  and 
SAR  do  not  exceed  the  upper  limits  for  the  First  Level  Controlled  Operating 
Mode.  This  control  shall  be  independent  of  Operator  input  as  the  patient  size, 
weight  or  position. 

Second  Level  Controlled  Operating  Mode 

Magnetic  Resonance  equipment  which  allows  values  of  dB/dt  or  SAR  within  the 
Second  Level  Controlled  Operating  Mode  shall  comply  with  the  requirements 
above  and  in  addition: 

a)  Shall  include  specific  security  measures  that  allow  entry  in  the  Second 
Level  Controlled  Operating  Mode.  These  measure  shall  include  an  indication  to 
the  Operator  that  the  operating  conditions  are  potentially  hazardous  and  that  these 
conditions  should  be  not  applied  for  normal  clinical  use. 

b)  The  specific  security  measures  shall  be  so  constructed  that  it  can  be  made 
accessible  only  under  the  authorization  of  the  medically  responsible  person. 

c)  The  specific  security  measures  shall  be  effective  for  each  scan  for  which  the 
predicted  values  of  dB/dt  or  SAR  are  above  the  limits  for  the  First  Level 
Controlled  Operating  Mode. 

d)  The  specific  security  measures  shall  involve  a  key-lock,  a  combination 
local,  a  software  password  or  other  protective  devices. 
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Checking  Compliance 

The  accuracy  of  the  predicted  SAR  with  the  specification  shall  be  checked  by 
measurement  at  the  maximum  SAR  of  which  the  Magnetic  Resonance  equipment  is 
capable  or  at  a  SAR  of  3W/kg,  whichever  is  lower. 

The  upper  Umit  of  the  First  Level  Controlled  Operating  Mode  shall  be  measured  to  check 
that  it  complies  with  the  maximum  SAR  of  which  the  Magnetic  Resonance  equipment  is 
capable  or  at  a  SAR  of  4  W/kg,  whichever  is  lower. 

The  requirement  to  prevent  exceeding  the  upper  limit  of  the  First  Level  Controlled 
Operating  Mode  for  dB/dt  shall  be  demonstrated  by  measurement  using  the  maximum 
gradient  rate. 


dB/dt  [T/s) 


Limits  for  dB/dt 


Range  of  Second  Level  Controlled  Operating  Mode 


Range  orFirst  Level  Controlled  Operating  Mode 


Figure  103  -  Limits  for  dB/dt 

Compliance  is  checked  by  measurement  of  the  maximum  trapezoidal  dB/dt  or  the 
maximum  sinusoidal  dB/dt  component,  that  the  system  is  capable  of  generating  under  any 
conditions  allowed  for  clinical  use  in  the  First  Level  Controlled  Operating  Mode. 

Method  of  measurement  of  dB/dt  to  demonstrate  compliance  with  the  given  requirements: 

Method  of  measurement  of  maximum  trapezoidal  dB/dt: 


a)         Auxiliary  definitions: 
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search  coil  a  small  diameter  coil  used  to  measure  gradient  field  rates  of 

change, 
magnet  isocentre  normally  defined  by  magnet  geometry,  in  this  document  it  is 

defined  as  the  zero  crossing  for  all  three  gradients. 
Patient  space  radius     the  distance  fi^om  magnet  isocentre  to  the  patient  support  (without 

pads)  along  the  lines  y  =  ±  x  (see  figure  103). 
dB/dt  pulse  width,  T:  the  Full  Width  at  Half  Maximum  of  a  single  dB/dt  pulse  or  the  half 
period  of  a  sinusoidal  dB/dt  waveform  (see  figure  104) 


Search  coil 
position  2 


Search  coil 
position  1 


Magnet 
isocentre 


Figure  103  -  Search  coil  placement  for  measurement  of  dB/dt 
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a:  magnetic  field  change  waveforms 
b:  dB/dt  waveforms 
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Figure  104  -  Magnetic  field  change  waveforms  and  dB/dt  waveforms. 

b)        Test  hardware: 


1) 


Search  coil 


The  search  coil  shall  be  circular  and  shall  be  small  with  respect  to  the  gradient  coils 
under  test  to  ensure  accuracy.  The  search  coil  consists  of  n  turns  of  wire  with  a 
radius  of  r  meters.  The  axial  length  of  the  coil  shall  be  less  than  20%  of  its 
diameter. 

The  search  coil  shall  be  no  more  than  50  mm  in  diameter.  The  response  of  the 
search  coil  shall  be  determined  by  calculation.  The  instantaneous  magnitude  of  the 
component  of  dB/dt  coaxial  with  the  search  coil  shall  be  determined  from  the  peak 
voltage,  V,  induced  in  the  coil  by  the  time  varying  magnetic  flux: 


dB 
dt 

- 

V 

nr^ 
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For  example  a  typical  coil  would  consist  of  1 5  turns  of  copper  wire  of  0.8  mm 
diameter  on  a  form  of  50  mm  diameter  (r  =  25  mm)  resulting  in  a  circular  coil 
approximately  12  mm  long.  An  induced  voltage  of  200  mV  would  then  result  in  a 
dB/dt  =  6.79  T/s  coaxial  with  the  search  coil. 

2)  Voltage  measurement  device 

The  device  used  to  measure  voltage  induced  in  the  search  coil  shall  have  a  high 
input  impedance  and  suflBcient  bandv^dth  to  prevent  signal  attenuation,  e.g.  a 
storage  oscilloscope. 

The  voltage  measuring  device  (storage  oscilloscope)  shall  be  placed  at  a  location 
where  it  is  accurate  and  not  affected  by  magnetic  fields. 

The  search  coil  shall  be  connected  to  the  vohage  measurement  device  by  means  of 
a  low  inductance  cable,  e.g.  a  twisted  pair  to  avoid  ringing  on  the  waveform  that 
may  be  experienced  with  coaxial  cables. 

3)  Positing  device 

A  means  shall  be  provided  for  positioning  and  aligning  the  search  coil  in  the 
magnet  in  a  stable  and  reproducible  manner.  The  device  shall  permit  positioning  of 
the  search  coil  throughout  the  region  normally  occupied  by  the  patient  while 
maintaining  the  search  coil  coaxial  with  the  static  magnetic  field. 

c)         Maximum  gradient  rate  measurement  for  trapezoidal  gradients: 

Turn  off  or  maximally  attenuate  the  radio  fi^equency  (RF)  transmitter  to  prevent 
mterference. 

Position  the  search  coil  int  he  Magnetic  Resonance  equipment  using  the  positioning  device 
so  that  the  search  coil  axis  is  aligned  with  that  of  the  static  magnetic  field,  regardless  of 
patient  orientation. 

The  amplitudes,  A^,^  and  A.,^  are  defined  as  the  maximum  gradient  strength  that  can  be 
generated  by  the  systems  under  normal  scan  conditions.  The  transition  time  ti,  tj,  tj  and  ti 
shall  be  the  minimum  that  the  system  can  generate  using  the  standard  pulse  sequence 
generation  software  and  hardware  used  by  the  manufacturer.  Begin  simuhaneous  pulsing 
of  all  three  gradients  (x,  y  and  z)  using  the  maximum  amplitude  bipolar  gradient 
waveform. 

Move  the  search  coil  to  locate  the  position  of  the  maximum  induced  voltage  that  can  be 
found  within  the  region  of  the  Magnetic  Resonance  equipment  normally  occupied  by  a 
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patient.  Once  this  position  is  established,  the  coil  shall  be  immobilized  to  prevent 
movement-induced  voltage  from  interfering  with  the  measurement. 

The  following  procedure  locates  the  approximate  position  of  maximum  dB^dt  in  the 
patient  space  of  the  Magnetic  Resonance  equipment  with  an  axial  field  magnet  and  whole- 
body  gradient  coils.  Consider  a  cartesian  coordinate  system  with  the  Z-axis  pointing  in  the 
direction  of  the  static  field  and  the  Y-axis  pointing  vertically  (e.g.  as  in  figure  103).  For 
the  maximum  gradient  rate  measurement  for  trapezoidal  gradients  procedure,  the  points  of 
maximum  dB/dt  will  occur  where  each  of  the  gradients  is  near  maximum  and  all  the 
gradient  fields  add.  Typically,  this  is  near  the  maximum  vending  density  for  the  z-coil  and 
near  the  wall  of  the  bore  for  X-  and  Y-gradients. 

Assume  that  the  distance,  p,  from  the  magnet  isocentre  to  the  patient  support  along  the 
lines  y  =  ±  X  is  the  same  as  the  maximum  distance  the  patient  may  extend  above  the 
magnet  isocentre  along  the  lines  y  =  ±  x,  then  the  maximum  dB/dt  value  is  located  at  the 
points: 


^2    ^2 


±-^,^,Z. 


m 


Where  Z„  is  the  Z  coordinate  of  maximum  dB/dt. 

Assuming  all  three  gradients  are  applied  simultaneously  and  with  equal  amplitudes,  the 
highest  induced  voltages  will  be  found  near  the  walls  of  the  magnet  shroud  near  the  ends 
of  the  gradient  coils  along  the  lines  y  =  ±  x  (see  figure  103). 

Search  coil  position  1  in  figure  103  shall  be  located  at: 


Y    =    X    =    -P- 
Similarly,  search  coil  position  2  shall  be  located  at: 


Y    =    -X    =    -P- 

>/2 

Place  the  search  coil  coaxial  vsdth  the  static  magnetic  field  at  position  1  and  constrain  the 
search  coil  so  it  may  move  only  in  the  Z  direction.  Move  the  coil  in  the  Z  direction  to  the 
position  where  the  maximum  voltage  is  induced  in  the  search  coil.  Repeat  this  procedure 
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for  position  2.  The  results  from  the  position  which  yields  the  greatest  search  coil  voltage 
shall  be  recorded. 

Measure  the  peak  voltages  which  correspond  to  the  maximum  dB/dt  for  transitions  tj,  t2, 
tj  and  t4  (see  figure  101). 

Record  the  pulse  width,  T,  of  the  dB/dt  waveform  corresponding  to  the  tj,  tj,  t,  and  t4  (see 
figure  101). 
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Figure  101  -  Waveform  for  performing  measurements  of  acoustic  noise  and  dB/dt. 

d)        Maximum  sinusoidal  dB/dt: 

Step  1 :     turn  off  or  maximally  attenuate  the  radio  frequency  (RF)  transmitter  to 
prevent  interference. 

Step  2:     position  the  search  coil  in  the  Magnetic  Resonance  scarmer  using  the 
positioning  device  so  that  the  search  coil  axis  is  aligned  with  that  of  the  static 
magnetic  field. 
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Step  3:     tumoff  all  gradients  other  than  the  sinusoidal  gradient  of  interest.  Apply 
the  maximum  amplitude  and  frequency  waveform  that  the  system  can  generate  using 
the  standard  pulse  sequence  generation  software  and  hardware  used  by  the 
manufacturer. 

Step  4:     move  the  search  coil  to  locate  the  position  of  the  maximum  induced 
vohage  that  can  be  found  within  the  region  of  the  machine  normally  occupied  by  the 
patient.  Once  the  position  is  estabHshed,  the  coil  shall  be  immobilized  to  prevent 
movement-induced  voltage  from  interfering  with  the  measurement. 

Step  5:     for  most  coils  a  suitable  search  procedure  is  as  follows: 

for  X  gradients,  fix  the  coil  along  X  =  p,  Y  =  0 
for  y  gradients,  fix  the  coil  along  X  =  0,  Y  =  p 
for  z  gradients,  fix  the  coil  along  X  =  0,  Y  =  p 

Step  6:     then  search  along  Z  to  locate  the  position  of  the  maximum  induced  voltage 
that  can  be  fond  within  the  region  of  the  Magnetic  Resonance  equipment  normally 
occupied  by  the  patient. 

Step  7:     measure  the  peak  voltage  which  corresponds  to  the  maximum  dB/dt. 

Record  the  pulse  width,  T,  of  the  dB/dt  waveform. 

Repeat  the  above  procedure  for  any  other  sinusoidally  driven  gradient. 

a)     Report  of  results: 

The  following  data  shall  be  recorded  for  the  maximum  gradient  rate  measurement  for 
trapezoidal  gradients: 


Parameter 

Dimension 

-  maximum  gradient  pulse  amplitude 

mT/m 

-  patient  space  radius,  p 

M 

-ti 

JIS 

-t2 

^is 

-t, 

US 

-t4 

\is 

-  measured  voltage 

voltage  indicating  T/s 

The  following  data  shall  be  recorded  for  the  maximum  sinusoidal  dB/dt  method: 
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Parameter  Dimension 

-  maximum  gradient  pulse  amplitude  mT/m 

-  patient  space  radius,  p  M 

-  pulse  width,  T  fis 

-  measured  voltage  voltage  indicating  T/s 
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Protection  against  high  Radio  Frequency  energy 

NOTE:  The  requirements  regarding  SAR  from  RF  power  assumes  that  the  manufacturers 
recommendations  regarding  the  control  of  temperature,  humidity  and  other  environmental 
conditions  are  met.  •       .       ... 

SAR  Limits 

The  following  requirements  on  whole-body  SAR  should  be  used  only  where  the 
environmental  temperature  does  not  exceed  24°C  and  the  relative  humidity  does  not 
exceed  60%.  For  temperatures  above  24°C  or  relative  humidity  above  60%,  the  limits  in 
the  definitions  of  the  operating  modes  for  whole-body  SAR  shall  be  reduced  as  described 
in  the  rationale,  possibility  changing  the  operating  mode  for  a  given  value  of  SAR. 

NOTE:    The  requirements  regarding  SAR  from  RF  power  fiarther  assume  that  the 
SAR  averaged  over  any  10  s  period  does  not  exceed  5  times  the  stated  temporal 
average  SAR  limit. 

a)  Whole-body  SAR 

Whole-body  SAR  is  the  value  of  the  SAR  averaged  over  the  entire  body  of  the  patient. 

Three  operating  modes  with  regard  to  whole-body  SAR  are  defined,  namely  Normal 
Operating  Mode,  First  Level  Controlled  Operating  Mode  and  Second  Level  Controlled 
Operating  Mode: 

1)  the  Normal  Operating  Mode  comprises  values  of  whole-body  SAR  not 
higher  than  1.5  W/kg  averaged  over  any  period  of  15  min: 

2)  the  First  Level  Controlled  Operating  Mode  comprises  values  of  whole- 
body  SAR  not  higher  than  4  W/kg  averaged  over  any  period  of  1 5  min: 

3)  the  Second  Level  Controlled  Operating  Mode  comprises  values  of  whole- 
body  SAR  than  may  exceed  4  W/kg  averaged  over  any  period  of  15  min. 

b)  Head  SAR 

Head  SAR  is  the  value  of  the  SAR  averaged  over  the  head  of  the  patient. 

Two  operating  modes  with  regard  to  head  SAR  are  defined,  namely  Normal  Operating 
Mode  and  Second  Level  Controlled  Operating  Mode: 
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1)  the  Normal  Operating  Mode  comprises  values  of  head  SAR  nor  higher  than 
3  W/kg  averaged  over  any  period  of  10  min: 

2)  the  Second  Level  Controlled  Operating  Mode  comprises  values  of  head 
SAR  that  may  exceed  3  W/kg  averaged  over  any  period  of  10  min. 

c)     Local  tissue  SAR  using  e.g.  transmitting  special  coils 

Local  tissue  SAR  is  the  value  of  the  SAR  averaged  over  any  localized  mass  of  tissue. 

Two  operating  modes  with  regard  to  local  tissue  SAR  are  defined,  namely  Normal 
Operating  Mode  and  Second  Level  Controlled  Operating  Mode: 

1)  the  Normal  Operating  Mode  comprises  values  of  local  tissue  SAR  in  any 
gram  of  tissue  not  exceeding  8  W/kg  in  the  head  or  torso  of  12  W/kg  in  the 
extremities  averaged  over  any  period  of  5  min; 

2)  the  Second  Level  Controlled  Operating  Mode  comprises  values  of  local 
tissue  SAR  that  may  exceed  the  upper  limit  for  the  Normal  Operating  Mode. 

Control  of  SAR 

a)  Transmission  of  body  coil: 

-  when  the  body  coil  is  used  to  image  the  torso,  the  system  shall  control  the  whole- 
body  SAR, 

-  when  the  body  coil  is  used  to  image  the  head,  the  system  shall  control  either  the 
head  SAR  or  the  SAR  averaged  over  the  part  of  the  body  that  is  being  irradiated.  In 
the  latter  case  the  limits  that  apply  are  numerically  equal  to  those  for  the  whole-body 
SAR; 

-  when  the  body  coil  is  used  to  image  extremities,  the  system  shall  control  either  the 
local  tissue  SAR  or  the  SAR  averaged  over  the  part  of  the  body  that  is  being 
irradiated.  In  the  latter  case  the  limits  that  apply  are  numerically  equal  to  those  for 
the  whole-body  SAR. 

b)  Transmission  by  head  coil: 

-  when  the  head  coil  is  used  to  image  the  head,  the  system  shall  control  the  head 
SAR; 
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-  when  the  head  coil  is  used  to  image  extremities,  the  system  shall  control  either  the 
local  tissue  SAR  or  the  SAR  averaged  over  the  part  of  the  body  that  is  being 
irradiated.  In  the  latter  case  the  limits  that  apply  are  numerically  equal  to  those  for 
the  head  SAR. 

c)      Other  combinations  of  coil  and  body  areas: 


-  regardless  of  the  coil  used  for  RF  excitation,  the  system  shall  control  whole  body 
SAR,  and  local  tissue  SAR  appropriately. 

Display  of  SAR 

The  value  of  the  SAR  that  is  displayed  shall  be  the  value  that  is  controlled  by  the  system. 

Compliance  shall  be  tested  by  measurement  of  the  SAR  according  to  one  of  the  methods 
specified  below. 

Method  of  measurement  of  SAR  to  demonstrate  compUance. 

The  methods  employed  may  be  known  as  "pulse  energy  method"  and  "calorimetric 
method". 


a)     Auxiliary  definitions: 
tip  angle: 
landmark: 


fbrwanl2° 


reflected  2- 


other  2' 


object- 
cod- 


'  fivwnd  1* 


Angle  through  which  macroscopic  magnetization  vector  is 

moved  by  a  radio  firequency  pulse 

Region  of  the  body  to  be  placed  in  the  center  of  the  RF 

body  coil 

Forward  peak  RF  power  going  into  the  terminals  of  the  RF 

coil  when  loaded  with  test  object  2 

Peak  RF  power  reflected  fi^om  the  terminals  of  the  RF  coil 

when  loaded  with  test  object  2 

Peak  RF  power  absorbed  or  dissipated  elsewhere  in  the 

system  when  the  coil  is  loaded  with  test  object  2.  For 

example,  in  a  quadrature  system  fed  by  a  four-port  splitter, 

this  is  the  peak  RF  power  delivered  to  the  dummy  load. 

Net  peak  RF  power  delivered  to  test  object  2 

Net  peak  RF  power  absorbed  by  the  RF  transmit  coil  for  the 

tip  angle  specified 

Forward  peak  RF  power  going  into  the  terminals  of  the  RF 

coil  when  loaded  with  test  object  1 
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Prekcted  1^  Peak  radio  RF  reflected  from  the  terminals  of  the  RF  coil 

when  loaded  with  test  object  1 

Pofter  1^  Peak  RF  power  absorbed  or  dissipated  elsewhere  in  the 

system  when  the  RF  coil  is  nearly  unloaded.  For  example, 
in  a  quadrature  system  fed  by  a  four-port  splitte,  this  is  the 
peak  power  delivered  to  the  dummy  load  when  the  radio 
frequency  coil  is  loaded  with  test  object  1 

B 1 :  Magnetic  induction  of  the  RF  magnetic  field  used  in 

Magnetic  Resonance  imagmg  and  spectroscopy 

b)  Test  hardware 

Test  object  1 

Test  object  1  shall  be  a  physically  small  object  used  during  the  calibration  of  tip  angles  and 
during  the  determination  of  radio  frequency  coil  losses  with  the  pulse-energy  method. 
Test  object  1  shall  be  filled  with  a  material  which  can  be  imaged  by  the  Magnetic 
Resonance  equipment,  but  whole  conductivity  shall  be  less  than  0.003  S/m.  This  ensures 
that  the  radio  frequency  losses  of  the  object  are  negUgible,  thus  permitting  the 
measurement  of  coil  losses  while  contributing  little  to  system  losses.  Test  object  1  shall  be 
physically  small  and  compact  in  shape.  The  test  object  volume  shall  be  less  than  250  ml. 
The  longest  dimension  of  the  test  object  shall  be  not  more  than  three  times  the  shortest 
dimension. 

Test  object  2 

Test  object  2  shall  be  the  test  object  whole  SAR  shall  be  determined  using  the  pulse- 
energy  method.  This  test  object  should  have  coil  loading  characteristics  similar  to  patient 
loading.  These  loading  characteristics  will  differ  for  the  case  of  whole-body  SAR,  head 
SAR  and  local  tissue  SAR. 

c)  Test  arrangements 

Peak  radio  frequency  power  measurements  shall  be  made  either  with  the  use  of  an 
oscilloscope  whole  bandwidth  exceeds  the  desired  frequency  by  a  factor  of  5  or  more  and 
dual  directional  couplers  whole  directivity  exceeds  30  dB,  or  with  attenuators. 

Connect  the  directional  coupler,  oscilloscope  and  other  power  measurement  equipment 
with  50  Q  coaxial  transmissions  lines  as  shown  in  figure  105  or  figure  106.  Ensure  that 
the  oscilloscope  input  is  terminated  in  a  50  Q  impedance.  The  oscilloscope  shall  be  placed 
at  a  location  where  it  is  accurate  and  not  affected  by  magnetic  fields. 


100 


RF 

power 

amplifier 


Oscilloscope 


Dual 

directional 

coupler 


Pforward 


Preflected 


Transmission  line 


Dummy  load  port 


Pother 


Quadrature 
RF  splitter 


RF 
coil 


NOTE  -  Not  shown  are  the  50  H 
terminations  on  the  oscilloscope. 


Figure  105  -  Hardware  setup  for  pulse-energy  method  for  the  measuremtn  of  SAR 
with  a  quadrature  coil. 


101 


RF 

power 

amplifier 


Dual 

directional 

coupler 


Oscilloscope 


Pfoiward 


Preflected 


NOTE  -  Not  shown  are  the  50  fl 
terminations  on  the  oscilloscope. 


Figure  106  -  Hardware  setup  for  pulse-energy  method  for  the  measurement 
of  SAR  with  a  linear  coil. 


d)     SAR  measurement  procedure,  pulse-energy  method 

Step  1 :     position  test  object  1  in  the  scanner  at  isocentre. 

Step  2:     calibrate  the  tip  angle.  The  standard  method  (hardware  and/or  software) 
employed  by  each  manufacturer  shall  be  used  to  set  the  tip  angle  to  the  desired  value 
for  the  type  of  scan  selected.  The  same  tip  angle  shall  be  used  for  test  object  1  and 
test  object  2  power  measurements. 

Step  3:     start  the  scan  sequence  chosen  for  the  SAR  measurement.  Ensure  that  the 
scan  time  is  long  enough  to  allow  the  completion  of  the  desired  measurements,  or  set 
the  scanner  to  repeat  the  pulse  sequence  indefinitely. 

Step  4:     measure  the  peak-to-peak  radio  fi^equency  envelope  voltages  (forward, 
reflected  and  other)  for  each  radio  frequency  pulse  in  the  sequence,  or  use  other 
means  to  measure  the  instantaneous,  radio  frequency  peak  power  levels  (forward, 
reflected  and  other  if  applicable). 
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Step  5:  measure  the  power  attenuation  of  the  transmission  line  between  the  body- 
coil  and  the  coupler  and  the  transmission  line  between  the  dummy  load  port  and  the 
oscilloscope. 

Step  6:  convert  the  peak-to-peak  voltage  measurements  into  peak  measurements. 
The  peak  power  measured  at  the  oscilloscope  is: 


PP 


8Z 

o 


where 

Vpp  is  the  measured  peak-to-peak  voltage: 

Zq  is  the  characteristic  impedance  of  the  transmission  line  (and  also  the 

input  impedance  of  the  oscilloscope). 

The  calculated  peak  power  at  the  oscilloscope  shall  be  converted  into  peak  power 
referenced  to  the  input  terminals  of  the  radio  frequency  coil  by  considering  transmission 
Une  attenuation  and  directional  coupler  coupling  constants. 

For  a  quadrature  radio  frequency  coil: 

Pfotwaid  ~  Pfin  ^  '"l 
'^reflected  ~  ^m  '  ^l 
Pother  ~  Pom  '  **2 

where 

Pfi„       is  the  forward  power  measured  at  the  terminals  of  the  oscilloscope 

coupler; 
P^       is  the  reflected  power  measured  at  the  terminals  of  the  oscilloscope; 
Po„       is  the  "other"  power  measured  at  the  terminals  of  the  oscilloscope, 
aj         is  the  fractional  power  attenuation  (not  in  dB)  of  the  transmission 

line  between  the  body  coil  and  the  coupler;  and 
Oj        is  the  attenuation  between  the  dummy  load  part  and  the 

oscilloscope. 

Note  that  for  unmatched  linear  coils,  transmission  line  attenuation  increases  with  the 
mismatch: 

a=  a„X(R'+l)/(2XR) 
where 
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a         is  the  resultant  attenuation; 

cLq        is  the  attenuation  under  matched  conditions; 

R         is  the  standing  wave  ratio. 

Calculate  the  peak  power  absorbed  by  the  unloaded  coil,  P^. 
P     =  P  -  P  -  P 

^  cofl       ^ forward  1      ^reflected!      '^ other  1 

Note  that  Pfon^anj  1.  Preflected  1  ^nd  P^j^  j  are  from  the  experiment  using  test  object  1. 

After  initial  measurements  are  made  on  test  object  1,  they  are  repeated  using  test  object  2. 

Remove  test  object  1  from  the  coil,  place  test  object  2  in  the  scanner  at  isocentre  and 
repeat  steps  4  and  6.  The  same  tip  angle  shall  be  used  for  test  object  2  as  was  used  for 
test  object  1. 

Calculate  the  peak  power  absorbed  in  test  object  2,  Pobject- 
P       =P  -P  -P         p 

'^  objecl       'forward  2      '^reflected  2      •'^oflier2-^cofl 

Note  that  Pfon^arfj,  Preflected2  ^nd  pofl,„2  ^Tc  all  obtained  from  the  measurements  on  test  object 
2,  while  Pjoa  is  obtained  from  calculations  usmg  test  object  1. 

Calculate  the  energy  per  pulse  absorbed  by  the  object:  integrate  each  of  the  radio 
frequency  waveforms  in  the  pulse  sequence  and  find  the  pulse  durations  of  equivalent  (in 
energy)  rectangular  pulses  which  have  the  same  peak  power  levels  as  the  pulse  sequence 
being  evaluated.  This  step  shall  be  repeated  for  each  type  of  radio  frequency  pulse  in  the 
pulse  sequence  (e.g.  90°  pulses  as  well  as  180°  pulses).  The  equivalent  rectangular  pulses 
shall  resuh  in  both  the  same  peak  power  and  in  the  same  energy  deposition  as  the  pulse 
being  evaluated.  The  pulse  duration,  W^,^,^  shall  be  calculated  from  the  radio  frequency 
pulse  waveform,  Bl(t)  normalized  to  its  peak  value  Blp^  as  follows: 

For  example,  a  sine  (or  sin0/6)  pulse  is  approximately  equivalent  in  energy  to  a 
rectangular  pulse  vwth  the  same  peak  power  and  half  the  width  of  the  main  lobe  of  the  sine 
pulse. 

Calculate  the  energy  per  pulse,  Ej ,  as: 

^  ~  "object  ^   "rectangle 


.  *»  •  -v 
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Calculate  the  total  average  power,  P^,  absorbed  by  the  object.  The  total  average  power 
for  sequences  with  well-defined  pulse  repetition  times  (TR)  is  the  sum  of  all  the  energies  in 
all  the  pulses  per  repetition  divided  by  TR; 


p     =  E 

ave  i-^ 


E^ 


I 


.  TR> 


For  sequences  where  TR  is  undefined  (such  as  single-shot  sequences),  the  total  average 
power  shall  be  calculated  fi-om  the  sum  of  all  energies  in  all  pulses  used  in  the  scan  divided 
by  the  total  scan  time,  S: 


p     =  E 


'  E^ 


Calculate  the  SAR  for  test  object  2,  whose  mass  is  M: 


TAS    = 


ave 


M 


e)     Report  of  results,  pulse  energy  method 

The  following  data  shall  be  recorded: 

Measured  values  of  SAR  for  each  object  and  pulse  sequence  used 

Parameter 
Dimension 

-  Peak  forward  power 
W 

-  Peak  reflected  power  W 

-  Peak  coil  loss  power  W 

-  Other  peak  power  lost  W 

-  Radio  fi-equency  waveform 

-  Pulse  parameters 

5Tipangle(s) 

-  Accuracy  of  power  measurement  equipment  % 
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-  Equivalent  rectangular  pulse  duration  (W^,3^J  ^s 
-TRorS 

ms 

-  Mass  of  test  object  2  kg 

-  Any  other  parameters  required  to  ensure  repeatability 


f)      Test  hardware,  calorimetry  method  (whole-body  SAR  only) 

Test  object  1 

Test  object  1  shall  be  a  physically  small  object  used  during  the  calibration  of  tip  angles  and 
during  the  determination  of  radio  frequency  coil  losses  with  the  pulse-energy  method. 
Test  object  1  shall  be  filled  with  a  material  which  can  be  imaged  by  the  Magnetic 
Resonance  equipment,  but  whole  conductivity  shall  be  less  than  0.003  S/m.  This  ensures 
that  the  radio  frequency  losses  of  the  object  are  negligible,  thus  permitting  the 
measurement  of  coil  losses  while  contributing  little  to  system  losses.  Test  object  1  shall  be 
physically  small  and  compact  in  shape.  The  test  object  volume  shall  be  less  than  250  ml. 
The  longest  dimension  of  the  test  object  shall  be  not  more  than  three  times  the  shortest 
dimension. 

Test  object  3 

Test  object  3  shall  be  an  annular  test  object  designed  to  provide  radio  frequency  coil 
loading  equivalent  to  that  of  a  human.  It  shall  be  filled  with  an  dqueous  sodium  chloride 
solution  whose  concentration  is  adjusted  to  yield  a  conductivity  which  produces  the 
loading  specified  below.  The  filling  solution  shall  also  contain  a  doping  material  such  as 
manganese  chloride  to  shorten  the  relaxation  properties  of  the  solution  to  render  it 
invisible  in  Magnetic  Resonance  images.  The  total  volume  of  the  test  object  shall  be  less 
than  25  L.  The  test  object  shall  be  designed  to  load  the  same  as  a  50  kg  to  90  kg  human, 
positioned  for  a  scan  with  the  stated  tolerances: 

Coil  3  dB  bandwidth  ±  15% 

Coil  impedance:  ±  20%  phase 

Coil  center  frequency  shift:  ±  1%  of  centre  frequency 

Annular  test  object  thermal  properties  shall  be  suflBcient  to  maintain  a  2  "C  temperature 
rise  over  ambient  temperature  (in  the  filler  material)  to  within  5%  for  1  h. 

Temperature  measurement 

A  temperature  measurement  system  accurate  to  within  0. 1  °C  over  the  range  of  1 5  °C  to 
45  °C  shall  be  provided.  The  system  will  be  used  to  make  initial  and  final  measurements  of 
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temperature.  It  is  not  necessary  for  the  system  to  be  immune  to  radio  frequency 
interference  from  a  Magnetic  Resonance  equipment  nor  to  temperature  artifacts  induced 
by  the  static  magnetic  field.  Acceptable  examples  are  fiber-optic  temperature  sensors, 
thermister-based  temperature  sensors  and  the  thermocouple  temperature  sensors. 

g)     Test  arrangement,  calorimetry  method 

Place  the  temperature  measurement  hardware  in  a  convenient  location  remote  from  the 
magnet  to  permit  measurement  of  test  object  filler  temperature  quickly  before  and  after 
scanning.  This  will  prevent  inaccuracies  due  to  the  static  magnetic  field. 

Position  test  object  3  at  isocentre  in  the  radio  frequency  body  coil  before  the  scan. 

The  number  of  acquisitions  for  the  scan  shall  be  adjusted  to  yield  a  scan  suflBciently  long 
to  create  a  temperature  rise  at  least  20  times  greater  than  the  error  of  the  temperature 
measurement  system.  For  an  aqueous  solution,  an  exposure  of  one  hour  to  an  SAR  of  1 
W/kg  will  resuh  in  an  average  temperature  rise  of  0.86  "C  (0.0143  °C/min),  assuming  no 
thermal  losses.  Higher  temperature  rises  may  be  measured  to  more  significant  places  than 
small  temperature  rises.  However,  higher  temperatures  rises  may  increase  heat  loss  errors. 

h)     SAR  measurement  procedure,  calorimetry  method 

Step  1 :     equihbrium  test  object  3  in  the  magnet  bore  with  air  flow  turned  off. 

Step  2:     measure  the  initial  temperature  of  the  bore  of  the  Magnetic  Resonance 
equipment  and  the  initial  temperature  of  the  scan  room. 

Step  3 :     remove  the  test  object  from  the  bore  and  invert  three  times  to  stir  the  filler 
material. 

Step  4:     measure  the  initial  temperature,  Tj,  of  the  test  object  by  inserting  a 
temperature  sensor  into  the  filler  material.  The  initial  bore  temperature  should  not 
differ  from  Tj  by  more  than  1  °C. 

Step  5:     quickly  replace  the  test  object  3  in  the  Magnetic  Resonance  scanner  at 
isocentre  and  place  test  object  1  centrally  within  it. 

Step  6:     calibrate  the  tip  angle.  The  standard  methods  (hardware  and/or  soflware) 
employed  by  each  manufacturer  shall  be  used  to  set  the  tip  angle  to  the  desired  value 
for  the  type  of  scan  selected. 

Step  7:     scan  with  enough  acquisitions  to  ensure  the  desired  measurement  (see  item 
g  above). 
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Remove  the  test  object  from  the  bore  and  measure  the  final  temperature,  Tj,  of  the  filler 
solution  of  test  object  3  after  inverting  the  test  object  three  times. 

If  Tf  -  Tj  <  2  °C,  the  experiment  can  be  continued  by  replacing  the  test  object  in  the  bore 
and  restarting  the  scan.  The  measurement  delay  time  shall  be  subtracted  from  the  total 
time  to  assure  accuracy. 

Measurement  the  final  bore  and  the  final  scan  room  temperatures.  If  the  scan  room 
temperatures  of  bore  temperature  changes  by  more  than  0.5  "C/h  then  the  experiment 
should  be  repeated  under  more  stable  conditions. 

Calculate  the  energy,  E  (in  J),  absorbed  by  the  test  object  3  in  terms  of  the  mass,  M,  of  the 
filler  (in  kg)  and  of  the  specific  heat,  c,  of  the  filler  fluid  (in  J/kg  °C): 

E  =  Mxcx(Tf-Ti) 
Calculate  the  average  power  P  (in  W)  during  the  scan  time  T  (in  s): 

P  =  E/t 
Calculate  the  SAR  of  the  test  object  filler  material  from: 

SAR  =  P/M 
Calculate  the  apparent  SAR,  SAR^pp  using: 

where  M^,,  is  the  mass  of  a  body  which  test  object  3  simulates  (assuming  that  the 
validity  of  the  simulation  has  been  independently  verified). 

Report  of  results,  calorimetry  method 

Parameter  Dimension 

-  Initial  test  object  3  temperature  °C 

-  Final  test  object  3  temperature  °C 

-  Initial  bore  temperature  °C 

-  Final  bore  temperature  °C 

-  Initial  scan  room  temperature  °C 

-  Final  scan  room  temperature  *^C 

-  Test  object  3  mass  kg 

-  Test  object  3  simulated  mass  (M^  kg 
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Test  object  3  volume  1 

Test  object  3  dimensions  m 
Radio  frequency  waveform  description 

Repetition  time  (TR)  ms 

Echo  time  (TE)  ms 

Number  of  slices  — 

Number  of  echoes  — 

Total  scan  time  s 

Magnetic  Resonance  tip  angle  employed  in  study  degree 
Any  other  parameters  required  to  ensure  repeatability 
The  accuracy,  type  of  location  (during  the  measurement) 
of  the  temperature  measurement 
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Definitions: 

Scan  Duration  -  This  term  is  often  referred  to  as  acquisition  time  and  is  a  variable  needed 
in  algorithms  that  predict  the  energy  exposure  to  radio  fi-equency  power.  Scan  duration 
for  acquisition  gated  to  period  functions  of  the  patient  may  be  diflBcult  to  assess 
accurately.  A  conservative  calculation  of  scan  duration  in  such  cases  can  be  assumed  to 
be  based  on  the  shortest  period  that  can  be  realized  under  the  circumstances. 

Magnetic  Resonance  Examination  -  A  magnetic  resonance  examination  usually  is 
comprised  of  more  than  one  acquisition  of  data,  each  with  a  predictable  scan  duration 
during  one  session. 

Magnetic  Resonance  -  The  phenomenon  of  magnetic  resonance  occurs  when  the 
fi-equency  of  electro-magnetic  radiation  equals  the  Larmor  precession  fi-equency  of  the 
nuclear  or  electron  magnetic  moments. 

Specific  Absorption  Rate  (SAR)  -  The  SAR  is  a  fianction  of  the  fi-equency  (increasing 
approximately  as  the  square  of  the  fi-equency),  the  type  and  number  of  radio  fi-equency 
pulses,  the  duration  and  repetition  rate  of  pulses  and  the  type  of  coil  used  for  transmission. 
The  important  biological  factors  are:  conductivity  of  tissue,  specific  gravity  of  the  tissue, 
anatomical  region  examined,  tissue  type  (e.g.  the  degree  of  perfusion)  and  mass  of  the 
patient. 

Time  Rate  of  Change  of  Magnetic  Field  (dB/dt)  -  The  rate  of  change  of  the  magnetic  flux 
density  with  time  is  a  potential  concern  for  safety  limits,  since  changing  magnetic  fields 
induce  eddy  currents  within  the  patient. 

Quench  -  a  quench  is  due  to  excessive  heating  of  the  wires  of  the  magnet  immersed  in 
liquid  helium. 

Superconducting  magnet  can  vent  up  to  several  hundred  liters  of  cryogenic  gases  per 
hour  during  normal  operation.  During  the  quench,  approximately  l**  L  to  l'**  L  of 
gas  at  atmospheric  pressure  can  be  vented  within  a  few  minutes. 

Usually,  a  quench  occurs  when  the  quantity  of  liquid  helium  becomes  insufficient  to 
cool  the  superconducting  coil.  Due  to  the  increase  in  the  temperature  of  the  coil,  the 
superconducting  wire  exhibits  normal  conductivity  and  an  excessive  boil-off  starts. 

If  proper  venting  is  not  used,  three  effects  can  occur  during  rapid  boil-off  at  a 
quench: 

-  excessively  cooled  gases  will  freeze  water  molecules  in  the  area  adjacent  to  the 
magneti  causing  a  dense  white  fog. 
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-  oxygen  in  the  room  will  be  displaced  by  helium  making  it  difficuh,  if  not 
impossible,  to  breathe; 

-  the  helium  gas  that  escapes  during  a  quench  will  freeze  any  object  it  contacts. 

Electromagnetic  compatibility 

Within  the  controlled  access  area  the  magnetic  resonance  system  will  generally  not  comply 
with  the  current  requirements  for  radio  frequency  emission.  These  requirements  are 
primarily  intended  to  protect  radio  communication  and  are  laid  down  in  International 
Standards  such  as  CISPR  1 1 .  Permissible  limits  for  radio  frequency  emission  are  in  the 
range  of  30  dBjiV/m  to  50  dB^V/m.  It  is  proposed  that  in  coming  standards  emission 
levels  of  magnetic  resonance,  computed  tomography  and  large  x-ray  systems  are  to  be 
measured  at  the  boundaries  of  the  room  or  of  the  building. 

EEC  standards  on  EMC  are  in  preparation  that  deal  with  the  immunity  for  radio  frequency 
fields  of  medical  equipment.  It  is  expected  that  immunity  will  be  required  for  the  general 
case  of  fields  up  to  1  V/m  or  3  V/m  and  for  special  cases,  such  as  life  sustaining  and  some 
patient  monitoring  equipment,  at  fields  up  to  10  V/m  or  100  V/m. 

Actually,  within  the  controlled  access  area  around  the  magnetic  resonance  equipment, 
radio  frequency  field  strengths  may  easily  exceed  these  limits  and  can  have  values  in 
excess  of  100  V/m.  It  is  clear,  that  peripheral  equipment  used  in  the  controlled  access 
area  might  be  disturbed  by  this  field. 

Protection  against  high  dB/dt  associated  with  gradient  coils 

Time-varying  magnetic  fields  induce  electrical  potentials  and  eddy  currents  in  accord  with 
Faraday's  law.  In  the  simple  case  of  a  circular  loop  intersected  by  a  uniform  time- varying 
magnetic  field,  the  magnitude  of  the  average  induced  electric  field  depends  directly  on  the 
radius  of  the  loop  and  the  magnitude  of  dB/dt.  The  larger  the  radius  and  dB/dt,  the  larger 
is  the  induced  current. 

The  dB/dt  waveform  may  be  pulsed  or  sinusoidal.  The  pulse  width,  T,  is  defined  as  the 
fiill  width  at  half  maximum  of  a  single  pulse  of  dB/dt.  It  should  be  noted,  that  "pulse"  and 
"waveform"  refer  to  the  time  fimction  of  dB/dt  and  not  to  the  current  pulse  in  the  gradient 
coil.  When  the  waveform  is  sinusoidal  rather  than  pulsed,  T  is  the  half  period  of  the 
sinusoid. 

When  a  time-varying  magnetic  field  impinges  on  the  human  body,  currents  are  induced  in 
the  body  in  planes  perpendicular  to  the  magnetic  field.  These  induced  currents  can  affect 
excitable  tissue  with  thresholds  that  are  frequency  dependent.  At  higher  frequencies, 
higher  currents  are  required  for  an  effect.  Although  induced  currents  will  also  cause 
heating  according  to  Ohm's  law,  heating  effects  become  the  limiting  consideration  only  for 
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frequencies  above  200  kHz,  or  pulse  widths  shorter  than  2.5  ^s.  Heating  effects  are 
covered  in  the  sections  on  radio  frequency  energy. 

The  primary  concerns  regarding  dB/dt  are  cardiac  fibrillation  and  the  stimulation  of 
cardiac  and  peripheral  nerve  and  muscle  (twitching),  fibrillation  being  most  serious  in 
nature  and  induced  by  a  larger  dB/dt  value  and  peripheral  muscle  stimulation  being  of  least 
concern  and  induced  by  a  smaller  dB/dt  value.  Protection  against  peripheral  nerve 
stimulation  should  serve  to  protect  also  against  cardiac  stimulation. 

Most  empuical  data  on  thresholds  for  cardiac  and  neuromuscular  responses  are  based  on 
animal  studies  or  a  limited  number  of  observations  on  humans  during  appUcation  of  50  Hz 
or  60  Hz  sinusoidal  currents  using  electrodes.  It  has  been  reported  that  in  animal  studies 
time-varying  magnetic  fields  induce  skeletal  muscle  contractions  at  dB/dt  values  at  which 
no  disturbances  in  cardiac  fijnction  occur. 

Reports  describe  peripheral  nerve  and  phrenic  nerve  stimulation  at  dB/dt  values  and 
frequencies  that  are  similar  to  those  predicted  are  described  below.  A  more  recent  report 
describes  the  threshold  dB/dt  for  localized  excitation  of  the  brachial  plexus  in  human 
subjects  using  a  train  of  trapezoidal  pulsed  magnetic  fields  for  read-out  gradient  pulses  of 
"echo-planar"  imaging  systems.  These  later  resuUs  are  also  consistent  with  the 
predications  outlined. 

Cardiac  fibrillation  thresholds  from  animal  experiments  where  electrodes  were  placed 
directly  on  the  heart  show  a  log-normal  distribution  with  the  one-percentile  value  at  about 
half  the  medium  value.  Combining  this  relationship  with  human  data  which  showed  an 
average  fibrillation  threshold  of  about  5  A/m^  the  one-percentOe  level  for  fibrillation 
would  be  about  2.5  A/m^  At  60  Hz  the  cardiac  excitation  thresholds  for  humans  are 
about  40%  of  the  60  Hz  fibrillation  values,  and  this  would  indicate  a  threshold  for 
excitation  of  about  1  A/m^. 

A  model  study  compared  the  expected  thresholds  for  cardiac  and  peripheral  nerve 
stimulation  as  a  fiinction  oft.  For  x  shorter  than  0.4  ms  (frequency  greater  than  1200 
Hz),  the  threshold  for  peripheral  nerve  stimulation  is  at  least  10  times  lower  than  for 
cardiac  stimulation.  For  x  longer  than  5  ms  (frequency  below  100  Hz),  however,  the 
modeled  threshold  for  peripheral  nerve  stimulation  (about  60  T/s  for  the  main  magnetic 
field  parallel  to  the  body)  is  only  about  15%  below  that  for  cardiac  stimulation  (about  70 
T/s). 

"Echo-planar"  imaging  techniques  use  sinusoidal  excitation  at  about  1200  Hz  to  1500  Hz 
on  one  or  more  of  the  gradient  coils.  For  frequencies  above  1200  Hz,  peripheral  nerve 
stimulation  should  have  a  10  times  lower  value  than  that  for  cardiac  stimulation  threshold 
and  supra-threshold  peripheral  stimulation  should  provide  an  early  warning  for  cardiac 
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stimulation  with  an  adequate  margin  of  safety.  Cardiac  fibrillation  has  a  threshold  5  to  10 
times  higher  than  cardiac  stimulation  at  these  frequencies. 

Extensive  calculations  suggest  that  peripheral  nerve  and  muscle  (including  cardiac  muscle) 
stimulation  will  be  avoided  if  induced  current  density  is  restricted  to  less  than  400  mA/m^ 
at  any  body  location  for  induced  current  pulse  widths  exceeding  120  |as.  For  shorter 
current  pulse  widths,  the  induced  current  density  should  be  restricted  according  to  the 
relationship  J  x  t  <  48  x  10^,  where  J  is  the  value  of  the  induced  current  density  expressed 
in  mA/m^  and  t  is  the  duration  of  the  induced  current  pulse  in  jis.  When  only  cardiac 
stimulation  is  to  be  avoided,  i.e.  when  the  risk  of  peripheral  nerve  stimulation  is  not  taken 
into  account,  relaxation  of  this  restriction  is  allowed  for  periods  of  magnetic  field  change 
of  less  than  3  ms  according  to  the  relationship  J  x  t  <  1.2  x  10*. 

Making  certain  assumptions,  limits  expressed  in  dB/dt  can  be  derived  fi-om  these 
restrictions  on  induced  current  density.  It  can  be  calculated  that  peripheral  nerve  and 
cardiac  stimulation  will  be  avoided  by  an  adequate  margin  of  safety  if  the  time- varying 
magnetic  field  does  not  exceed  20  T/s  peak  value  for  x  greater  than  120  us.  As  the  value 
of  T  decreases  between  120  ^s  and  2.5  \is,  the  safe  level  increases  linearly  to  a  value  of 
960  T/s,  above  which  heating  fi-om  eddy  currents  will  dominate  (greater  than  1  W/kg),  and 
the  radio  fi-equency  SAR  limits  are  more  appropriate  limiting  factors  than  dB/dt. 

Recent  reports  cite  peripheral  and  cardiac  nerve  stimulation  at  dB/dt  values  and 
fi-equencies  which  are  remarkably  similar  to  those  predicted  and  lend  fiirther  credence  to 
the  dB/dt  limits  stated  in  5 1 . 1 02. 

The  values  suggested  in  5 1 .  102  are  intended  to  avoid  cardiac  stimulation  by  an  adequate 
margin  of  safety.  Operation  of  the  system  beyond  those  levels  entails  the  users  and/or 
operators  direct  responsibility. 

For  t  greater  than  5  ms  (fi-equency  less  than  100  Hz),  cardiac  stimulation  thresholds  are 
only  slightly  higher  than  peripheral  nerve  stimulation  thresholds.  In  this  range  peripheral 
nerve  stimulation  thresholds.  In  this  range  peripheral  nerve  stimulation  should  not  be  used 
as  a  surrogate  or  "early  warning"  for  cardiac  stimulation.  However  it  is  not  expected  that 
Magnetic  Resonance  equipment  will  operate  with  such  fi-equencies  or  pulse  widths. 

The  derived  restrictions  on  dB/dt  can  be  exceeded  if  it  can  be  shown  that  the  basic 
restriction  on  induced  current  density  is  not  exceeded. 

Protection  against  high  radio  frequency  energy 

Heating  is  a  major  consequence  of  exposure  to  the  RF  fi-equencies  used  in  Magnetic 
Resonance  (usually  greater  than  1  MHz).  Many  of  the  biological  effects  of  acute  exposure 
to  RF  are  consistent  with  responses  to  induced  heating  resuUing  fi-om  rises  in  tissue  or 
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body  temperature  of  about  TC  or  more,  or  with  responses  for  minimizing  the  total  heat 
load. 

The  RF  induced  heat  load  can  be  directly  related  to  the  SAE  or  to  the  SAR.  Other 
important  factors  which  influence  patient  response  to  a  given  heat  load  include  the  air 
temperature,  relative  humidity,  air  flow  rate  and  the  degree  of  patient  insulation. 

Localized  regions  of  heating  or  "hot  spots"  may  cause  local  temperature  rises.  In  practice, 
the  local  tissue  SAE  is  always  higher  than  the  average  whole-body  value:  only  part  of  the 
body  is  exposed  to  RF,  power  absorption  increases  in  proportion  to  the  square  of  the 
radial  distance  from  the  center  of  the  body  and  the  electrical  inhomogenety  of  the  body 
alters  current  flow  and  local  energy  absorption.  Studies  using  spherical  models  predict 
that  worst  case  "hot  spots"  would  be  produced  by  a  sphere  of  low  conductivity  (bone  or 
fat)  located  on  the  outer  edge  of  a  large  conductive  sphere  (e.g.  muscle);  a  "hot  spot"  can 
occur  with  an  SAR  of  up  to  2.5  times  the  local  average  value.  Recent  calculations  based 
on  a  heterogeneous  mathematical  model  of  the  human  body  suggest  that  localized  tissue 
SAR  within  the  body  could  be  up  to  10  or  even  70  times  greater  over  small  volumes 
compared  to  the  whole  body  average  value.  However,  these  relative  increases  are  reduced 
to  a  factor  of  about  2  to  4  when  averaged  over  complete  individual  body  organs.  The 
probability  of  high  localized  tissue  SAR  is  reduced  by  quadrature  excitation  and  its  effects 
are  moderated  by  thermal  diffusion  and  blood  flow. 

Thermal  effects  arise  because  of  the  temperature  sensitivity  of  most  biological  processes. 
The  primary  concerns  regarding  RF  exposure  are  to  avoid  excessive  physiological  •  ^    '•  •  ; 
responses  to  a  marked  rise  in  body  temperature  and  to  avoid  raising  tissue  temperature  to 
a  level  which  incur  some  degree  of  harm. 

The  most  sensitive  responses  of  humans  to  acute  whole-body  heating  by  radio  frequency 
are  probably  those  concerned  with  thermoregulation  and  include  increased  cardiac  output 
and  skin  blood  flow  coupled  with  a  slight  drop  in  arterial  blood  pressure.  These  responses 
become  maximal,  even  in  subjects  lying  passively  in  ambient  temperatures,  as  body 
temperatures  rise  by  more  than  about  2T.  Adverse  health  effects  are  not  expected  in 
people  with  unimpaired  thermoregulatory  and  cardiovascular  functions  if  the  increase  in 
body  temperature  does  not  exceed  \°C;  some  specific  exceptions  are  given  below. 

The  health  of  some  patients  might  be  comprised  under  these  conditions.  Physiological  and 
physical  correlates  which  are  associated  vAth  a  decreased  ability  to  adapt  to  an  increased 
heat  load  include  old  age,  obesity  and  hypertension.  Various  drugs  such  as  diuretics, 
tranquillizers  and  sedatives,  vasodilators  and  others  decrease  heat  tolerance.  In  addition, 
the  thermoregulatory  ability  of  infants  is  not  well  developed;  pregnant  women  may  also  be 
comprised  in  their  ability  to  dissipate  heat.  In  this  context,  it  is  worth  noting  that  heat  loss 
from  the  embryo  and  fetus  across  the  placental  barrier  may  be  less  efficient  than  heat 
dissipation  in  other  well  vascularized  tissue.  Raised  body  temperature  is  known  to  be 
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teratogenic  to  a  number  of  mammalian  species  including  primates,  and  has  been  implicated 
in  central  nervous  system  and  facial  defects  in  children  hose  mothers  developed  prolonged 
severe  hyperthermia  (>  39°C),  especially  during  the  first  trimester  of  pregnancy.  In  these 
cases  it  is  desirable  to  limit  rises  in  body  temperature  to  less  than  0.5  °C  (IRPA/INIRC 
1991) 

The  whole-body  thermoregulatory  response  of  humans  exposed  to  Magnetic  Resonance 
imaging  has  been  mathematically  modeled.  The  body  temperature  of  a  lightly  clothed 
patient  whose  thermoregulatory  ability  was  unimpaired  would  rise  by  up  to  0.6  °C, 
depending  on  environmental  conditions,  during  exposure  at  whole-body  SARs  of  up  to  4 
W/kg.  The  results  of  these  calculations  are  in  reasonable  agreement  with  studies  of 
volunteers  exposed  to  up  to  4  W/kg  for  20  min  to  30  min.  In  addition,  there  have  been  no 
reports  of  adverse  health  effects  fi-om  excessive  local  or  systematic  heating  in  patients 
exposed  at  whole-body  SARs  of  up  to  about  1.5  W/kg. 

Limits  on  whole-body  S  AR  are  presented  in  5 1 . 1 03 .  A  limit  of  1 . 5  W/kg  for  the  Normal 
Operating  Mode  is  recommended  as  the  highest  whole-body  SAR  which  all  persons, 
regardless  of  health  status,  should  be  able  to  tolerate.  People  with  unimpaired 
thermoregulatory  and  cardiovascular  system  should  tolerate  higher  whole-body  SARs;  an 
upper  limit  of  4  W/kg  for  the  First  Level  Controlled  Operating  Mode  is  recommended. 
Individual  tolerance  to  raised  body  temperature  is  highly  variable,  however,  even  in 
healthy  people  (NIOSH,  1986),  therefore,  medical  supervision  is  required. ' 

The  limits  described  above  were  developed  assuming  temperatures  in  the  Magnetic 
Resonance  equipment  examination  room  of  less  than  24  °C,  a  relative  humidity  of  less  than 
60%,  and  minimal  air  flow;  in  addition,  the  patient  is  assumed  to  be  lightly  clothed. 
Calculations  can  be  used  to  derive  correction  factors  for  environments  which  restrict  heat 
loss.  It  can  be  estimated  that,  for  each  °C  of  ambient  temperature  above  24  °C,  First  Level 
Controlled  Operating  Mode  restrictions  on  whole-body  SAR  should  be  lowered  by  0.25 
W/kg.  Similarly,  for  each  10%  of  excess  relative  humidity,  the  limits  would  be  lowered  by 
0. 1  W/kg.  For  the  Normal  Operating  Mode,  reductions  in  the  SAR  restrictions  should  be 
half  these  amounts,  except  that  it  is  not  necessary  to  reduce  the  limit  below  1  W/kg  for 
ambient  temperatures  less  than  30  °C.  Individual  responses  will,  however,  be  variable; 
increased  attention  at  the  appropriate  level  of  supervision  should  be  observed  when  using 
these  correction  factors. 

Some  regions  of  the  body,  such  as  the  head,  may  be  particularly  vulnerable  to  raised 
temperature.  The  developing  embryo  or  fetus  should  be  regarded  as  particularly  sensitive 
to  raised  temperatures;  however,  tissues  in  the  trunk  and  limbs  are  considered  less 
sensitive.  It  has  been  suggested  that  localized  temperatures  of  about  38  °C  in  the  head,  39 
"C  in  the  trunk  and  40  °C  in  the  extremities  are  unlikely  to  produce  adverse  effects. 
Simple  calculation  relating  localized  heating  in  the  eye  to  SAR  in  the  head  suggests  that 
exposure  resulting  in  3  W/kg  to  the  head  is  unlikely  to  raise  the  temperature  of  the  eye  by 
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more  than  1 .6  °C;  brain  temperatures  are  unlikely  to  rise  by  more  than  1  °C  under  these 
conditions.  It  seems  probable  that  higher  SARs  are  acceptable  for  the  trunk  (except 
during  pregnancy)  and  limbs  providing  that  whole-body  heat  loads  are  not  exceeded. 
However,  restrictions  on  local  tissue  SAR  which  will  appropriately  limit  the  temperature 
rise  in  these  tissues  are  not  well-defined;  there  is,  in  addition,  little  information  available  at 
present  to  justify  a  First  Level  ControUed  Operating  Mode  standards. 

The  rise  in  temperature  in  any  object  that  occurs  in  response  to  the  intermittent  absorption 
of  radio  fi-equency  energy  fi-om,  e.g.  a  pulsed  radio  fi-equency  source  can  be  equated  with 
the  SAR  value  averaged  over  the  50%  thermal  equilibrium  time,  the  time  taken  for  the 
temperature  increase  at  the  center  of  a  heated  region  to  rise  by  50%  of  the  maximum  value 
after  the  application  of  a  heat  source.  The  thermal  equilibrium  time  for  the  body  is  not 
known  accurately,  but  can  be  estimated  as  around  15  min  to  30  min,  whereas  smaller 
masses  of  tissues,  such  as  the  eye,  have  a  thermal  equilibrium  time  of  about  5  min. 

The  control  of  whole-body  SAR  alone  in  the  case  of  body-coil  imaging,  as  prescribed  by 
requirement  51.103  d)  1)  offers  a  simplification  of  the  provisions  for  control  of  the  SAR  in 
Magnetic  Resonance  systems.  The  simplification  enables  the  safe  use  of  a  relatively  large 
amount  of  radio  fi-equency  power  for  this  type  of  operation.  The  restriction  of  this 
simplified  operation  to  cases  where  the  transmit  coil  is  of  suflBcient  size  and  homogeneity 
(as  can  be  assumed  to  be  the  case  for  body  coils)  is  chosen  in  order  to  limit  the  ratio 
between  whole-body  average  SAR  on  the  one  hand  and  head  SAR  or  local  tissue  SAR  on 
the  other.  A  similar  line  of  argument  has  been  used  in  requirement  5 1 .  103  d)  2),  where 
the  control  of  the  system  of  the  head  SAR  alone  is  permitted  in  cases  where  a  head  coil  is 
used  as  the  transmit  coil. 


APPENDIX  B 
FDA  GUIDANCE 


Guidance  For  Magnetic  Resonance  Diagnostic  Devices  -  Criteria  For  Significant 
Risk  Investigations 


This  document  is  intended  to  provide  guidance  in  the  preparation  of  a  regulatory 
submission.  It  does  not  bind  the  FDA  or  the  regulated  industry  in  any  manner. 

Computed  Imaging  Devices  Branch 
Division  of  Reproductive,  Abdominal, 
Ear,  Nose,  Throat  and  Radiological  Devices 
OfiBce  of  Device  Evaluation 

Document  Issued  on  September  29,  1997 

While  this  guidance  document  represents  a  final  document,  comments  and 
suggestions  may  be  submitted  at  any  time  for  Agency  consideration  by 
writing  Robert  A.  Phillips,  Ph.D.,  Chief,  Computed  Imaging  Devices  Branch, 
Center  for  Devices  and  Radiological  Health,  HFZ-470,  9200  Corporate  Blvd., 
Rockville,  MD  20850.  For  questions  regarding  the  use  or  mterpretation  of 
this  guidance,  contact  Loren  A.  Zaremba,  Ph.D.,  at  301-594-1212. 

U.S.  Department  of  Health  and  Human  Services 
Food  and  Drug  Administration 
Center  for  Devices  and  Radiological  Health 


Magnetic  Resonance  Diagnostic  Devices 

Criteria  for  Significant  Risk  Investigations 

Patient  studies  utilizing  magnetic  resonance  diagnostic  devices  which  are 
conducted  under  any  one  of  the  following  operating  conditions  are 
considered  significant  risk  investigations,  and  require  approval  of  an 


116 


117 


investigational  device  exemption  (IDE)  by  the  Food  and  Drug  Administration 
(FDA)  Center  for  Devices  and  Radiological  Health  (CDRH): 

1 .  Main  static  magnetic  field  greater  than  4  tesla; 

2.  Specific  absorption  rate  (SAR)  greater  than 

a.  4  W/kg  averaged  over  the  whole  body  for  any  period  of  1 5  minutes, 
or 

b.  3  W/kg  averaged  over  the  head  for  any  period  of  10  minutes;  or 

c.  8  W/kg  in  any  gram  of  tissue  in  the  head  or  torso,  or  12  W/kg  in 
any  gram  of  tissue  in  the  extremities,  for  any  period  of  5 
minutes; 

3.  Time  rate  of  change  of  gradient  fields  (dB/dt)  sufficient  to  produce 
severe  discomfort  or  painfijl  nerve  stimulation;  or 

4.  Peak  unweighted  sound  pressure  level  greater  than  140  dB  or  A-weighted 
r.m.s.  sound  pressure  level  greater  than  99  dBA  with  hearing 
protection  in  place. 

These  criteria  apply  only  to  device  operating  conditions.  Other  aspects  of 
the  study  may  involve  significant  risks  and  the  study  may  therefore  require 
IDE  approval. 


APPENDIX  C 
PHANTOM  PROTOCOL  -  THERMOCOUPLE 


Protocol:  PHANTOM 


PATIFNT  POSl  HON 

SCANNING  RANGE 

Patient  Entry:              Head  First 

Field  of  View: 

24  cm 

Patient  Position:          Supine 

Scan  Thickness: 

5.0  mm 

Axial/Sag  Landmark: 

InterScan  Spacing: 

1.5  mm 

Coil  Type:                  Head 

Start  Location:  (I/S) 

135 

Scan  Plane:                 Axial 

End  Location:  (I/S) 

S36.5 

No.  Of  Scan  Locations: 

12 

IMAGING  PARAMETERS 

FOV  Center: 

Image  Mode:              2D 

(L/R): 

0 

Pulse  Sequence:          Spin  Echo 

(P/A): 

0 

SCAN  LIMING 

Number  of  Echoes: 

1 

IE: 

17 

rE2: 

153 

Rep  Time  (TR): 

500 

SCAN  SET-UP 

Prescan  Options: 

Autoshim 
Auto  CF:  Water 

Receive  bandwidth  (+/-)      1 6. 00 

kHz 


ACQUISITION  TIME 

Acq  Matrix:  freq  256 

Acl  Matrix:  phase  128 

Frequency  Direction:  A/P 
Phase  FOV:  24  cm 

Imaging  Time:  30  NEX 

48:25 
Contrast:  No 
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APPENDIX  D 
PHANTOM  PROTOCOL  -  FIBER  OPTIC 


ProtocoUS:  GEORGE  HARDER  TEST 

2.  Series:  Head, Axial,2D,  Spin  Echo,None 

Scan  1 :  HEAD,Axial,2D,Spin  Echo,  niemp,None 


PATIENT  POsrnoN 

SCANNING  RANGE 

Patient  Entry: 

Head  First 

Field  of  View: 

24  cm 

Patient  Position: 

Supine 

Scan  Thickness: 

5.0  mm 

Axial/Sag  Landmark: 

InterScan  Spacing: 

1.5  mm 

Coil  Type: 

Head 

Start  Location:  (I/S) 

122.5 

Scan  Plane: 

Axial 

End  Location:  (I/S) 

S22.5 

No.  Of  Scan  Locations: 

8 

IMAGING  PARAMETERS 

FOV  Center: 

Image  Mode: 

2D 

(L/R): 

0 

Pulse  Sequence: 

spin  Echo 

(P/A): 

0 

Imaging  Option: 

None 

Enter  PSD  Filename: 

ACQUISniON  TEVfF 

Acq  Matrix: 

freq  256 

SCAN  TIMING 

Acl  Matrix: 

phase  128 

Number  of  Echoes: 

1 

Phase  FOV: 

24  cm 

TE: 

30 

Frequency  Direction:A/P 

Rep  Time  (TR): 

500 

Imaging  Time:  1 

NEX1:12 

Contrast: 

No 

SCAN  SET-UP 

Prescan  Options: 

Autoshim 

Auto  CF: 

Water 

Receive  bandwidth  (+/-) 

16.00  kHz 
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APPENDEX  E 
HEAT  STRESS  INDEX 


The  focus  of  this  dissertation  was  to  derive  a  dosimetry  method  applicable 
to  the  human  head  and  the  phantom  head.  One  method  of  dosimetry  researched 
that  was  applicable  only  to  the  human  head  was  Heat  Stress  Index  (HSI). 
Industrial  Hygjenists  utilize  a  HSI  to  monitor  a  worker's  condition.  Since  the 
results  of  RF  energy  deposition  was  similar  to  that  of  energy  deposition  resuhing 
from  work,  a  HSI  was  originally  considered  as  an  option  for  RF  dosimetry.  This 
method  would  provide  a  qualitative  dose  reference  relying  on  an  individual's  own 
metabolism.  Physiological  parameters  such  as  heart  rate,  breathing  rate, 
temperature  and  perspiration  would  be  the  indicators  to  determine  an  individual's 
dose  limits.  Originally  this  method  seemed  plausible,  however,  reliability  would  be 
diflScult  to  ensure  for  all  patients.  This  appendix  is  included  to  provide  information 
about  heat  stress  indicators  and  demonstrate  the  reasoning  behind  their  unusability. 

The  body's  temperature-regulating  center  lies  in  the  hypothalamus.  The 
anterior  portion  of  the  hypothalamus  responds  to  its  own  increase  in  temperature 
and  to  nerve  impulses  from  receptors  in  the  skin.  The  hypothalamus  contains  a 
group  of  neurons  in  the  anterior  portion  referred  to  as  the  preoptic  area.  Increased 
blood  temperature  increases  the  rate  of  neuron  impulses  while  decreased  blood 
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temperature  decreases  the  rate  of  neuron  impulses  (Tortora  1981).  The  blood 
circulates  and  temperatures  can  fluctuate  between  35.5  °C  to  40  °C  with  it 
remaining  mostly  at  37  T.  During  MRI,  as  RF  energy  heats  local  tissue,  neuron 
impulses  increase.  The  hypothalamus  responds  by  triggering  increased  blood  flow 
to  transport  heat  to  the  skin  and  then  trigger  sweating  if  needed.  The  heat 
transported  to  the  skin  and  any  energy  deposition  in  the  skin  can  be  released  to  the 
environment  through  conduction,  convection,  radiation  and  evaporation  of  sweat 
(Adair  1992). 

The  most  common  indicators  of  heat  stress  used  by  Industrial  Hygienists 
are  sweat  rate,  skin  temperature,  heart  rate  and  core  temperature  (Patty's  1978). 
Sweat  rate  can  be  an  indicator  over  the  entire  range  of  heat  stress.  Its  rise  parallels 
a  core  temperature  rise  and  an  increase  m  heart  rate.  Sweat  rates  are  not 
instantaneous  but  are  a  time-weighted  average  of  heat  stress.  The  rate  of 
evaporation  is  the  difference  between  the  vapor  pressure  of  the  wetted  skm  and  the 
partial  pressure  of  water  vapor  in  the  ambient  air.  The  vapor  pressure  of  the 
wetted  skin  can  also  be  affected  by  the  skin  temperature.  Increased  skin 
temperature  increases  vapor  pressure  and  encourages  evaporation.  The 
endothermic  action  of  evaporation  then  reduces  skin  temperatures.  The  body's 
first  and  most  prevalent  response  in  maintaining  a  thermal  balance  without  a  rise  in 
core  temperature  was  to  increase  the  heart  rate.  This  increased  blood  flow 
channels  more  blood  past  heated  tissue  and  transports  the  heat  to  the  body  surface. 


122 
Heart  rate  has  the  dual  function  of  responding  to  increased  work  of  muscles  to  re- 
oxygenate  them  and  to  remove  heat. 

During  a  MRI  of  the  head,  there  are  two  areas  of  concern.  First,  increased 
heating  to  the  head  may  excite  the  hypothalamus  thereby  signaling  the  body  to 
respond  to  a  heat  stress  that  may  not  be  present  throughout  the  entire  body.  An 
increase  in  heart  rate  could  generate  an  unnecessary  increase  of  heat  loss  (Polk 
1986).  Second,  RF  signals  may  be  of  such  a  frequency  and/or  intensity  as  to 
generate  or  suppress  neuron  impulses.  This  may  also  affect  the  hypothalamus  into 
misinterpreting  signals.  Gradient  magnetic  fields  are  known  to  generate  false 
impulses  in  the  optic  nerve  manifesting  in  magnetophosphenes.  However,  there 
was  no  information  addressing  whether  RF  signals  can  generate  or  suppress  the 
neuron  impulses  to  the  hypothalamus. 

The  HSI  combines  the  heat-exchange  components  of  radiation  and 
convection  with  the  metabolic  heat  generation  of  the  body.  The  resultant 
expression  was  a  function  of  stress  in  terms  of  the  required  sweat  evaporation. 
The  physiological  strain  of  an  individual  can  then  be  expressed  as  a  percentage  of 
the  required  sweat  evaporation  to  the  maximum  evaporation  possible  by  an 
individual.  Stated  algebraically: 


HSI    =    100 


metabolic  heat  ±  radiation  ±  convection 


maximum  evaporation 
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A  HSI  greater  than  100  indicates  the  body  stores  heat  and  a  HSI  less  than 
100  indicates  the  body  loses  heat. 

Convection 

Convection  was  a  fonction  of  the  difference  in  temperature  between  the 
skin  and  the  surrounding  air  with  an  additional  effect  from  the  movement  of  air 
over  the  skin  surface. 

C=1.0V"*(TA-TS) 

where  C  =  convection  (kcal/hr) 
V  =  air  speed  (m/min) 
TA  =  air  temperature  (°C) 
TS  =  skin  surface  temperature  (°C) 

Convection  increased  when  the  air  temperature  (TA)  was  greater  than  the 
skin  surface  temperature  (TS)  and  the  individual  stores  heat.  The  exponent,  0.6, 
applies  to  forced  convection  over  vertical  cylinders  which  was  empirically  derived 
from  a  simplistic  geometric  configuration  of  a  man.  The  coefficient  1.0  was  a 
fijnction  of  the  surface  area  for  an  "average"  man  (1.8  m^).  A  non-average 
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individual  would  adjust  this  coefficient  as  a  ratio  of  their  estimated  skin  area  to  1 .8 
m^.  An  exact  estimate  was  not  required. 

Radiation 

Radiation  was  a  function  of  the  difference  in  temperature  between  the  skin 
and  the  mean  radiant  temperature  (MRT)  of  the  solid  in  contact  with  the  skin. 
Tests  have  shown  that  radiative  heat  exchange  was  a  function  of  the  fourth  power 
of  absolute  temperature  (Patty '  s  1 978).    Within  the  small  temperature  range  of  the 
equipment  that  a  patient  would  be  in  contact  with  during  an  MRI  a  first-order 
approximation  was  sufficiently  accurate  for  estimating  radiation. 

R  =11.3  (MRT -TS) 

where  R  =  radiation  (kcal/hr) 

MRT  =  mean  radiant  temperature  of  the  solid  surroundings  (°C) 
TS  =  skin  surface  temperature  (°C) 

The  coefficient  1 1.3  was  an  empirically  derived  value.  From  just  radiation 
and  convection,  the  body  loses  heat  wherever  air  touches  the  skin  and  a  cooler 
item  touches  the  skin.  Radiation  was  reduced  when  the  patient  was  placed  on  a 
cushion.  The  thermal  properties  ofthe  cushion  insulate  the  body.  Additionally, 
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since  the  MRI  chamber  was  kept  cool,  the  patient  may  also  be  provided  a  blanket. 

This  prevents  convection  and,  after  the  inside  of  the  blanket  has  reached 

equilibrium  with  the  skin  temperature,  also  prevents  radiation.  The  body's  natural 

heat  loss  avenues  are  now  restricted  to  evaporation. 

Evaporation 

Evaporation  was  a  function  of  the  difference  between  the  vapor  pressure  of 
perspiration  on  the  skin  (vapor  pressure  of  water  at  skin  temperature)  and  the 
partial  pressure  of  water  in  the  surrounding  air.  There  was  an  additional  effect  of 
air  speed  over  the  skin  which  can  reduce  the  partial  pressure  of  water  in  the 
immediate  air  and  increase  evaporation.  Evaporative  heat  loss  in  humid 
environments  was  reduced  as  the  capacity  of  the  ambient  air  to  accept  additional 
moisture  decreases.  In  a  MRI  chamber,  humidity  was  kept  low  thereby  increasing 
the  difference  in  the  vapor  pressure  of  perspiration  and  the  partial  pressure  of 
water. 

E„^  =  2  V'  (VPS  -  PPA) 

where  E^  =  maximum  evaporative  capacity  (kcal/he) 
V  =  air  speed  (ra/min) 
VPS  =  vapor  pressure  of  water  on  the  skin  (mm  Hg) 
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PPA  =  partial  pressure  of  water  in  air  (mm  Hg) 

In  dryer  environments,  such  as  the  MRI  chamber,  the  maximum  heat 
exchange  from  evaporation  becomes  limited  by  the  amount  of  perspiration  that 
can  be  produced  by  the  individual.  Maximum  sweating  was  sUghtly  above  one  liter 
per  hour  which  was  equivalent  to  an  evaporative  heat  loss  of  580  kcal/hr.  A  MRI 
patient  will  not  reach  this  limit  since  a  large  portion  of  skin  area  was  in  contact 
with  the  cloths  and  the  carriage.  These  pores  are  unable  to  remove  heat  from  the 
body. 

Example  information  concerning  heat  loss  from  patients  was  not  available, 
but  information  was  available  from  laborers.  For  a  laborer  and  under  favorable 
conditions,  25  percent  of  metabolic  heat  was  lost  through  convection  to  cooler  air, 
50  percent  by  radiation  to  cooler  surfaces  and  25  percent  by  evaporation.  In 
comparison,  a  patient  undergoing  an  MRI  has  a  larger  portion  of  their  body  in 
contact  v^th  the  carriage  which  would  account  for  a  large  percentage  of  heat  loss 
through  radiation.  If  the  patient  were  on  a  cushion  there  would  be  a  reduced  heat 
loss  through  radiation.  Additionally,  the  MRI  area  maintains  a  relative  humidity 
below  60%  which  becomes  favorable  for  heat  transfer  to  the  environment  through 
convection.  As  long  as  the  patient  has  the  ability  to  sweat,  heat  transfer  away  from 
a  patient  undergoing  an  MRI  examination  would  be  accounted  for  by  a  higher  heat 
loss  from  evaporation. 
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Body  temperatures  range  from  35.5  to  40  T  and  are  the  result  of  circadian 
variation,  vigorous  exercise,  variations  in  ambient  temperature,  sequelae  of  food 
intake,  age  factors,  cyclical  variation  in  females  and  emotional  factors.  In  the 
steady  state  the  body  obeys  the  first  law  of  thermodynamics  where  the  heat 
produced  in  the  body  was  balanced  by  the  heat  lost  to  the  environment. 
Quantitatively,  up  to  25  kcal/hr  of  energy  deposition  may  not  be  noticed  by  an 
individual.  For  reference,  greater  than  80  kcal/hr  represents  the  usual  upper  limit 
from  exercise  and  work.  At  double  this  value  (160  kcal/hr)  there  was  a  50%  risk 
of  collapse  and  triple  this  value  (240  kcal/hr)  becomes  intolerable.  Evaporation  of 
water  can  remove  0.58  kcal  per  gram.  The  body's  maximum  sweat  rate  was 
approximately  1  liter  per  hour  or  580  kcal  per  hour.  The  SAR  Umit  for  RF 
exposure  to  the  head  is  currently  at  3.2  W/kg.  This  equates  to  13.4  kcal/hr.  This 
was  below  the  25  kcal/hr  threshold  for  an  individual  to  notice  the  heat,  however, 
the  heat  may  be  deposited  in  a  small  region  such  as  the  vitreous  humor. 

Measurement  of  skin  temperature,  heart  rate  and  core  temperature  in  high 
magnetic  fields  was  possible.  However  using  HSI  for  dosimetry  would  be  difficult. 
A  baseline  of  responses  to  the  deposition  of  RF  energy  and  subsequent  increases  in 
thermal  loads  would  need  to  be  established  for  each  patient.  Each  individual 
responds  differently  as  a  fimction  of  their  age  and  conditioning.  Establishing  one 
universal  baseline  to  apply  to  all  patients  would  vary  as  a  fimction  of  the  following 
variables: 


128 
1  As  heat  was  deposited  within  the  body,  blood  flow  increases  to  transport 

the  heat  to  the  skin.  This  response  was  a  fiinction  of  an  individual's  cardiovascular 
fitness  and  was  unique  to  each  individual.  It  was  this  uniqueness  that  makes 
developing  a  uniform  standard  using  heart  rate  for  identifying  heat  stress  difficult. 
A  given  heart  rate  signifies  approximately  the  same  level  of  strain  regardless  of  the 
degree  of  an  individual's  fitness,  but  a  person  with  high  cardiovascular  fitness  will 
achieve  a  given  heart  rate  at  a  higher  level  of  thermal  stress  than  an  unfit  person. 
A  corollary  to  cardiovascular  fitness  was  the  effect  an  illness  imposes.  An  illness 
weakens  autonomic  functions.  This  would  prevent  developing  a  baseline  thermal 
stress  indicator  for  an  individual  and  using  it  for  subsequent  MRI  sessions. 

2.  Another  variable  of  heat  tolerance  was  obesity.  Obese  or  stocky 
individuals  have  a  lower  surface  area  to  weight  ratio.  An  individual  can  become 
acclimatized  over  time,  but  in  general  obese  individuals  are  at  a  greater  risk  of 
succumbing  to  heat  strain.  This  factor  was  slightly  mitigated  by  the  bore  size  of 
the  MRI  unit.  Some  obese  of  stocky  individuals  are  physically  limited  from  access 
to  higher  field  MRI  units. 

3.  An  older  individual  was  at  a  distinct  disadvantage  compared  Avith  younger 
individuals.  An  older  individual  compensates  for  increased  heat  less  effectively 
than  do  younger  individuals,  as  indicated  by  higher  core  temperature  and 
peripheral  blood  flow.  This  occurs  because  there  was  a  delay  in  the  onset  of 
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sweating  and  a  lower  sweat  rate  capacity  in  older  individuals  resulting  in  greater 
heat  storage  and  longer  time  for  heat  recovery.  Biological  eflfects  and  hazards  are 
a  function  of  the  dose  rate  and  not  the  total  dose.  Delays  in  reacting  to  a  thermal 
insult  would  cause  the  most  damage.  Acute  increases  in  core  temperatures  would 
not  trigger  heat  stress  indicators  to  provide  sufficient  warnings  of  danger.  Again, 
illnesses  would  further  alter  any  established  baseline  responses. 

4.         Drugs  can  adversely  affect  a  patient's  tolerance  from  thermal  insults.  Dogs 
have  displayed  a  greater  susceptibility  to  microwave  heating  after  administration  of 
pentobarbital  sodium,  morphine  sulfate  or  chlorpromazine,  indicating  that 
mechanisms  of  heat  loss  may  be  compromised  by  treatment  with  certain  drugs 
(Polk  1986).  RF  exposure  produced  greater  increases  in  core  temperature  than 
occurred  in  untreated  control  ammals  indicating  an  interference  with  the  body's 
natural  tendency  to  reduce  excess  heat.  The  variety  of  medications  administered  to 
mdividuals  would  need  to  be  tested  in  order  to  determine  their  effect  on 
thermoregulation.  Again,  individual  responses  to  medications  would  make  this 
difficult. 

While  the  concept  of  a  heat  stress  index  as  a  MRI  dosimeter  appears 
feasible,  the  practical  application  makes  its  effectiveness  tentative  (Shellock  1992). 
Without  extensive  testing  on  each  individual  to  determine  their  response  to  thermal 
stress,  quantifying  the  amount  of  energy  deposition  would  not  be  possible. 
However,  monitoring  heart  rate  and  skin  temperature  could  be  used  to  ensure 
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patients  remain  comfortable  and  are  not  in  the  early  stages  of  thermal  stress.  This 

would  become  valid  at  the  higher  field  strength  magnets  where  more  energy 

deposition  was  expected.  While  a  better,  more  quantitative  method  of  SAR  was 

required,  HSI  would  still  provide  a  good  secondary  indicator. 


APPENDIX  F 
EXPERIMENTAL  RESULTS 
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APPENDIX  G 
ARTIFACT  TUNING 


Ghost  Tuning.  ^  t  ^     . 

5.  Log  onto  the  InstaScan  Computer. 

6.  Remove  Ghost  images. 

a.  Go  to  the  ANMR  console. 

1 .  Select  a  New  Exam . 

2.  Select  Site  Protocol  List. 

3.  Select  Protocol  16.  Clinical  Heads. 

4.  Select  View  Protocols. 

5.  Select  Next. 

6.  Select  4.  Series  EPI  Sequences 

7.  Select  Scan  1:  H£AD^ial^D4(esearch,ismemp,None. 

8.  Select  Scan  Ops. 

9.  Select  Auto  Prescan. 

NOTE:  There  is  a  possibility  of  getting  a  message  "Auto  Prescan 
incomplete..."  In  this  case,  select  Manual  Prescan.  Once  an 
image  is  seen  on  the  display  screen,  select  Backup. 

b.  Go  to  the  InstaScan  computer. 

L         Select  Open  Database  and  select  desired  database  or  create  a  new 
database. 

2.  Select  Options  |  Save  Raw  Data. 

3.  Select  Scan  Setup. 

NOTE:  The  InstaScan  computer  will  direct  the  operator  to  begin 
the  scan  on  the  ANMR  console. 

c.  Go  to  the  ANMR  console. 
1.         Select  Scan. 

NOTE:  The  EPI  technique  creates  a  single  image  in  less  than  a 
second.  The  image  will  be  displayed  on  the  InstaScan  computer. 

d.  Go  to  the  InstaScan  computer. 

1.  Select  Utility  |  Autotuning 

2.  Tag  one  of  the  square  boxes  and  move  it  over  a  ghost  image.  Do 
not  place  the  box  on  an  actual  (desired)  image. 

3.  Repeat  with  the  second  box. 
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4.         Repeat  step  2.d.2.  through  2.d.4.  as  needed  to  remove  the  ghost 
images. 

NOTE:  If  the  desired  diffusion  technique  is  applicable  to  only  a 
single  slice  direction  then  removing  ghosts  from  the  other  two 
perpendicular  directions  is  not  required.  Proceed  to  step  3, 
otherwise  continue. 
e.  Go  to  the  ANMR  console. 

1.  Select  CANCEL. 

2.  Select  RxAhd  Exam. 

3.  Select  CANCEL. 

4.  Select  RxAhd  Series. 

NOTE:  This  methods  retains  the  positioning  information  of  the 
phantom  head.  If  this  is  not  used  then  the  phantom  head  will 
require  a  new  landmark  position. 

5.  Select  Acq.  Time 

6.  Change  Frequency  Direction  to  R/L. 

7.  Select  Scan  Ops. 

8.  Select  Auto  Prescan. 

NOTE:  There  is  a  possibility  of  getting  a  message  "Auto  Prescan 
incomplete  ..."  In  this  case,  select  Manual  Prescan.  Once  an 
image  is  seen  on  the  display  screen,  select  Backup. 

f  Repeat  2.b.2.  through  2.d.5. 

g.  Go  to  the  ANMR  console. 

1.  Select  Patient  Position. 

2.  Select  Scan  Plane  =  Sagital. 

3.  Select  Scanning  Range. 

4.  Enter  for  Start  (L/R):  LO. 

5.  Enter  for  End  (L/R):  RO. 

6.  Select  Scan  Ops. 

7.  Select  Auto  Prescan. 

NOTE:  There  is  a  possibility  of  getting  a  message  "Auto  Prescan 
incomplete  ..."  In  this  case,  select  Manual  Prescan.  Once  an 
image  is  seen  on  the  display  screen,  select  Backup. 
h.  Repeat  2.b.2.  through  2.d.5. 


APPENDIX  H 
DIFFUSION  IMAGING 


8.      Go  to  the  ANMR  console. 

i.  If  ghost  tuning  has  just  been  performed,  do  the  following. 

1.  Select  CANCEL. 

2.  Select  RxAhd  Exam 

3.  Select  CANCEL. 

4.  Select  RxAhd  Series. 

NOTE:  This  methods  retains  the  positioning  information  of  the 
phantom  head.  If  this  is  not  used  then  the  phantom  head  will 
require  a  new  landmark  position, 
j.  Select  the  difiusion  technique. 

1 .  Select  List/Sel  Protocol. 

2.  Select  Protocol  45. 

3.  Select  View  Protocol. 

4.  Select  Series  1. 

5.  Select  Scan  Ops. 

6.  Select  Auto  Prescan. 

k.  Go  to  the  InstaScan  computer. 

1.         Select  Scan  Setup. 

NOTE:  The  InstaScan  computer  will  direct  the  operator  to  begin 
the  scan  on  the  ANMR  console. 
1.  Go  to  the  ANMR.  console. 
1.         Select  Scan. 

NOTE:  The  difiusion  technique  creates  seven  images  within  44 
seconds.  The  image  will  be  displayed  on  the  InstaScan  computer. 
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APPENDIX  I 
DIFFUSION  IMAGING  PROTOCOL 


Protocol  45:  GEORGE  HARDER  TEST 

1.  Series:  Head,Axial,2D,Spin  Echo,None 

Scan  1:  HEAD,Axial,2D,Research,isclifif3T,None 


PATTFNT  POSITION 

SCAN  SET-UP 

Patient  Entry: 

Head  First 

AutoCF:                  ,   v: 

Water 

Patient  Position: 

Supine 

Axial/Sag  Landmark: 

SCANNING  RANGE 

Coil  Type: 

Head 

Field  of  View: 

40  cm 

Scan  Plane: 

Axial 

Scan  Thickness: 

5  mm 

InterScan  Spacing: 

5mm 

IMAGING  PARAME  lERS 

Start  Location:  (I/S) 

0 

Image  Mode: 

2D 

End  Location:  (I/S) 

0 

Pulse  Sequence: 

Spin  Echo 

No.  Of  Scan  Locations: 

1 

Imaging  Option: 

None 

FOV  Center:    (L/R): 

0 

Enter  PSD  Filename: 

isdiff3T 

(P/A): 

0 

USER  CV 

ACQUISITION  TEVIF, 

Image/Slice 

1.000000 

Acq  Matrix: 

freq 

Enhance 

1.000000 

256 

#  of  Shots 

1.000000 

Acl  Matrix: 

phase 

Catch  and  Hold 

0.000000 

128 

#  of  diffusion  steps 

2.000000 

Phase  FOV: 

'/2  20 

diffusion  weight  (s/mm2) 

900.000000 

cm 

Frequency  Direction: 

A/P 

SCAN  TIMING 

Imaging  Time: 

INEX 

TE: 

Minimum 

0:44 

FuU 

Contrast: 

No 

Rep  Time  (TR): 

Other  4000 
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APPENDIX  J 
IMAGE  TRANSFER 


2.         Go  to  the  IJF3T  Sun  workstation, 
m.        Log  onto  the  UF3T  system. 
n.         Change  directory  to  the  desired  location  for  placing  the  images. 

Use  the  Unix  command  CD  {path}, 
o.         Enter  xferx  at  the  prompt.  This  program  allows  transfer  of  images 

from  the  ANMR  to  a  temporary  directory  for  fUrther  processing. 
p.         Under  List  Exa  select  the  desired  exam  from  the  ANMR  system, 
q.         Under  List  Ser  select  the  desired  series  from  the  selected  exam. 
r.  Under  List  Ima  select  the  desired  images  from  the  select  series. 

s.  Hit  Get  From  Remote  to  transfer  the  images. 

t.  Hit  Quit  to  exit  the  xferx  program. 

u.         Enter  signa2unc  {path  only}  at  the  prompt.  Ensure  that  the  entire 

path  and  only  the  path  is  entered  after  the  signa2unc  command. 

This  program  converts  the  images  to  a  format  that  can  be 

recognized  by  the  Unix  computer. 
V.         Repeat  the  above  until  all  the  images  have  been  converted, 
w.        Enter  xidbmmaint  at  the  prompt.  This  is  the  only  acceptable 

method  for  deleting  images  retrieved  from  the  ANMR  system. 
X.         Select  connect. 
y.         For  SYS:  select  VA3. 
z.         For  PAT:  select  patient. 
aa.        Delete  images 
bb.       Quit  xidbmmaint. 
cc.        Enter  XV  at  the  prompt.  This  program  displays  the  images  and 

provides  a  method  of  transferring  the  images  from  a  Unix 

workstation  to  a  PC. 
dd.        Select  the  desired  image, 
ee.        Save  the  image  in  *.PGM  (ASCII)  format. 
S.        Repeat  the  conversion  for  all  images. 
gg.        Select  Quit. 
hh.        The  images,  in  a  *.PGM  (portable  gray  map)  format  can  now  be 

transferred  to  disc, 
ii.         Log  off  the  UF3  T  system. 
3.         Go  to  PC  system. 
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a.  Access  an  image  using  a  drawing  program,  such  as  Paint  Shop  Pro, 
by  JASC,  Inc. 

b.  Save  the  image  in  a  *.BMP  (bitmap)  format. 


1 ;" 


APPENDIX  K 
DIFFUSION  IMAGE  CONVERSION  PROGRAM 


Begin  Program. 

Private  Declare  Function  GetFocus  Lib  "user32"  Q  As  Long 

Private  Declare  Function  GetPixel  Lib  "gdi32"  (ByVal  hdc  As  Long,  ByVal  X  As 

Long,  ByVal  Y  As  Long)  As  Long 
Private  Declare  Function  SetPixel  Lib  "gdi32"  (ByVal  hdc  As  Long,  ByVal  X  As 

Long,  ByVal  Y  As  Long,  ByVal  crColor  As  Long)  As  Long 
Private  Declare  Function  GetBkColor  Lib  "gdi32"  (ByVal  hdc  As  Long)  As  Long 
Private  Declare  Function  GetColorSpace  Lib  "gdi32"  (ByVal  hdc  As  Long)  As 
Long 

Private  Declare  Function  LoadBitmap  Lib  "user32"  Alias  "LoadBitmapA"  (ByVal 

hinstance  As  Long,  ByVal  IpBitmapName  As  String)  As  Long 
Private  Declare  Function  GetDC  Lib  "user32"  (ByVal  hwnd  As  Long)  As  Long 
Private  Declare  Function  EmptyClipboard  Lib  "user32"  ()  As  Long 
Private  Declare  Function  CloseClipboard  Lib  "user32"  ()  As  Long 
Private  Declare  Function  OpenClipboard  Lib  "user32"  ^ByVal  hwnd  As  Long)  As 

Long 
Private  Declare  Function  SetClipboardData  Lib  "user32"  Alias 

"SetClipboardDataA"  (ByVal  wFormat  As  Long,  ByVal  hMem  As  Long) 
As  Long 

Private  Declare  Function  SetClipboardViewer  Lib  "user32"  (ByVal  hwnd  As 

Long)As  Long 
Private  Declare  Function  SelectObject  Lib  "user32"  (ByVal  hdc  As  Long,  ByVal 

hObject  As  Long)  As  Long 
Private  Declare  Function  CreateCompatibleBitmap  Lib  "gdi32"  (ByVal  hdc  As 

Long,  ByVal  nWidth  As  Long,  ByVal  nHeight  As  Long)  As  Long 
Private  Declare  Function  CreateCompatibleDC  Lib  "gdi32"  (ByVal  hdc  As  Long) 

As  Long 

Dim  Myhwd  As  Long 
Dim  MyDb  As  Database 
Dim  MySet  As  Recordset 
Dim  MySetA  As  Recordset 
Dim  MySetB  As  Recordset 
Dim  MyRay 
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Dim  Whichlmage  As  Integer 

Private  Const  CFBITMAP  =  2 

Private  Const  CFDIB  =  8 

Private  Const  CFDff  =  5 

Private  Const  CFDSPBITMAP  =  &H82 

Private  Const  CFDSPMETAFILEPICT  =  &H83 

Private  Sub  FormLoadQ 

Set  MyDb  =  OpenDatabase(App.Path  +  "\diflf.mdb") 
For  i  =  0  To  27 

Picl(i).  Visible  =  True 
Picl(i).ScaIe  (210,  0)-(550,  385) 
Next 
End  Sub 


Public  Sub  LoadlmageO 

MousePointer  =11 

MySet.MoveFirst 

Do  While  Not  MySet.EOF 

Pic2.PSet  (MySetC'X"),  MySet("Y")),  MySet("dwColor") 
MySet.MoveNext 

Loop 

Pic2.  AutoRedraw  =  False 

MousePointer  =  0 
End  Sub 

Public  Sub  DecipherIniage(ByRef  TableName  As  String) 
MousePointer  =  1 1     hourglass 
Set  MySet  =  MyDb.OpenRecordset(TableName) 
Label6  =  TableName 
Picl(0).SetFocus 
hwd2&  =  GetFocusO 
Myhwd  =  GetDC(hwd2&) 
For  i  =  0  To  385 
Label  1  =i 
Gaugel.Value  =  i 
Forj  =  210To550 

dwColor&  =  GetPixel(Myhwd,  j,  i) 
MySet.AddNew 
MySetC'X")  =j 
MySetC'Y")  =  i 
MySet("dwColor")  =  dwColor& 
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MySet.Update 
DoEvents 
Next 
Next 

Loadlmage 
MySet.Close 
Duration  =  Timer  +  3 
Do  While  Timer  <  Duration 

DoEvents 
Loop 

MousePointer  =  0 
End  Sub 

Private  Sub  Menu21_ClickO 

On  Error  Resume  Next  * 

MySet.Close 

Set  MySet  =  MyDb.OpenRecordset("002  Image  1") 

On  Error  GoTo  0 

Loadlmage 
End  Sub 


Private  Sub  Menu22_ClickO 

On  Error  Resume  Next 

MySet.Close 

Set  MySet  =  MyDb.OpenRecordset("002  Image  2") 

On  Error  GoTo  0 

Loadlmage 
End  Sub 


Private  Sub  Menu23_ClickO 

On  Error  Resume  Next 

MySet.Close 

Set  MySet  =  MyDb.OpenRecordset("002  Image  3") 

On  Error  GoTo  0 

Loadlmage 
End  Sub 


Private  Sub  Menu24_Click() 
On  Error  Resume  Next 
MySet.Close 


-5 
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Set  MySet  =  MyDb.OpenRecordset("002  Image  4") 
On  Error  GoTo  0 
Loadlmage 
End  Sub 


Private  Sub  Menu25_ClickO 

On  Error  Resume  Next 

MySet.Close 

Set  MySet  =  MyDb.OpenRecordset("002  Image  5") 

On  Error  GoTo  0 

Loadlmage 
End  Sub 


Private  Sub  Menu26_ClickO 

On  Error  Resume  Next 

MySet.Close 

Set  MySet  =  MyDb.OpenRecordset("002  Image  6") 

On  Error  GoTo  0 

Loadlmage 
End  Sub 


Private  Sub  Menu27_ClickO 

On  Error  Resume  Next 

MySet.Close 

Set  MySet  =  MyDb.OpenRecordset("002  Image  7") 

On  Error  GoTo  0 

Loadlmage 
End  Sub 


Private  Sub  Menu41_ClickO 

On  Error  Resume  Next 

MySet.Close 

Set  MySet  =  MyDb.OpenRecordset("004  Image  1") 

On  Error  GoTo  0 

Loadlmage 
End  Sub 


Private  Sub  Menu42_ClickO 
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On  Error  Resume  Next 
MySet.Close 

Set  MySet  =  MyDb.OpenRecordset("004  Image  2") 
On  Error  GoTo  0 
Loadlmage 
End  Sub 


Private  Sub  Menu43_ClickO 

On  Error  Resume  Next 

MySet.Close 

Set  MySet  =  MyDb.OpenRecordset("004  Image  3") 

On  Error  GoTo  0 

Loadlmage 
End  Sub 

Private  Sub  Menu44_ClickO 

On  Error  Resume  Next 

MySet.Close 

Set  MySet  =  MyDb.OpenRecordset("004  Image  4") 

On  Error  GoTo  0 

Loadlmage 
End  Sub 


Private  Sub  Menu45_ClickO 

On  Error  Resume  Next 

MySet.Close 

Set  MySet  =  MyDb.OpenRecordset("004  Image  5") 

On  Error  GoTo  0 

Loadlmage 
End  Sub 


Private  Sub  Menu46_ClickO 

On  Error  Resume  Next 

MySet.Close 

Set  MySet  =  MyDb.OpenRecordset("004  Image  6") 

On  Error  GoTo  0 

Loadlmage 
End  Sub 
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Private  Sub  Menu47_ClickO 

On  Error  Resume  Next 

MySet.Close 

Set  MySet  =  MyDb.OpenRecordset("004  Image  7") 

On  Error  GoTo  0 

Loadlmage 
End  Sub 


Private  Sub  Menu51_Click() 

On  Error  Resume  Next 

MySet.Close 

Set  MySet  =  MyDb.OpenRecordset("005  Image  1") 

On  Error  GoTo  0 

Loadlmage 
End  Sub 


Private  Sub  Menu52_ClickO 

On  Error  Resume  Next 

MySet.Close 

Set  MySet  =  MyDb.OpenRecordset("005  Image  2") 

On  Error  GoTo  0 

Loadlmage 
End  Sub 


Private  Sub  Menu53_Click() 

On  Error  Resume  Next 

MySet.Close 

Set  MySet  =  MyDb.OpenRecordset("005  Image  3") 

On  Error  GoTo  0 

Loadlmage 
End  Sub 


Private  Sub  Menu54_ClickO 

On  Error  Resume  Next 

MySet.Close 

Set  MySet  =  MyDb.OpenRecordset("005  Image  4") 

On  Error  GoTo  0 

Loadlmage 
End  Sub 
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Private  Sub  MenuSSCIickQ 

On  Error  Resume  Next 

MySet.  Close 

Set  MySet  =  MyDb.OpenRecordset("005  Image  5") 

On  Error  GoTo  0 

Loadlmage 
End  Sub 


Private  Sub  Menu56_Click() 

On  Error  Resume  Next 

MySet.  Close 

Set  MySet  =  MyDb.OpenRecordset("005  Image  6") 

On  Error  GoTo  0 

Loadlmage 
End  Sub 


Private  Sub  Menu57_ClickO 

On  Error  Resume  Next 

MySet.Close 

Set  MySet  =  MyDb.OpenRecordset("005  Image  7") 

On  Error  GoTo  0 

Loadlmage 
End  Sub 


Private  Sub  MenuQuitClickQ 

End 
End  Sub 


Private  Sub  MenuShowAllClickQ 
On  Error  Resume  Next 
Label  1.  Visible  =  False 
Label2.  Visible  =  False 
LabeW.  Visible  =  False 
Labels.  Visible  =  False 
Label6.Visible  =  False 
Gauge  1.  Visible  =  False 
'use  red  to  contrast  images 
For  k  =  14  To  27 
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Labels  =  "Contrast  box  "  &  Format(k,  "#") 

Picl(k).Cls 

Picl(k).  AutoRedraw  =  True 

Fori  =  210  To  550 
Forj  =  0To385 

Picl(k).PSet  (i,  j),  RGB(255,  0,  0) 
Next 

Next 

Pic  l(k).  AutoRedraw  =  False 
Next 

'display  images 
For  i  =  0  To  27 

Select  Case  i 
Case  0:  TableName  =  "002  Image  1" 
Case  1:  TableName  =  "002  Image  2" 
Case  2:  TableName  =  "002  Image  3" 
Case  3:  TableName  =  "002  Image  4" 
Case  4:  TableName  =  "002  Image  5" 
Case  5:  TableName  =  "002  Image  6" 
Case  6:  TableName  =  "002  Image  7" 
Case  7:  TableName  =  "004  Image  1" 
Case  8:  TableName  -  "004  Image  2" 
Case  9:  TableName  =  "004  Image  3" 
Case  10:  TableName  =  "004  Image  4" 
Case  11:  TableName  =  "004  Image  5" 
Case  12:  TableName  =  "004  Image  6" 
Case  13:  TableName  =  "004  Image  7" 
Case  14:  TableName  =  "Image  1" 
Case  15:  TableName  =  "Image  2" 
Case  16:  TableName  =  "Image  3" 
Case  17:  TableName  =  "Image  4" 
Case  18:  TableName  =  "Image  5" 
Case  19:  TableName  =  "Image  6" 
Case  20:  TableName  =  "Image  7" 
Case  21:  TableName  =  "Image  1-45" 
Case  22:  TableName  =  "Image  2-45" 
Case  23:  TableName  =  "Image  3-45" 
Case  24:  TableName  =  "Image  4-45" 
Case  25:  TableName  -  "Image  5-45" 
Case  26:  TableName  =  "Image  6-45" 
Case  27:  TableName  =  "Image  7-45" 
End  Select 
Label3  =  "Drawing  "  +  TableName 
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DoEvents 

MySet.  Close 

Set  MySet  -  MyDb.OpenRecordset(TableName) 

MousePointer  =11 

Picl(i).Visible  =  True 

Picl(i).Cls 

Picl(i).AutoRedraw  =  True 

MySet.MoveFirst 

Do  While  Not  MySet.EOF 

Picl(i).PSet  (MySetC'X"),  MySet("Y")),  Abs(MySet("dwColor")) 
MySet.MoveNext 

Loop 

Picl(i).AutoRedraw  =  False 

MousePointer  =  0 

DoEvents 
Next 

Labels.  Visible  =  False 
On  Error  GoTo  0 
End  Sub 

Private  Sub  MenuStartClickQ 

Picl(0).Cls 

Picl(O).  AutoRedraw  =  True 
Pic2.Scale  (235,  0)-(600,  301) 
ForImageNr=  12  To  21 
Select  Case  ImageNr 
Case  1:  PictureName  =  "T002_l.bmp" 

TableName  =  "002  Image  1" 
Case  2:  PictureName  =  "T002_2.bmp" 

TableName  =  "002  Image  2" 
Case  3:  PictureName  =  "T002_3.bmp" 

TableName  =  "002  Image  3" 
Case  4:  PictureName  =  "T002_4.bmp" 

TableName  =  "002  Image  4" 
Case  5:  PictureName  =  "T002_5.bmp" 

TableName  =  "002  Image  5" 
Case  6:  PictureName  =  "T002_6.bmp" 

TableName  =  "002  Image  6" 
Case  7:  PictureName  =  "T002_7.bmp" 

TableName  =  "002  Image  7" 

Case  8:  PictureName  =  "T004_l.bmp" 

TableName  =  "004  Image  1" 
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Case  9:  PictureName  =  "T004_2.bmp" 

TableName  =  "004  Image  2" 
Case  10:  PictureName  =  "T004_3.bmp" 

TableName  =  "004  Image  3" 
Case  11.  PictureName  =  "T004_4.bmp" 

TableName  =  "004  Image  4" 
Case  12:  PictureName  =  "T004_5.bmp" 

TableName  =  "004  Image  5" 
Case  13:  PictureName  =  "T004_6.bmp" 

TableName  =  "004  Image  6" 
Case  14:  PictureName  =  "T004_7.bmp" 

TableName  =  "004  Image  7" 
Case  15:  PictureName  =  "T005_l.bmp" 

TableName  =  "005  Image  1" 
Case  16:  PictureName  =  "T005_2.bmp" 

TableName  =  "005  Image  2" 
Case  17:  PictureName  =  "T005_3.bmp" 

TableName  =  "005  Image  3" 
Case  18:  PictureName  =  "T005_4.bmp" 

TableName  =  "005  Image  4" 
Case  19:  PictureName  =  "T005_5.bmp" 

TableName  =  "005  Image  5" 
Case  20:  PictureName  =  "T005_6.bmp" 

TableName  -  "005  Image  6" 
Case  21:  PictureName  =  "T005_7.bmp" 
TableName  =  "005  Image  7" 
End  Select 

If  ImageNr  =  1  Or  ImageNr  =  2  Or  ImageNr  =  1 1  Or  ImageNr  =  12  Or 
ImageNr  =  13  Or  ImageNr  =  1 5  Or  ImageNr  =  19  Or  ImageNr  =  20  Or  ImageNr 
=  21  Then 

Label3  =  PictureName 

Picl(0).Picture  =  LoadPicture(App.Path  +  "\"  +  PictureName) 
Pic2.Cls 

Decipherlmage  (TableName) 
End  If 
Next 

'Analyze  Image  Sets 

ReDim  MyRay(210  To  550,  0  To  385,  0  To  2) 
On  Error  Resume  Next 
MySet.  Close 
For  m  =  2  To  7 
Whichlmage  =  m 
Select  Case  m 
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Case  1:  Set  MySetA  =  MyDb.OpenRecordset("005  Image  1") 
Set  MySetB  -  MyDb.OpenRecordset("004  Image  1") 
Set  MySet  =  MyDb.OpenRecordset("Image  1-45") 
Case  2:  Set  MySetA  =  MyDb.OpenRecordset("005  Image  2") 
Set  MySetB  =  MyDb.OpenRecordset("004  Image  2") 
Set  MySet  =  MyDb.OpenRecordset("  Image  2-45") 
Case  3:  Set  MySetA  =  MyDb.OpenRecordset("005  Image  3") 
Set  MySetB  =  MyDb.OpenRecordset("004  Image  3") 
Set  MySet  =  MyDb.OpenRecordset("Image  3-45") 
Case  4:  Set  MySetA  =  MyDb.OpenRecordset("005  Image  4") 
Set  MySetB  =  MyDb.OpenRecordset("004  Image  4") 
Set  MySet  =  MyDb.OpenRecordset("Image  4-45") 
Case  5:  Set  MySetA  =  MyDb.OpenRecordset("005  Image  5") 
Set  MySetB  -  MyDb.OpenRecordset("004  Image  5") 
Set  MySet  =  MyDb.OpenRecordset("Image  5-45") 
Case  6:  Set  MySetA  =  MyDb.OpenRecordset("005  Image  6") 
Set  MySetB  =  MyDb.OpenRecordset("004  Image  6") 
Set  MySet  =  MyDb.OpenRecordset("Image  6-45") 
Case  7:  Set  MySetA  =  MyDb.OpeiiRecordset("005  Image  7") 
Set  MySetB  =  MyDb.OpenRecordset("004  Image  7") 
Set  MySet  =  MyDb.OpenRecordset("Image  7-45") 
End  Select 
Analyzelmage 
MySet.Close 
Next 

MyDb.Close 
End  Sub 


Private  Sub  Picl_GotFocus(Index  As  Integer) 
Select  Case  Index 
Case  0:  TableName  =  "002  Image  1" 
Case  1:  TableName  =  "002  Image  2" 
Case  2:  TableName  =  "002  Image  3" 
Case  3:  TableName  =  "002  Image  4" 
Case  4:  TableName  =  "002  Image  5" 
Case  5:  TableName  =  "002  Image  6" 
Case  6:  TableName  =  "002  Image  7" 
Case  7:  TableName  =  "004  Image  1" 
Case  8:  TableName  =  "004  Image  2" 
Case  9:  TableName  =  "004  Image  3" 
Case  10;  TableName  =  "004  Image  4" 
Case  11:  TableName  =  "004  Image  5" 
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Case  12;  TableName  =  "004  Image  6" 
Case  13:  TableName  =  "004  Image  7" 
Case  14:  TableName  =  "005  Image  1" 
Case  15:  TableName  =  "005  Image  2" 
Case  16:  TableName  =  "005  Image  3" 
Case  17:  TableName  -  "005  Image  4" 
Case  18:  TableName  =  "005  Image  5" 
Case  19:  TableName  =  "005  Image  6" 
Case  20:  TableName  =  "005  Image  7" 

End  Select 

LabeI6  =  TableName 
End  Sub 


Private  Sub  Picl_MouseMove(Index  As  Integer,  Button  As  Integer,  Shift  As 
Integer,  X  As  Single,  Y  As  Single) 

LabeW  =  Format(X,  "###") 

Label5  =  Format(Y,  "###") 
End  Sub 


Public  Sub  AnalyzelmageO 
MousePointer  =11 
'load  a  set  of  images  into  the  array 
MySetA.MoveFirst 

Label3  =  "Loading  A  - "  &  Format(WhichImage,  "#") 
DoEvents 
Do  While  Not  MySetAEOF 

MyRay(MySetA("X"),  MySetA("Y"),  0)  =  MySetA("dwColor") 

MySetA.MoveNext 
Loop 

MySetAClose 
MySetB.MoveFirst 

Label3  =  "Loading  B  - "  &  Format(WhichImage,  "#") 
DoEvents 
Do  While  Not  MySetB.EOF 

MyRay(MySetB("X"),  MySetB("Y"),  1)  =  MySetB("dwColor") 

My  SetB  MoveNext 
Loop 

MySetB.  Close 
DoEvents 
Label3  =  "Subtracting  Images" 
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Gaugel.Min  =  210 
Gaugel.Max  =  550 
DoEvents 

Label6  =  MySet.Name 
For  nx  =  210  To  550 
Gaugel.Value  =  nx 
DoEvents 
For  ny  =  0  To  385 
MySet.AddNew 
MySetC'X")  =  nx 
MySetC'Y")  =  ny 

MySet("dwColor")  =  MyRay(nx,  ny,  0)  -  MyRay(nx,  ny,  1) 
MySet.Update 
Next 
Next 

MousePointer  =  0 
DoEvents 
End  Sub 
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